On Tue, Nov 24, 2009 at 11:51 AM, Roy Stogner wrote:
>
> On Tue, 24 Nov 2009, John Peterson wrote:
>
>> On Tue, Nov 24, 2009 at 11:10 AM, Roy Stogner
>> wrote:
>>>
>>> It looks like the cholesky back substitution never reuses b(i) outside
>>> of (or at later iterations in) the i loop; surely we c
On Tue, 24 Nov 2009, John Peterson wrote:
On Tue, Nov 24, 2009 at 11:10 AM, Roy Stogner wrote:
It looks like the cholesky back substitution never reuses b(i) outside
of (or at later iterations in) the i loop; surely we can replace that
with a single scalar on the stack. The LU back substitut
On Tue, Nov 24, 2009 at 11:10 AM, Roy Stogner wrote:
>
> On Mon, 23 Nov 2009, John Peterson wrote:
>
>> {lu,cholesky}_solve(rhs, soln);
>>
>> This seems fairly reasonable on the surface, but there is a gotcha in
>> that both routines modify the rhs *as well as* the soln vector
>> internally.
>
> T
On Mon, 23 Nov 2009, John Peterson wrote:
> {lu,cholesky}_solve(rhs, soln);
>
> This seems fairly reasonable on the surface, but there is a gotcha in
> that both routines modify the rhs *as well as* the soln vector
> internally.
That's surprising. Do we need to do this? On first glance:
It lo
2009/11/24 Martin Lüthi :
> Hi
>
> Here is another publication for which where was used:
>
> @Article{Luethi2009,
> author = {M. P. L{\"u}thi},
> title = {Transient response of idealized glaciers to climate
> variations},
> journal = {Journal of Glaciology},
> year =
Hi
Here is another publication for which where was used:
@Article{Luethi2009,
author = {M. P. L{\"u}thi},
title ={Transient response of idealized glaciers to climate
variations},
journal = {Journal of Glaciology},
year = 2009,
number = 193,
volume =
