Hi everyone, As part of some work I was doing, I had to develop a method for efficiently evaluating integrals of multidimensional, piece-wise-continuous functions based on a priori knowledge of where the jumps would happen. The resulting code takes a Sympy expression and processes it into a SciPy numerical integration call. The interesting bit is that the resulting integration call exactly fits the discontinuity(s), so the integrator only has to deal with the smooth regions, rather than trying to progressively refine a polynomial fit to a step function. Accomplishing that relies heavily on symbolic math manipulations and the code generation tools in Sympy.
Anyway, I would like to package this up in a way that would be publicly useful, but I'm not sure where it fits. Sympy seemed a likely guess, but the SciPy dependence is problematic. Alternatively, the discontinuity-processing could be decoupled and used with any iterated integrator that supports manually specified points of discontinuity. If Sympy isn't a good fit for this, I would appreciate suggestions about other places that might be better. Thanks for the help, Nathan Woods -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at http://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/a9223f76-fc9c-4182-b9c3-ca4f7aa95db4%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
