Re: Convection coefficients for vector equations

2014-01-16 Thread Charles Reid
Hi there, I've managed to resolve the issues I was having. It came down to a problem with the type of the convection coefficient object. I was able to create a convection coefficient matrix by creating a list of size (1,Nvariables,Nvariables), and passing that as the coefficient. (Note that the

Re: Mass and stiffness matrix

2014-01-16 Thread Daniel Wheeler
On Fri, Jan 10, 2014 at 10:27 AM, Benjamin Hepp benjamin.h...@bsse.ethz.ch wrote: Hello, I have a diffusion problem: u_t = D * u_{xx} Is there any way to extract the mass and stiffness matrices from FiPy after everything has been setup? I'm not sure what the mass and stiffness matrices are

Re: Taking Dot Product of faceGrad

2014-01-16 Thread Jane Hung
I'm using PySparse On Jan 15, 2014, at 10:27 PM, Guyer, Jonathan E. Dr. jonathan.gu...@nist.gov wrote: On Jan 14, 2014, at 2:05 PM, Jane Hung jyh...@mit.edu wrote: The vector formulation is still rather slow (much much slower than with just one equation). Is that expected? I'm not sure.

Residuals with vector equations - patch

2014-01-16 Thread Charles Reid
Hi, I noticed that if I am solving a vector equation, running sweep() instead of solve() raises an error from the pysparseMatrix class's multiplication method. I looked deeper into the error and found that the problem was that a square matrix and a vector, both with the same number of total

Re: Residuals with vector equations - patch

2014-01-16 Thread Charles Reid
Hi all, FYI, I just noticed this was also reported as an issue on Trac ( http://matforge.org/fipy/ticket/658) and solved at a better/more upstream location ( http://matforge.org/fipy/changeset/2f015d7d9acd1953462a93a52e03c18303c8e2a6/fipy). But it doesn't look like that's made it to the main git

Re: Vertex vs Cell-centered finite volume

2014-01-16 Thread Daniel Wheeler
On Mon, Jan 13, 2014 at 10:33 PM, boyfarr...@gmail.com boyfarr...@gmail.com wrote: Dear list, When reading through the FiPy docs, http://www.ctcms.nist.gov/fipy/documentation/numerical/discret.html I noticed that discussion of vertex centred (CV) and cell centred (CC) meshes. I'm just

Re: Convection coefficients for vector equations

2014-01-16 Thread Daniel Wheeler
On Thu, Jan 16, 2014 at 1:17 PM, Charles Reid charlesre...@gmail.com wrote: Hi there, I've managed to resolve the issues I was having. It came down to a problem with the type of the convection coefficient object. I was able to create a convection coefficient matrix by creating a list of size

Re: Computation of finite element integrals

2014-01-16 Thread Daniel Wheeler
On Wed, Jan 15, 2014 at 5:15 PM, Benjamin Hepp benjamin.h...@bsse.ethz.ch wrote: Hello, my question from a few days ago might have been a bit imprecise. I reformulate it and hope this helps: I could not find any details in the FiPy documentation so I assume it is using first order Lagrange

Re: Taking Dot Product of faceGrad

2014-01-16 Thread Guyer, Jonathan E. Dr.
My tests show the coupled variant is actually faster (with Trilinos). As mentioned in ticket 658, the PySparse LU solver uses a lot of memory for this problem (it always uses a lot, but it's particularly high here). It's generally quite fast, but if it needs a substantial fraction of the RAM

Re: Vertex vs Cell-centered finite volume

2014-01-16 Thread boyfarr...@gmail.com
Hello Daniel, Thanks for the insight, having a large project which is still light weight and transparent is a challenge so I understand that approach. My interest in the CV approach came about because my 1D FVM code fails when trying to solve a particular electrostatics problem. I coded up