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 the

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 on

Re: Residuals with vector equations - patch

2014-01-16 Thread Guyer, Jonathan E. Dr.
Oops. I missed that I was on the hook for code review. I've merged and pushed to develop. Thanks for pointing this out. On Jan 16, 2014, at 2:44 PM, Charles Reid wrote: > Hi all, > > FYI, I just noticed this was also reported as an issue on Trac > (http://matforge.org/fipy/ticket/658) and so

Re: Computation of finite element integrals

2014-01-16 Thread Daniel Wheeler
On Wed, Jan 15, 2014 at 5:15 PM, Benjamin Hepp 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 elements psi_i, wher

Re: Convection coefficients for vector equations

2014-01-16 Thread Daniel Wheeler
On Thu, Jan 16, 2014 at 1:17 PM, Charles Reid 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 > (1,Nvariables,Nv

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 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 curious, why did

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 b

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 elemen

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 wrote: > >> The vector formulation is still rather slow (much much slower than with just one equation). Is that expected? > > > I'm not sure. We'll d

Re: Mass and stiffness matrix

2014-01-16 Thread Daniel Wheeler
On Fri, Jan 10, 2014 at 10:27 AM, Benjamin Hepp 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 in the Finite Volume me

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 sa