[LORENE] Corrected bug in angular Laplacian

[Jerome Novak]

Dear All,

I have just corrected an error in one of our functions calculating the 
angular Laplace operator (d2 / dtheta2 + 1/tan(theta) d / dtheta + 
1/sin^2(theta) d2 / dphi2); namely, in the file op_lapang.C, the 
function calculating the action of this operator on a Mtbl_cf was doing 
a wrong job for the following case:
_ base T_LEG_II (T_SIN_P + P_COSSIN_I) and,
_ m = 1 / sin(phi) part

In practice, this is a *rather non-standard* base in the sense that it 
corresponds to the super-symmetric case (symmetry with respect to the 
equatorial plane + pi-symmetry (x,y) -> (-x, -y) ), AND it is used only 
for some of the Cartesian components of a rank-2 tensor (xz and yz). In 
particular, this is of *NO CONCERN FOR BINARY STARS/BLACK HOLE 
CONFIGURATIONS*, *NOR FOR ROTATING NEUTRON STARS*.

I implemented this function (without testing it enough: mea culpa, mea 
maxima culpa!) in September 2003 and we did not use it until very 
recently. Now everything should work fine with this function ...
I hope this error did not introduce any problem in your studies.

Bon code,

-- 
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