[SciPy-User] Parallelizing a for loop

Askey, Scott A Capt USAF AETC AFIT/ENY Scott.Askey at afit.edu
Sat Dec 26 14:09:40 EST 2009


Chuck 
For a "small" problem 10 elements (420 dof total) runs in 2.2 second when fsolve is is given a jacobian function the problem solves in 2.2 seconds jerval=3 and feval =23 , without the jacobian it takes 20+ seconds.  

Scott
-
 The reason I asked is that there are a couple of finite element packages
 available for python. Two of them are SfePy <http://sfepy.org/> and
 FiPy<http://www.ctcms.nist.gov/fipy/>.
 I don't know if either is directly applicable to your problem. As to fsolve,
 are you directly computing the Jacobian or do you let the routine compute it
 numerically?

Chuck





More information about the SciPy-User mailing list