Dear FiPy developpers and users,
I'm trying to solve the following system of PDE and Non-Liner Algebraic
Equations (NLAE):
the PDEs read:
\frac{\partial{\phi}}{\partial t} = \nabla ( D(y) \nabla y)
where both phi and y are arrays of cellVariable and the number of PDE
equals the length of the \phi array
and the set of NLAE
F(\phi, y)= 0 at every node in the domain
where y is the local root of F, i.e. \frac{\partial F}{\partial y} is a
square full matrix.
This local problem is non-linear and I use SciPy.optimize.fsolve as
solver.
My problem:
I don't know an explict way of writing y = G{\phi} (and its derivatives)
as suggested for example here:
http://www.ctcms.nist.gov/fipy/examples/phase/generated/examples.phase.simple.html?highlight=expansion%20source
Therefore I wrote a local solver for F(\phi, y) (+ numerical derivatives)
and then use a series of setValue() getValue() to interact with FiPy
variables and solve the PDE sequentially with sweeps....and the whole
process becomes very slow as you mentionned in a recent post:
http://comments.gmane.org/gmane.comp.python.fipy/2568
Note that the process is slow even if I use the local solver at each time
iteration and not at each sweep.
My question:
Do you have any suggestion for improving the performance of this
resolution? and avoiding the use of setValue getValue for this kind of
coupled problem?
Thanks a lot for your help !!!!!!
Bruno HUET
Lafarge Centre de Recherche
"This e-mail is confidential and may contain legally privileged information.
If you are not the intended recipient, you should not copy, distribute,
disclose or use the information it contains.
Please e-mail the sender immediately and delete this message from your system.
E-mails are susceptible to corruption, interception and unauthorised amendment;
we do not accept liability for any such changes, or for their consequences.
You should be aware, that the company may monitor your emails and their
content."
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
