Hi there, I found a bug in Fipy's Pysparse class that raises an error when sweeping toward a solution with a vector PDE. Some notes on the problem are here:
http://charlesmartinreid.com/wiki/Category:Fipy#Pysparse_TypeError_Using_Vector_Equations In short, the problem results from trying to treat a square matrix (number of variables x number of gridpoints) as a vector. This can be easily resolved by adding a "ravel()" call, which turns the matrix into a vector. I have a patch file here: http://files.charlesmartinreid.com/pysparseMatrix.patch Charles
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
