On Dec 19, 2013, at 12:59 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote:
I tried this and the residuals decrease quite a bit. Would you expect
the residuals to approach zero at every time step? This seldom happens
with non-linear equations. The absolute magnitude of the residuals
isn't
.)
Thanks again for the follow-up links!
Charles
On Thu, Dec 12, 2013 at 7:05 AM, Guyer, Jonathan E. Dr.
jonathan.gu...@nist.gov wrote:
On Dec 11, 2013, at 5:00 PM, Charles Reid charlesre...@gmail.com wrote:
Answered my own question - it looks like some modules went missing between
On Dec 5, 2013, at 6:19 AM, Richard Gillilan r...@cornell.edu wrote:
I'm finally writing up my results. Can anyone suggest a good paper to cite
crediting FiPy? Also, is there specific text you all would like to
see in the acknowledgment section (grant number and agency supporting FiPy
and
On Nov 12, 2013, at 12:15 PM, Daniel Wheeler daniel.wheel...@gmail.com wrote:
On Mon, Nov 11, 2013 at 4:30 PM, Richard Gillilan r...@cornell.edu wrote:
I'm rather new to FiPY too, but my impression is that viewers in FiPy are
very basic ... really intended only for monitoring results during
On Oct 18, 2013, at 6:50 PM, James Snyder jbsny...@fanplastic.org wrote:
If you're still interested:
http://ediacaran.mech.northwestern.edu/~jsnyder/podisoplot_conductor_animation.mp4
http://ediacaran.mech.northwestern.edu/~jsnyder/podisoplot_insulator_animation.mp4
These are done with the
For a FaceVariable, you'll need to use
X = mesh.faceCenters[0]
Y = mesh.faceCenters[1]
to get the shape right.
X and Y (or x and y) can be used to set the value because MeshVariables can be
cast to numpy arrays.
In this case, I'd be inclined to write:
vel = X * Y * [[1],[0]] + (Y-X) * [[0],
my iPhone
On Oct 17, 2013, at 12:07 PM, Guyer, Jonathan E. Dr.
jonathan.gu...@nist.gov wrote:
On Oct 17, 2013, at 9:43 AM, Richard Gillilan r...@cornell.edu wrote:
Congratulations on being able to open up again!
Thanks!
___
fipy
On Sep 30, 2013, at 4:38 AM, Damian Kösters koesters.dam...@gmail.com wrote:
I try now to work with a Constraint object. This doesn't give me an error
message but I also do not see the rectangle in the plot.
x, y = mesh.cellCenters[0], mesh.cellCenters[1]
X, Y = mesh.faceCenters[0],
On Sep 27, 2013, at 5:12 AM, Damian Kösters koesters.dam...@gmail.com wrote:
these are the errors I ran into:
1)
Traceback (most recent call last):
File damian_mail.py, line 47, in module
phi.constrain(phi.value[bacteria5], where=bacteria2)
IndexError: unsupported iterator index
On Sep 27, 2013, at 8:26 AM, Jon Guyer gu...@nist.gov
wrote:
2) I included
phi.constrain(phi[bacteria5.value], where=bacteria2)
before
viewer = Viewer(vars=phi, datamin=0., datamax=2.)
and never got beyond this line. (It didnt give me an error within two
minutes…)
That needs to be
On Sep 25, 2013, at 4:59 AM, Damian Kösters koesters.dam...@gmail.com wrote:
I changed the x,y section as you recommended.
I then tried inside the if step10-clause:
phi[bacteria2.value] = phi[bacteria5.value]
phi.constraints[2].value = phi.value[bacteria5]
It
On Sep 24, 2013, at 3:22 PM, Damian Kösters koesters.dam...@gmail.com wrote:
* Why is my selection bacteria of type type 'numpy.ndarray', whereas the
facesTopLeft is class 'fipy.variables.binaryOperatorVariable.binOp' ?
Because of a bug: http://matforge.org/fipy/ticket/447
Use:
x, y =
On Sep 23, 2013, at 3:47 AM, Damian Kösters koesters.dam...@gmail.com wrote:
I have got a problem I can't solve using COMSOL+Livelink to Matlab.
http://www.comsol.com/community/forums/general/thread/39531/
I just need diffusion calculation and reactions with the substances.
Would you
On Aug 28, 2013, at 5:39 PM, John Assael iass...@gmail.com wrote:
Moreover, what would be a good place to start reading about the error metrics
in order to do the refinement for the mesh i described?
Here a couple of articles.
Jasak and Gosman. Automatic resolution control for the
On Aug 26, 2013, at 5:35 PM, John Assael iass...@gmail.com wrote:
For example my smallest dimension in the mesh is a circle with diameter 6μ
(E-06). (where in order to get good results I have to define cellsizes of
1E-07)
I am completely confused. Previously you said that your 1 cm x 2 cm
If you'd like to post .geo file and the .msh file to pastebin or gist, I'll
take a look.
On Aug 13, 2013, at 9:08 PM, John Assael iass...@gmail.com wrote:
I dont know why but I get an empty array []... I am on MacOSX.
On Tue, Aug 13, 2013 at 11:02 PM, Guyer, Jonathan E. Dr.
jonathan.gu
On Aug 13, 2013, at 9:29 AM, John Assael iass...@gmail.com wrote:
Hi,
I have a problem when I am saving my mesh from gmsh (.msh file), fipy does
not recognise the Physical Groups, that I have defined.
More specifically, although my .geo file contains entries like: Physical
Surface(100) =
continuous help!
On Tue, Aug 13, 2013 at 8:28 PM, Guyer, Jonathan E. Dr.
jonathan.gu...@nist.gov wrote:
On Aug 13, 2013, at 9:29 AM, John Assael iass...@gmail.com wrote:
Hi,
I have a problem when I am saving my mesh from gmsh (.msh file), fipy does
not recognise the Physical Groups
Maybe you already did? http://stackoverflow.com/questions/tagged/fipy
Any posting from http://stackoverflow.com/search?q=fipy will do, though.
On Aug 2, 2013, at 3:03 PM, Phil Miller mille...@illinois.edu wrote:
Yeah, I can. Link to a relevant question, and I'll mark it as such.
On Fri, Aug
On Jul 18, 2013, at 4:58 AM, Edwin Sze Lun Khoo edwin...@mit.edu wrote:
J_face.setValue(-1.0/2.0*(c.harmonicFaceValue-rho/2.0)*phi.faceGrad)
In this case, you are getting a boundary value that depends on the
phi.faceGrad.constrain()
J_cell.setValue(-1.0/2.0*(c-rho/2.0)*phi.grad)
whereas
On Jul 18, 2013, at 3:02 PM, Edwin Sze Lun Khoo edwin...@mit.edu wrote:
I ran the code using 3 different methods:
• Without [] and with printing cL.faceGrad[0][0]
• With [] and without printing cL.faceGrad[0][0]
• With the mesh._orientedFaceNormals method you proposed and
On Mar 8, 2013, at 9:25 AM, Jonathan Guyer wrote:
On Mar 7, 2013, at 2:53 PM, Adrian Jacobo wrote:
On a related comment, I'm able to compile your example code using the openmp
wrapper for gcc:
openmpicc trivialmpi.c -o trivialmpi
and it runs just fine using openmpirun.
This is the
Thanks for those links, Dan. The wiki discussion is quite old. In addition to
packages listed there, of which PyTrilinos is probably the simplest route to
multigrid, I discovered recently that the authors of PyAMG have put together a
FiPy solver class that wraps their adaptive multigrid
I'm at a workshop right now, and my webmail client is unpleasant, so I'll
reserve a more complete answer until later, but I can confirm your results.
I've run examples.phase.anisotropy under a couple of different configurations
on my 10.5.8 MacBook Pro.
At 100x100, PySparse runs 100 steps in
It's just been pointed out to me that the reference I sent was ambiguous. The
issue is May/June. I have no idea where month = {Mar} came from.
@article{,
author = {J. E. Guyer and D. Wheeler and J. A. Warren},
journal = {Computing in Science and Engineering},
title = {FiPy:
301 - 325 of 325 matches
Mail list logo