Hi Andre, On Wed, Nov 17, 2010 at 8:29 PM, freevryheid <[email protected]> wrote: > Hi > > I’m on a mission to code an open source elastic multi-layer analysis > system to calculate stress, strain and displacement response in > layered systems to surface loads. As proof of concept I opted to use > sympy and coded a simple two layer system under a single circular load > [1]. The vision is to eventually solve a system with many layers and/ > or loads. Most of the solution time appears to be in simplifying the > functions before integration - this step is necessary to avoid > division by zero errors (for example). Will performance improve using > Cython? Any other suggestions?
I tried your script and it hangs for me in the first simplify function print. I just tried to print the expression and it takes forever. What is the structure of the expressions? I tried print A0, and it's many pages. Can this be simplified? Can some other method be used instead of the LU solve? It seems that the expression is very big, so even things like expand() take forever. What kind of simplification is necessary --- maybe it's just expand? If so, maybe things can be expanded before they grow big (either use a different method than LU, or expand things before you pass it in). Ondrej -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to [email protected]. To unsubscribe from this group, send email to [email protected]. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
