Thread Previous • Date Previous • Date Next • Thread Next |
On 27/05/10 15:58, Anders Logg wrote:
Wouldn't it work to just to this: v = v + (h/2.0*vnorm)*dot(velocity, grad(v)) U = 0.5*(u0 + u) F = v*(u - u0)*dx + 0.5*k*(v*dot(velocity, grad(U))*dx + c*dot(grad(v), grad(U))*dx) - k*v*f*dx a = lhs(F) L = rhs(F) That's my understanding of SUPG, but perhaps I've missed something?
Looks like the same thing (I'd have to check). I prefer to expressing the method as weighting the residual, e.g. Galerkin + weighted residual.
The code I added can be made more compact using lhs/rhs, but I was having some trouble and initially wanted to make a bare bones implementation.
Garth
-- Anders _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@xxxxxxxxxxxxxxxxxxx Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
Thread Previous • Date Previous • Date Next • Thread Next |