Hi everyone, I have a question that is related to my other thread "Axisymmetric geometry", but the point is different.
I need to compute "\int dtheta sin(theta) cos(theta) n(r,theta)" where n is a fipy cellvariable that is solution to my equation. Right now, I use n( (x,y,z) ) to evaluate the solution at a point in space. Are there better solutions to this problem? I have thought of dropping interpolation and using cellcenter values by selecting cellcenters such that r0 < r < r+dr where r=sqrt(x**2+y**2+z**2), and summing with some weigth. Actually, my interrogation comes from the fact that plotting n(r,theta) against theta gives a very noisy result (png file attached). Here, n(r,theta) is actually n( (x,y,z) ) with corresponding values of r,theta because my data is cartesian 3D. Thanks for any suggestion, Pierre _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
