|
|
Re: volume integration
Posted:
Sep 9, 2011 3:16 PM
|
|
Thanks so much Bruno. You blew my mind.
One last question,
For every node, actually I need to evaluate the integral 4 times. Using the same notation, I have JA JB JM JN
The integration is triple integral ( (JAJM -JAJN -JBM +JBN))
So, if I change the last part to something like this
JAJM = dot(JA,JM,3); JAJM = mean(JAJM,2);
JAJN = dot(JA,JN,3); JAJN = mean(JAJN,2);
JBJM = dot(JB,JM,3); JBJM = mean(JBJM,2);
JBJN = dot(JB,JN,3); JBJN = mean(JBJN,2);
J_ALL=(JAJM-JAJN-JBJM+JBJN)
jac2 = J_ALL.*abs(D)/6;
it should be alright.
But is possible to write this in the following way
J_ALL=dot(JA-JB,JM-JN,3) J_ALL=mean(J_ALL,2)
thanks in advance.
|
|