Re: [Octave forge] [bim package] bim1a_advection*

[JuanPi]

>  Instead of passing the vector field BETA directly one can pass a
>      piecewise linear conforming scalar function PHI as the last input.
>      In such case BETA = grad PHI is assumed.

Looking at the source code I see that actually it is BETA = diff PHI,
not grad PHI what is used (I see now why you were puzzled)
Another misleading documentation point, if I understand correctly.
Why not using the function gradient, instead of diff to compute these steps?

-----
“An article about computational result is advertising, not
scholarship. The actual scholarship is the full software environment,
code and data, that produced  the  result.”
- Buckheit and Donoho