On Fri, Aug 5, 2016 at 12:15 PM, Nils Becker
wrote:
>
> ok. it seems that you have not changed the iteration, only added print
> functions. so my code seems to be converging, yes?
Absolutely. Appears to be converging well.
> that's also what i thought. however testing the solution externally, it
hi daniel, thanks for looking into this!
> Hi Nils,
>
> It seems like the code is converging well as far as I can tell. For
> example, when I modify the sweep loop to be
>
>
> print
> print 'step: {0}'.format(tt)
> while res > 1e-9: #Solve.sweeps:
> res
Hi Nils,
It seems like the code is converging well as far as I can tell. For
example, when I modify the sweep loop to be
print
print 'step: {0}'.format(tt)
while res > 1e-9: #Solve.sweeps:
res = eq.sweep(var=rho, dt=Solve.step,
solver=Solve.customSolver)
>> I'm not sure if using ||lhs|| is a good way to normalize the residual.
>> Surely that could be very small. Is ||lhs - rhs|| getting smaller
>> without normalization?
if you don't see any obvious errors in building the solution in fipy, i
could continue poking around if i knew how to evaluate th
> I'm not sure if using ||lhs|| is a good way to normalize the residual.
> Surely that could be very small. Is ||lhs - rhs|| getting smaller
> without normalization?
ok, i was not very precise in my statement. actually i was doing this:
1. make a plot of lhs pointwise.
2. make a plot of lhs - rh
On Wed, Aug 3, 2016 at 9:27 AM, Nils Becker
wrote:
>> Could you possible post your code so I can play around with it? I
>> can't think offhand what could be the issue or exactly what the error
>> is that you're calculating, it seems like the normalized residual or
>> something.
>
> hi daniel, than
> Could you possible post your code so I can play around with it? I
> can't think offhand what could be the issue or exactly what the error
> is that you're calculating, it seems like the normalized residual or
> something.
hi daniel, thanks for your reply. the way i was calculating the residual
Hi Nils,
Could you possible post your code so I can play around with it? I
can't think offhand what could be the issue or exactly what the error
is that you're calculating, it seems like the normalized residual or
something. This is a linear equation so I'm not entirely sure why the
non-linear res
hi,
i scanned the list and the documentation but have not been able to make
this work properly yet:
i have a 2D pde of the form
\[
\partial_t \rho = \nabla \cdot (\rho \nabla \phi) + D \nabla^2 \rho +
\lambda \rho
\]
where phi is a potential field and lambda is a source field, both are
const
