Thanks I have already discover this miskate and fixed it.
Regards.
2008/8/10 Daniel Stiles <[EMAIL PROTECTED]>
> Hello,
>
> That error means that you are missing a quote somewhere in your program. I
> believe it is here:
>
> except importError:
> print "The Scipy library is not available to \
> test the solution to the transient diffusion equation
>
> You should have a " at the end of that last line
>
> Sincerely,
> Daniel
>
>
> On Sat, Aug 9, 2008 at 2:33 PM, franck kalala <[EMAIL PROTECTED]
> > wrote:
>
>> Thanks for the answer;
>>
>> I tryed the first example in the documentation.
>> I wanted just to write code for the first part as a python programme
>> but I have an error message, could any one correct this this programme for
>> me?
>>
>> ***************************************************************************************************************
>> nx=50
>> dx=1.
>> from fipy.grid1D import Grid1D
>> mesh=Grid1D(nx=nx, dx=dx)
>> from fipy.variables.cellVariables import CellVariable
>> phi=CellVariable(name="solution variable",mesh=mesh,value=0)
>> D=1.
>> valueLeft=1.
>> valueRight=0.
>> from fipy.boundaryConditions.fixedValue import FixedValue
>> BCs=(FixedValue(faces=mesh.getFacesRight(), value=valueLeft),
>> FixedValue(faces=mesh.getFacesLeft(), value=valueRight))
>> from fipy.terms.explicitDiffusionTerm import ExplicitDiffusionTerm
>> from fipy.terms.transientTerm import TransientTerm
>> eqX=TransientTerm()==ExplicitDiffusionTerm(coeff=D)
>> timeStepDuration=0.9*dx**2/(2*D)
>> steps=100
>> x=mesh.getCellCenters()[...,0]
>> t=timeStepDuration*steps
>> from fipy.tools.numerix import sqrt
>> phiAnalytical=CellVariable(name="analytical value", mesh=mesh)
>> try:
>> from scipy.special import erf
>> phiAnalytical.setValue(1-erf(x/(2*sqrt(D*t))))
>> except importError:
>> print "The Scipy library is not available to \
>> test the solution to the transient diffusion equation
>> if __name__ =='__main__':
>> from fipy import viewers
>> viewer=viewers.make(vars(phi, phiAnalytical), limits={'datamin':0,
>> 'datamax':1.})
>> viewer.plot()
>> for step in range(steps):
>> eqX.solve(var=phi,boundaryConditions=BCs,dt=timeStepDuration)
>> if __name__=='__main__':
>> viewer.plot()
>> print phi.allclose(phiAnalytical, atol=7e-4)
>> if __name__ == '__main__':
>> raw_input("Explicit transcient diffusion. Press <return> to
>> proceed...")
>>
>> ******************************************************************************************************************
>> here is the message error:
>>
>> There an error in your program:
>> EOL while scanning single-quoted string
>>
>> I dont understand this error message
>>
>> regards
>>
>> Franck.
>>
>>
>>
>>
>>
>>
>> 2008/8/9 Jonathan Guyer <[EMAIL PROTECTED]>
>>
>>
>>>
>>> On Aug 9, 2008, at 10:15 AM, franck kalala wrote:
>>>
>>> Thanks for the answer. Then how do I set this initial value
>>>>
>>>> phi(x)=0 for x<0
>>>> phi(x)=1 for x>0
>>>>
>>>
>>> In FiPy 1.2:
>>>
>>> >>> x = mesh.getCellCenters()[..., 0]
>>> >>> phi.setValue(0, where=x < 0)
>>> >>> phi.setValue(1, where=x > 0)
>>>
>>> You should also specify what happens at x == 0.
>>>
>>>
>>> If you're using trunk/ then the first line would be
>>>
>>> >>> x, = mesh.getCellCenters()
>>>
>>>
>>>
>>>
>>
>>
>> --
>> ***********************************************
>> ***********************************************
>> Franck Kalala Mutombo
>> [EMAIL PROTECTED]
>> African Institut for Mathematical Sciences, Muizenberg Cape Town, South
>> Africa
>> [EMAIL PROTECTED]
>> ************************************************
>> ************************************************
>>
>
>
>
> --
> -Daniel Stiles
>
--
***********************************************
***********************************************
Franck Kalala Mutombo
[EMAIL PROTECTED]
African Institut for Mathematical Sciences, Muizenberg Cape Town, South
Africa
[EMAIL PROTECTED]
************************************************
************************************************
