Attachment. On Fri, Mar 01, 2013 at 05:47:54PM +0100, Pierre de Buyl wrote: > 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 ]
-- ----------------------------------------------------------- Pierre de Buyl Université libre de Bruxelles Physique des Systèmes Complexes et Mécanique Statistique Membre du conseil d'administration - ULB EuroSciPy 2013 co-chair http://www.euroscipy.org/ tél: +32 (0)2 650 57 92 web: http://homepages.ulb.ac.be/~pdebuyl/ -----------------------------------------------------------
<<attachment: n_R0_theta.png>>
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
