[LORENE] compactified external domain

[Jerome Novak]

Dear All,

As you may know, in the compactified external domain, we are using a
flag (called dzpuis) that indicates the power of r by which one must
multiply the data stored in a Scalar to get the real (physical)
function.
e.g. if the data show a constant field equal to 1, but the dzpuis is set
to 2, then the physical field is 1/r^2 ...

I have just changed the behavior of this flag, when calculating
derivatives or 1/r type operators on a Scalar (the old Cmp are not
affected):
_ the derivatives act as follow : if the input dzpuis is 0, the output
is 2; if the input is anything but 0 (n), the output is n+1. This is
also true for the covariant derivatives and the divergence.
_ the mult_r, div_r, mult_r_ced, div_r_ced, div_rsint, mult_rsint
methods do not change the dzpuis
_ mult_r_dzpuis, div_r_dzpuis, etc... change the dzpuis flag accordingly
to the argument given to hem

As a general rule, the more you use inc/dec_dzpuis back and forth, the
bigger the numerical noise introduced.

Best wishes,
--
Jerome Novak
Laboratoire de l'Univers et de ses THeories (LUTH)
Observatoire de Meudon       92195 Meudon FRANCE
tel: (33) 01 45 07 75 82     fax: (33) 01 45 07 79 71