Re: derivative boundary condition in 2D problem

On Tue, Jan 24, 2017 at 4:33 AM, Francisco Vega Reyes wrote: > How do I implement a boundary condition \partial T /\partial y = 0 at > facesTop (for instance)? In that case, I don't believe that you need to add any constraints as that is the natural boundary condition for a

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 of type > 'fipy.matrices.scipyMatrix._ScipyMeshMatrix', creates a zero-dimensional > ndarray, with no shape. This isn’t

derivative boundary condition in 2D problem

How do I implement a boundary condition \partial T /\partial y = 0 at facesTop (for instance)? It seems to me that doing var.faceGrad.constrain(((0,),(0,)), where=mesh.facesTop) (the initial example in the manual) I would be like doing \partial T /\partial x + \partial T /\partial y = 0...