Dear all,
I'd like to make a slight change to transient_ex1 such that half of the
domain has a low diffusivity and the other half of the domain has a high
diffusivity. I wonder how I can get x coordinates of a quadrature. I
tried putting
if (qface_points[qp](0) < 0.0) { // subdivice the domain w.r.t x
= 0.0
diff_coeff = 1.0; // high diffusivity
} else {
diff_coeff = 0.5; // low diffusivity
}
in the quadrature for loop, but it seemed not working. Any adivce will be
appreicated.
Best regards,
K. Lee.
------------------------------------------------------------------------------
Own the Future-Intel® Level Up Game Demo Contest 2013
Rise to greatness in Intel's independent game demo contest.
Compete for recognition, cash, and the chance to get your game
on Steam. $5K grand prize plus 10 genre and skill prizes.
Submit your demo by 6/6/13. http://p.sf.net/sfu/intel_levelupd2d
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
