Hello, I'd like to calculate the integral of (grad u)^2 over my domain, after a calculation.
Currently, I assign two new variables, and set up new equations for them: Jx-du/dx=0, Jy-du/dy=0. Then I just square and add the results of calculate_norm (Jx and Jy, L2). It seems to work. But this adds two extra field variables and equations to the calculation, multiplying the memory usage and slowing things down a good bit. Alternatively, is it possible to calculate the L2 norm of the field derivative or gradient directly, through a method I've missed? Or if there's nothing quite that convenient, is there a way of calculating the derivatives and putting them into a separate vector in the system, decoupled from the original calculation? Thanks, -Adam -- GPG fingerprint: D54D 1AEE B11C CE9B A02B C5DD 526F 01E8 564E E4B6 Engineering consulting with open source tools http://www.opennovation.com/
signature.asc
Description: This is a digitally signed message part
------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
