[LORENE] Corrected bug in angular Laplacian
Jerome.Novak at obspm.fr
Wed May 18 11:19:33 CEST 2005
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.
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
More information about the Lorene.list