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

Dear Daniel, Jonathan, Thank you both very much for your helpful replies. Krishna and I can now access the 'L' discretisation matrix in Lx = b, as desired, using the code example at the bottom of this e-mail. We have some follow-up questions, though: 1) Applying numerix.array() to 'L', when

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

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 >

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

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

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

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