Nullfeldaufspaltung von Benzol und Naphthalin im ... - ScienceUp.de

A Quelltexte der C++-Programm 137
pow(rgd,2)/(alphaI + alphaJ + alphaK + alphaL);
result = sqrt(result);
/*FEHLERTEST*/
FehlerBehandlung.IstGroesser(result, MAX, "r0",
i, j, k, l);
FehlerBehandlung.SehrKlein(pow(result,4),
100.0*DBL_EPSILON,
"r0", i, j, k, l); //wegen 1/r0^4-Term
return result;
} //R0
/*********************************************************/
double RGammaDelta(double * RGamma, double * RDelta,
USHORT i, USHORT j, USHORT k, USHORT l)
//return rgd=Betrag des Abstandsvektors RGamma-RDelta
{
double result = 0.0;
USHORT count, NullSetzen[2];
Fehler FehlerBehandlung;
/*FEHLERTEST*/
//In Faellen, in denen rgd==0 sein sollte, ist wegen
// der Ungenauigkeit ab der
//16. Stelle rgd in der Groessenordnung von 1e-32, d.h
// rgd ist dann i.a. !=0 und < DBL_EPSILON.
//Da dies aber beim Test von rgd beruecksichtigt wird,
// muss der Fall hier nicht mehr beruecksichtigt
// werden. => NichtNull statt Ausloeschung
for (count = 0; count < 2; count++)
//z-Komponenten immer gleich Null
FehlerBehandlung.NichtNull(RGamma[count] -
RDelta[count], 10.0*DBL_EPSILON, 100.0*DBL_EPSILON,
"Eine Komponente von RGammaDelta", i, j, k, l);
for (count = 0; count < 2; count++) {
NullSetzen[count] = (fabs(RGamma[count] -
RDelta[count])