#L88
with return self.distanceVar._cellInterfaceFlag * 1e+20. You
probably want to construct the metal equation in your script to avoid
editing these files directly. That will make c=0 on the interface
(phi=0).
v= \nabla c
I think we've already discussed the above condition.
--
Daniel Wheeler
://matplotlib.sourceforge.net/faq/installing_faq.html#backends
--
Daniel Wheeler
D=0.
--
Daniel Wheeler
Sorry.
On Fri, Mar 4, 2011 at 5:39 PM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
You can check and see that the example you are trying to run don't not
match your installed version of FIPy by importing the Grid2D instance
and looking at the arguments it takes.
You can check and see
=.
$ python -c 'import fipy; print fipy.__version__'
What do you get? Is it different from 2.2-dev?
$ python examples/levelSet/electroChem/simpleTrenchSystem.py
Does this now work?
--
Daniel Wheeler
is in the following form:
\frac{\partial b}{\partial t} = D \frac{\partila ^2 b}{\partial x^2} +
rb(1-b/k)
where r and k are constants.
That is a different equation from above. Do you have a problem solving
it or posing it in FiPy?
--
Daniel Wheeler
to be more recent than
the version of fipy you are using.
--
Daniel Wheeler
express that
in terms of boundary conditions?
Sorry. Can't really help with that. We need equations and boundary
conditions to make progress.
--
Daniel Wheeler
()
phi.updateOld()
c.updateOld()
phi.extendVariable(c.getFaceGrad())
dt = cflNumber * cellSize / c.getFaceGrad().max()
advectionEquation.solve(phi, dt=dt)
ceq.solve(var=c, dt=dt,boundaryConditions=BCs)
--
Daniel Wheeler
for installing http://matforge.org/fipy/downloader/download/release/16.
--
Daniel Wheeler
matplolib.__version__
?
You can try editing line 124 of
C:\Python25\lib\site-packages\fipy\viewers\matplotlibViewer\matplotlib2DGridViewer.py
from
self.image.clim(vmax=datamax, vmin=datamin)
to
self.image.set_clim(vmax=datamax, vmin=datamin)
--
Daniel Wheeler
() has only been tested in 1D
cases, where this doesn't matter. Again, Wheeler will be happy to look into
this when he gets back... won't you Wheeler?
I'll change the examples to use vector values.
--
Daniel Wheeler
.getCellCenters()
v1 = CellVariable(mesh=m1, value=x1 * y1)
m0 = Grid2D(nx=2, ny=2)
v0 = CellVariable(mesh=m0, value=v1(m0.getCellCenters()))
I haven't checked the above, so possible typos.
--
Daniel Wheeler
value (ie not on a mesh ) ?
X and Y can have arbitrary values. They can also be arbitrary array values.
--
Daniel Wheeler
?
What happens when you print oldphi((xc,yc))? Is it actually non zero
and is it the same length as phi? You should also be able to do
newphi.setValue(oldPhi) if they are the same length.
--
Daniel Wheeler
phase[((x - cx)**2 + (y - cy)**2 + (z - cz)**2) r**2]= 0
print pore volume fraction=,phase.getCellVolumeAverage()
-
--
Daniel Wheeler
to symbol.
Am I missing something here?
Thanks
Ram.
--
Daniel Wheeler
is working fine for me now if I run it from its
own directory. I am using version 2.1. I was going to try 2.1.1
yesterday before sending mail but the site was broken for downloading
:-)
Thank you very much for your help.
Lucia
--
Daniel Wheeler
that fipy identifies are found in
C:\Python26\Lib\site-packages\. I know that EPD comes with mayavi and can be
found in C:\Python26\Lib\site-packages\enthought\. How do I make fipy
recognize mayavi which is in a different folder?
Thanks
Ram.
--
Daniel Wheeler
to know. Don't worry about
the tar.gz if you have already installed with the the .exe.
Cheers
--
Daniel Wheeler
these complicated boundaries in fipy?
You can't have interior boundary conditions in FiPy. Cast the boundary
condition to a volume integral and apply a source term based on location or
use Benny's suggestion, which can more accurately capture the boundary.
--
Daniel Wheeler
with a variable coefficient? If so, the answer
is yes. The coefficient can be any vector field that changes in space
and time, but the term can not be posed to FiPy in the form you gave.
Hope the above helps.
--
Daniel Wheeler
for ConvectionTerm needs to be a vector. Try,
cv = CellVariable(mesh = m, rank=1)
__ConvectionTerm(coeff = cv)
rank=1 makes the variable a vector value.
--
Daniel Wheeler
.
Best,
Erik
On Oct 5, 2010, at 3:09 PM, Daniel Wheeler wrote:
Hi Erik, I've attached a script that seems to run faster. It took
57.60s to do 200 steps as opposed to your original script, which took
183.50. In order to get some of this speed up, I had to hack the FiPy
code a little. The hack
the
residual to the initial residual, rather than some arbitrary number.
On Sep 30, 2010, at 12:55 PM, Erik Sherwood wrote:
Thanks for taking a look at the code. I'm using FiPy 2.1.
Erik
On Sep 30, 2010, at 11:21 AM, Daniel Wheeler wrote:
Eric, Thanks for your interest in FiPy. I'll
# Update the plots
if __name__ == '__main__':
Uviewer.plot()
Vviewer.plot()
Wviewer.plot()
#
toc = clock()
print Finished in %.3f secs.%(toc-tic)
--
Daniel Wheeler
%(step)
[snip]
# Update the plots
if __name__ == '__main__':
Uviewer.plot()
Vviewer.plot()
Wviewer.plot()
--
Daniel Wheeler
free to test it if you have an appropriate platform. Details for
building the deb and testing are here
http://matforge.org/fipy/blog/FiPyDeb. Cheers.
--
Daniel Wheeler
equation Daniel talks about
would not also solve the problem better, however, that fortran code I refer
to works apparently great without such a coupled solver.
Benny
--
Daniel Wheeler
-discussion
--
Daniel Wheeler
version.
Regards,
Adya
From: fipy@nist.gov [f...@nist.gov] on behalf of Daniel Wheeler
[daniel.wheel...@gmail.com]
Sent: Wednesday, July 28, 2010 9:41 AM
To: Multiple recipients of list
Subject: Re: Problem Installing Fipy
On Wed, Jul 28, 2010 at 10
On Wed, Jul 28, 2010 at 10:25 AM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
* Try deleting everything under C:\Pythonxy and install Enthough's
That should be C:\Python26.
--
Daniel Wheeler
in the near future. It uses a fully coupled matrix for
all the equations rather than the segregated approach used thus far in
FiPy. I'll let you know when that becomes live. Cheers.
On Mon, Jul 19, 2010 at 5:30 AM, Benny Malengier
benny.maleng...@gmail.com wrote:
2010/7/6 Daniel Wheeler daniel.wheel
prefer Imposing a
pressure value via a boundary condition on the pressure correction equation,
this seems to work great if one sets eg a boundary at a specific pressure to
a fixed value.
Okay, if it works that is much cleaner.
--
Daniel Wheeler
is causing it.
Benny
2010/7/15 Jonathan Guyer gu...@nist.gov
On Jul 15, 2010, at 3:02 AM, Benny Malengier wrote:
It seems ndarray has a dot function now in numpy trunk, updated my numpy
and fipy crashes like ...
Thanks for the alert
--
Daniel Wheeler
--
Daniel Wheeler
any
insight into this problem or would try this test and report what they see.
For example, did I forgot to disable the ENABLE_CATASTROPHIC_MEMORY_GUSHER
option in trilinos (Default = True)?
Thanks and regards,
+jtg+
--
Daniel Wheeler
counldn't figure out how to extract the data from
trilinos's local matrices, but eventually we'd like to not be using
pysparse when using the trilinos solvers. Just explaining why there is
an issue with pysparse even when running against trilinos.
On Fri, Jun 4, 2010 at 9:55 AM, Daniel Wheeler
Maybe try without any preconditioner or with only Jacobi.
--
Daniel Wheeler
On May 25, 2010, at 5:04 PM, Jonathan Guyer gu...@nist.gov wrote:
On May 24, 2010, at 11:49 PM, jtg wrote:
Continuing with Dr. Guyer's anisotropy example, here are
our times for his 3 mesh sizes (same 10 steps
there are no artifacts at all.
Which example was this? I'll take a look. Thanks for the report.
--
Daniel Wheeler
partitioning. The parallel changes to the code are independent
of the partitioning scheme so this shouldn't require major changes
just to the Gmsh mesh classes. Hopefully will be part of the next
release, but no guarantees right now.
--
Daniel Wheeler
some time to do it.
--
Daniel Wheeler
Make that branches/version-2_1 rather than trunk/.
On Tue, May 18, 2010 at 5:28 PM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
It might be easier to do a fresh checkout of trunk rather than copying
files around.
What happens if you put a print statement here
http://matforge.org
...@gmail.com
To: fipy@nist.gov
Subject: Re: Problem in solving Poisson equation with 1D Cylindrical mesh
Make that branches/version-2_1 rather than trunk/.
On Tue, May 18, 2010 at 5:28 PM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
It might be easier to do a fresh checkout of trunk rather
would appreciate any help!
Regards
Arvind Ajoy
PhD Student, Electrical Engineering
Indian Institute of Technology Madras
python.exe has encountered a problem and needs to close. We are sorry
for the inconvenience.
--
Daniel Wheeler
to incorporate those conditions
and
join these two meshes together.
Any help would be greatly appreciated.
Respectfully yours,
Dariusz Kedzierski
--
Daniel Wheeler
need to use the logarithmic discr. so after some
hours
of
sleep
I'll think about the reason.
Hotmail: Trusted email with powerful SPAM protection. Sign up
now.
--
Daniel Wheeler
the logarithmic discr. so after some hours
of
sleep
I'll think about the reason.
Hotmail: Trusted email with powerful SPAM protection. Sign up now.
--
Daniel Wheeler
Your E-mail and More
, I haven't really
looked into it, but it never causes problems so not particularly
motivated to take action.
--
Daniel Wheeler
and for his suggestions
and timing results(!) -- a significant speed
increase seems achingly close: I can see
it now, just over the rainbow...
--
Daniel Wheeler
On Mon, May 10, 2010 at 1:34 PM, jtg hi.jo...@gmail.com wrote:
I nominate this test for inclusion in the FiPy manual,
as a quick way to verify that Trilinos has been built
properly.
Seconded!
--
Daniel Wheeler
On Wed, May 5, 2010 at 11:14 AM, Alexander Holbach a.holb...@gmx.net wrote:
Is this possible?
It's often easier to reformulate the boundary condition as a source
term, this is probably what you need to do.
--
Daniel Wheeler
can formulate
this as a source term. This sometimes requires reformulating the
boundary condition from an area integral to a volume integral over a
grid width.
--
Daniel Wheeler
.
--
Daniel Wheeler
On Thu, May 6, 2010 at 10:55 AM, Igor gozdawa.a...@googlemail.com wrote:
In the problem, it should be updated each time-step.
As long as time is a Variable and the coefficient is somehow based
on time, it should just work.
--
Daniel Wheeler
automatically update if I use the same approach?
With some luck it should just work.
--
Daniel Wheeler
. Get it now.
--
Daniel Wheeler
with radius (external) = 0.03 m,
so quite little.
Just build it explicitly into the expression for the source term.
There is no need for a 2D grid. If you must have depth then just use a
2D grid with 1 grid cell. It is basically the same as using a 1D grid.
Hope this helps.
--
Daniel Wheeler
the Schroedinger equation,
I need to solve an eigenvalue problem to find allowed energies. Is there
a way to tell fipy to search for the required eigenvalues?
Many thanks,
John Gamble
--
Daniel Wheeler
/ticket/291 to track this.
--
Daniel Wheeler
other packages already
present at my PC and used by fipy v2.0.2, which is installed globally?
Thank you for the answers.
Cheers,
Igor.
--
Daniel Wheeler
=timeStepDuration,solver=LinearLUSolver()
,boundaryConditions=bcs)
Try this with and without an explicit source term as well.
--
Daniel Wheeler
--
Daniel Wheeler
, Mar 19, 2010 at 12:00 PM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
If alpha and beta are zero, which they seem to be in the equation that
you wrote down. Then you just have the three terms, which is easy to
set up. The source term will just have phi.getGrad()[0] -
phi.getOld().getGrad
(mesh.getExteriorFace(), N)
If you need more information please let me know, and I appreciate any help you
can give me.
Thanks,
--
Michael Johnston
812-290-1578
Purdue University
Materials Science Engineering
mljoh...@purdue.edu
--
Daniel Wheeler
, but as I said, you don't need to worry about it now.
Cheers
--
Daniel Wheeler
On Mon, Feb 1, 2010 at 10:48 AM, Daniel Wheeler
daniel.wheel...@gmail.com wrote:
Hi Igor,
I made two changes to get it to work. This is the diff:
22,27c22,23
convTerm = PowerLawConvectionTerm((1,0))
diffCoeff = FaceVariable(mesh=mesh_xy, rank=2)
diffCoeff[0] = 5.
Wrong! Should
the derivatives and this is often preferable.
--
Daniel Wheeler
line option -Wno-long-double
error: command 'gcc' failed with exit status 1
Computer: iMac intel
OS: MacOSX v10.6.2
--
Daniel Wheeler
%2F
Cheers
--
Daniel Wheeler
, Biswajit Ray r...@purdue.edu wrote:
Hi,
Is there any easy way to implement periodic boundary condition using fipy.
Can anyone please discuss on it as i did not find anything related in
the manual.
Thanks
--
Daniel Wheeler
, which is
not that bad. I am going to work on some examples out of the manual, but if
they don’t work as well, then I would hope to be able to get it fixed.
Thanks,
Kevin Buckman
--
Daniel Wheeler
On Tue, Nov 3, 2009 at 10:39 AM, Jonathan Guyer gu...@nist.gov wrote:
var2D[mask.getValue()] = z
Even better.
--
Daniel Wheeler
[index] = value
IndexError: unsupported iterator index
I don't know how to fix it. Please give me a hint. Many thanks for your
help.
Kundar
--
Daniel Wheeler
this.
var2D.setValue(var3D((x2D, y2D, numerix.ones * max(z3D
--
Daniel Wheeler
a way to do this without a loop without reshaping var3D
and collapsing along the z-axis and then you'd have to assume that x2D
and y2D somehow correspond to x3D and y3D, which I'd rather not do.
--
Daniel Wheeler
().getCellCenters()[:2], order=1))
4. It seems this just expand var2D into 3D, like the same value along
Z-axial. Is my understanding right?
Yes, is that what you want to do? I wasn't clear on what you actually needed.
--
Daniel Wheeler
if this problem has sharp changes in the variables
being solved for. In the end you
may need a higher order term still, but this could at least form a
basis to improve on.
--
Daniel Wheeler
into 3D mesh:
var3D.setValue(var2D(var3D.getMesh().getCellCenters()[:2], order=1))
I didn't test the above, but that is how it is advertised to work.
--
Daniel Wheeler
as a result of e-mail transmission.
NOTICE REGARDING PRIVACY AND CONFIDENTIALITY
Knight Capital Group may, at its discretion, monitor and review the content
of all e-mail communications.
http://www.knight.com
--
Daniel Wheeler
image001.jpg
as the diffusion coefficient in the end layer
of cells and B * A for gammaAdsorption.
--
Daniel Wheeler
On Mon, Sep 28, 2009 at 3:39 PM, Kun-Dar Li kun...@umich.edu wrote:
Another question is: can I describe the convection term as
convTerm=(convCoeff*var.getArithmeticFaceValue()).getDivergence()
convCoeff is a vector. What's the difference if I define the convection
term
by myself without
and the
getDivergence operator doesn't know about boundary conditions. You can
of course include boundary conditions by modifying the term
appropriately.
See, How do I represent a convective source? in the FAQ in the manual.
Cheers
--
Daniel Wheeler
(a)
resultOfConvolution = CellVariable(mesh=mesh)
eqn = someTerms + resultOfConvolution
while True:
resultOfConvolution,setValue(convolution(var))
eqn.solve(var)
Cheers
Thanks,
-Ranjit
--
Daniel Wheeler
?
numerix.array() or [:] is the most canonical, numerix.array will work
for all types though, [:] will break for stuff without length.
Thanks for pointing out this weird behavior of Variable's.
--
Daniel Wheeler
creates a new
array as does a copy. Should getValue() be creating a copy if it
doesn't do an evaluation?
--
Daniel Wheeler
]
A.setValue(x)
print A.getGrad().getValue()
tsv=TSVViewer(vars=[A, A.getGrad()])
tsv.plot('getgrad_test.txt')
--
Daniel Wheeler
both the mesh and solver or simply algebraic
multigrid and multilevel proconditioners, which can be though of as
simply a solution to the linear ststem?
--
Daniel Wheeler
--
Daniel Wheeler
being quite strongly coupled (due to the
matrix M(x)), so iterating over each equation does not work very well.
This is my experience hence the creation of the prototype coupled
equation class.
Cheers
--
Daniel Wheeler
algebraic
multigrid and multilevel proconditioners, which can be though of as
simply a solution to the linear ststem?
--
Daniel Wheeler
, but shouldn't be for periodic
boundaries.
Thank
you so much for any help.
--
Daniel Wheeler
it in a different way.
Cheers
--
Daniel Wheeler
On Tue, Jun 9, 2009 at 2:42 PM, William Gathrightgathri...@gmail.com wrote:
This seems very similar to a problem Daniel Wheeler reported in
February 2009. I was unable to find a resolution to that problem in
the mailing list.
I still haven't resolved this problem, but I haven't spent a lot
- http://matforge.org/fipy/ticket/191 Windows testing bug
Cheers
--
Daniel Wheeler
this anisotropy problem? Please give
me some hints. Thanks in advance. Have a good weekend.
Kundar
--
Daniel Wheeler
() (used instead of solve()) returns the unnormalized L2
residual. You can use this to determine number of sweeps and time
steps etc.
Cheers
--
Daniel Wheeler
to proceed...)
Could you help?
Thanks a lot
--
David Wende
home +972-8-9353488
work +972-2-5886116
mobile +972-54-234-6479
--
Daniel Wheeler
.
--
Daniel Wheeler
to be like that. I'd think that much better
performance could be had with array manipulations.
--
David Wende
home +972-8-9353488
work +972-2-5886116
mobile +972-54-234-6479
--
Daniel Wheeler
lists or arrays or the exact data types or something. I wrote it that way a
couple of years ago in r2055 and r2056, but my svn comments aren't enough of
a clue as to why it has to be like that. I'd think that much better
performance could be had with array manipulations.
--
Daniel Wheeler
401 - 500 of 628 matches
Mail list logo