CVS log for JSOC/proj/sharp/apps/sw_functions.c
	Help

Request diff between arbitrary revisions

changed the numerical derivative. here is a description of the problem that i fixed:

In computeJz(), the newly malloc'd jz array has its interior
points filled with values but its boundary points are left
alone:

    for (i = 1; i <= nx-2; i++)
    {
        for (j = 1; j <= ny-2; j++)
        {
            // calculate jz at all points
            ...

In all subsequent functions the array jz is used whole, including
the uninitialized border.

I ran your test case in the debugger.  On n02, jz was 0 on the
border; on solar3, I got wild values

(gdb) p jz[0]
$16 = 1.37765386e+10
(gdb) p jz[1]
$17 = 1.35754556e-19
(gdb) p jz[2]
$18 = 1.35631939e-19

If the code relies on the border being 0, then it was pure luck
that it ran correctly on x86_64.