Dear list, I see how integration can be included by converting to a sum (as explained here, http://thread.gmane.org/gmane.comp.python.fipy/2231/focus=2232),
x = mesh.getCellVolumes() area = (y * x).sum() However, I want to integrate a cell variable where the upper limit of integration is the current position. See here for the actual expression, http://i.imgur.com/4tNUt.jpg This can be solved really easily by integrating the numpy arrays over the specific range. For example, using scipy.quad function, import scipy import numpy numpy.exp( -scipy.quad( lambda x: phi(x), 0, x)[0] ) How can I perform this integral using Fipy? It is possible to incorporate this function call to scipy? Best wishes, Dan
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
