Time varying boundary conditions

Hi All, I have a problem where I would want to specify fixed value boundary conditions for times t=t0. The way I've implemented this is: a.constrain(a_0,(mesh.exteriorFaces)) while t = t0: a.faceGrad.constrain(0.,mesh.exteriorFaces) eq.solve(dt=timeStepDuration) t+=dt

Conda installation not working

Hi, For some reason my fipy conda installation stopped running. I've removed it and installed it again but I keep getting the same error: >>> import fipy Traceback (most recent call last): File "", line 1, in File

Steppers

Hi, Is there any example on how to use the different steppers? Any documentation on possible advantages/disadvantages of each one? Thanks, Adrian. ___ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY:

Re: Trilinos error in Anaconda installation in OsX

Jon Here is the error message I get when I try to run the Anaconda installation of FiPy using Trilinos on OsX. I'm using El Capitan: Traceback (most recent call last): File "CH_SqMesh.py", line 8, in import fipy as fp File

Re: Parallel speedup

) and worked for a couple of people I was at a workshop with. > If you've done this since then, I'd appreciate knowing what doesn't work. > > - Jon > >> On Feb 6, 2017, at 2:09 PM, Adrian Jacobo <ajac...@mail.rockefeller.edu> >> wrote: >> >> Hi, >> &g

Parallel speedup

Hi, I'm trying to speed up the Cahn-Hilliard example by running it in parallel but the running time is always the same regardless of how many processors I use. I'm using FiPy in Linux installed from Anaconda (I tried the same installation in OsX but trillinos doesn't work) following the

Scalar variable

Hi, I'm trying to solve a system where an equation for a field is coupled to an equation for a scalar. But I'm having problems writing this into FiPy. To illustrate the point here I show you a simplified version of the equations I'm trying to solve:

Adaptive time step

Hi, Is there a simple way to implement an adaptive time step in Fipy? More specifically, is there a way that, when calling the solve function store the time advanced fields in temporary variables?, so that the evolution can be recalculated using a different time step (and then comparing

Re: Saving in parallel

the available memory in the system and crashes. Cheers, Adrian. On Mon, Mar 18, 2013 at 1:22 PM, Adrian Jacobo ajac...@mail.rockefeller.edu mailto:ajac...@mail.rockefeller.edu wrote: Hi, I'm having problems when saving results of a calculation on a Gmsh grid when running

Saving in parallel

Hi, I'm having problems when saving results of a calculation on a Gmsh grid when running in parallel. I've attached an example of the problem, based on the diffusion/circle.py example. The code test_circle_parallel.py runs the diffusion equation on a circle and then saves the result to a

Re: Matplotlib Viewer

Sorry! I saw that was a bug with Numpy and I didn't pay that you where having problems with the viewer too. Have a look to my previous email about numpy 1.7 bug !!! ___ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST

Matplotlib Viewer

memory, and crashes). Thank you so much for all your help! sorry to bother you with so many questions! Adrian. On Tue, Mar 12, 2013 at 10:16 AM, Adrian Jacobo ajac...@mail.rockefeller.edu mailto:ajac...@mail.rockefeller.edu wrote: Hi, Thanks for your help. I gave up trying to make

Re: Trilinos installation

As you suspected, I've installed openmpi using macports, and this might be the source of the problem. Here I send you the result of 'otool -L' for the mpi4py and PyTrilinos libraries. otool -L

Re: Trilinos installation

MPI_C_COMPILER and MPI_CXX_COMPILER options when calling cmake, but it tells me is not using those values) Adrian. On Tue, Mar 5, 2013 at 6:38 PM, Adrian Jacobo ajac...@mail.rockefeller.edu mailto:ajac...@mail.rockefeller.edu wrote: Hi everyone, I know this might not be the right place

Trilinos installation

Hi everyone, I know this might not be the right place to ask, but may be someone can give me a hint. I'm trying to install Trilinos in OSX Mountain Lion, to use it with fipy. I tried following the instructions in fipy's manual, but the trilinos installation gives an error, saying that it

Re: Point Source

Thanks! That's what I was looking for, and it seems to work fine (although having to use an arbitrary large value to define the implicit source is not very elegant, but it does what I need to) Adrian. Hi Adrian, I'm not sure if this is what you're going for, but perhaps you could use an

Re: Point Source

Hi Pierre, Thank you for your suggestion, I now implemented a quick and dirty version fixing the value of the field at the required points, but what you suggest seems like a more elegant solution, I will try it later and see how it compares. Adrian. Hello, In general, defining a singular

Solving two coupled equations with nonlinear terms

Hi, I'm trying to solve the following equations: And I'm defining the system in this way: mesh = fp.Grid3D(dx=dx, dy=dy,dz=dz, nx=nx, ny=ny,nz=nz) C = fp.CellVariable(name = "Calcium",mesh = mesh, value = 0.) B

Re: Solving two coupled equations with nonlinear terms

Not for me. If you're seeing something else, then it may come down to either boundary conditions or parameter values. We probably need to see a minimal script that exhibits what you're seeing. I send you attached a minimal script that shows one of the problems I'm having. Now I'm solving

Re: Integrating field

Hi, Thank you for your answer. I'm solving an equation in 3d and I need to calculate the integral of the field. What is the most convenient/fast way of doing this? int = var.cellVolumeAverage * var.mesh.numberOfCells Just to clarify, as it might be useful for somebody else, this

Integrating field

Hi all, I'm solving an equation in 3d and I need to calculate the integral of the field. What is the most convenient/fast way of doing this? A related question: how can I access individual points of the field? Right now, to access the point C(x,y,z) I'm doing C(((x,),(y,),(z,))), where C is