Re: Question on accessing internal matrices of the system being solved

2017-02-09 Thread Guyer, Jonathan E. Dr. (Fed)
he system. >> >> Can FiPy return the system matrices and RHS vector by discretising only the >> spatial domain? >> >> Best Regards, >> >> Krishna >> >> From: fipy-boun...@nist.gov on behalf of >> Daniel Wheeler >> S

RE: Question on accessing internal matrices of the system being solved

2017-02-09 Thread Campbell, Ian
Krishna -Original Message- From: fipy-boun...@nist.gov [mailto:fipy-boun...@nist.gov] On Behalf Of Guyer, Jonathan E. Dr. (Fed) Sent: 26 January 2017 14:42 To: FIPY Subject: Re: Question on accessing internal matrices of the system being solved I would imagine you could subclass Tran

Re: Question on accessing internal matrices of the system being solved

2017-01-26 Thread Guyer, Jonathan E. Dr. (Fed)
3 PM > To: Multiple recipients of list > Subject: Re: Question on accessing internal matrices of the system being > solved > > Hi Ian, > > Sorry for the slow response. > > On Thu, Jan 12, 2017 at 12:20 PM, Campbell, Ian > wrote: > > > > 1) App

Re: Question on accessing internal matrices of the system being solved

2017-01-25 Thread Gopalakrishnan, Krishnakumar
, January 24, 2017 5:00:23 PM To: Multiple recipients of list Subject: Re: Question on accessing internal matrices of the system being solved Hi Ian, Sorry for the slow response. On Thu, Jan 12, 2017 at 12:20 PM, Campbell, Ian wrote: > > 1) Applying numerix.array() to ‘L’, when ’L’ is o

Re: Question on accessing internal matrices of the system being solved

2017-01-24 Thread Daniel Wheeler
Hi Ian, Sorry for the slow response. On Thu, Jan 12, 2017 at 12:20 PM, Campbell, Ian wrote: > > 1) Applying numerix.array() to ‘L’, when ’L’ is of type > 'fipy.matrices.scipyMatrix._ScipyMeshMatrix', creates a zero-dimensional > ndarray, with no shape. This isn’t what we expected because L has d

RE: Question on accessing internal matrices of the system being solved

2017-01-12 Thread Campbell, Ian
Right) v, L, b = DiffusionTerm(D)._buildMatrix(phi, SparseMatrix) print(type(L)) print(L) print("\n") numerixed_L = numerix.array(L) print(type(numerixed_L)) print(numerixed_L.shape) print(numerixed_L.ndim) print(numerixed_L) -Original Message- From: fipy-boun...@nist.gov [mailto:f

Re: Question on accessing internal matrices of the system being solved

2017-01-12 Thread Guyer, Jonathan E. Dr. (Fed)
By "display", I mean show a color map. It's not quantitative, but it can be useful to understand the stencil being generated. > On Jan 12, 2017, at 10:53 AM, Guyer, Jonathan E. Dr. (Fed) > wrote: > > In addition to this, if you define an environment variable > FIPY_DISPLAY_MATRIX, then FiPy w

Re: Question on accessing internal matrices of the system being solved

2017-01-12 Thread Guyer, Jonathan E. Dr. (Fed)
In addition to this, if you define an environment variable FIPY_DISPLAY_MATRIX, then FiPy will display the matrix and RHS vector. If FIPY_DISPLAY_MATRIX contains "terms", then the contribution of each term will be displayed before the aggregate matrix and vector are displayed. If FIPY_DISPLAY_M

Re: Question on accessing internal matrices of the system being solved

2017-01-12 Thread Daniel Wheeler
Hi Krishna, Yes you can do this, see https://github.com/usnistgov/fipy/blob/develop/fipy/terms/term.py#L332 The following should work, equation.cacheMatrix() equation.solve(...) the_matrix = equation.matrix The cacheMatrix is require because otherwise the reference to the matrix is

Question on accessing internal matrices of the system being solved

2017-01-11 Thread Gopalakrishnan, Krishnakumar
Hi, I was wondering if there is a way to access the internal matrices of the system of equations formulated by Fipy , given the coefficient terms of the general conservation equation. Let's say I have a simple diffusion equation (eg. Heat equation). I can hand-code the differentiation matrix a