One rule of the community board is not to post a problem that you have not yet taken effort to solve.
homogen, of {homogen,bc,ic} is not defined above. but it still hangs if defined.
I think the problem you've posted isn't very clear (for example, what u = g2/g1 entails as far as PDE solving method). I don't see clearly what is is the PDE, IC, BC to be solved, so knowing 'the problem' would require (me) to work some while just to un-cover what is getting solved. What is L? G? S? V? indicator[?] (this is ok if you cite a book page, but without that...) Is this and example from a book like a homework or test question (for credit)? Why does S depend on integrating f[] which is namingly v1[] which is Identity[x]=x, but also not the same f[] as in the PDE representing heat (normally u[] is heat in other books, but u and f are reversed)?. ic = {f[x, T - l] == y0[x]} would mean f[...] = f[...] which contains no new information so is not "an equation", which books call the bc, unless you "happen to know Identity[x] = x" and remember that f is heat and u is the non-homogenous part. I also see that using w1=Sin is possibly obstructive.
"normally the calculations are not complicated at this point". Could you provide an example of that? If you are saying that is so, then it would be nice to see such a similar example that was not a problem.
It appears to me this is a non-homogenous pde problem involving heat equation as solved by heat kernel simplification, convolution, and maybe generalized function, all which are examples to find in a PDE book and possibly use with a duhamel's chapter combined. But the above has linear IC and insulated 0 BC, which may need steady + transient consideration, not typical of example problems of these kind. The piecewise is not a generalized function, so that is out. A series of width 5 (in this situation) is pretty wide and should take minutes to solve always. And because u[] contains a heat kernel like form doesn't necessitate one to use a fourier/guass method: actually that method uses the integral to solve (not NDSolve at all), and your problem has a heat kernel "like" phrase in u[] (which again, most books use u for heat and f for function, you've done the reverse).
(use u as heat, f is the non-homogenous part and given, 0 BC, IC is a simple function (just x))
In all I would then call it ... "Nonhomogeneous PDE using Eigenfunction Expansions”, type problem (or superimposition / fourier technique). But u[] (and thus f[]), is periodic and needs to be solved "within 0 to Pi domain", u[] may not be expressed in a form for it to work for all x and t ... and what method works does depend on the interval (finite, infinite, one-sided) - I have to think a little harder if it needs to be polar, i would know infinite intervals may work easier with fourier. I think DSolve "does not normally easily solve this problem type". For certain similar problems, in polar form (and this problem might use it), "we can transfer the non-homogeneity from the PDE to the BC", so it's not totally clear without some pre-solving what must be done (if we are forced to consider the bc and u as separate or if infact they could be combined).
At any rate, the ways to solve a non-homogenous PDE with BC and IC are many. And that is why(for readers) it would be good to have a more straightforward question. u[] as heat, f[] as a function (not swapped up), clear pde ic bc to be solved (no mystery substitutions and functions unless explicitly called out). Hard enough to do these if we see what we are trying to solve.
summary: i am used to seeing this same kind of problem being solved in multiple steps "so that mathematica can solve them (it can't otherwise)", i cannot say if version 12.1 should be able to do them quickly