Κεφάλαιο 1 - Nemertes

ΑΝΑΠΤΥΞΗ ΚΩ∆ΙΚΑ ΜΟΡΙΑΚΗΣ ∆ΥΝΑΜΙΚΗΣ
Η συνάρτηση 66 που παρουσιάστηκε παραπάνω υλοποιείται εδώ:
/**
* Eq. 66
* @param coordinates
* @return
*/
public double getEnergyForPosition(double[] coordinates)
{
double numberOfAbsorbedPhotons;
double rz = r(coordinates[2]);
double x = coordinates[0];
double y = coordinates[1];
//Calculate the Number of photons at this depth using the Beer - Lambert Law
double firstTerm = getLaserEnergy() *
MDParameters.timeStep;//*getΛ()/(Constants.h*Constants.c);
numberOfAbsorbedPhotons = firstTerm;
//commented out - maybe will be added again
// double secondTerm = (rz * rz) / (rf0 * rf0);
// numberOfAbsorbedPhotons *= secondTerm;
double thirdTerm = Math.exp(-β * coordinates[2]);
numberOfAbsorbedPhotons *= thirdTerm;
//apply the bivariate gaussian distribution
double probability;
double d = ac * (x + material.particleRadius) / (2 * radius);
double e = ac * (x - material.particleRadius) / (2 * radius);
probability = Φ(d) - Φ(e);
double e2 = ac * (y + material.particleRadius) / (2 * radius);
double f = ac * (y - material.particleRadius) / (2 * radius);
probability *= Φ(e2) - Φ(f);
numberOfAbsorbedPhotons *= probability;
numberOfAbsorbedPhotons *= material.reflectivity;
if (("" + numberOfAbsorbedPhotons).equals("NaN"))
{
System.out.println("Energy is NaN!!!");
System.out.println("Coordinates: ( " + coordinates[0] + " , " + coordinates[1] + " , "
+ coordinates[2] + " )");
System.out.println("Rz = " + rz);
System.out.println("Probability: " + probability);
System.out.println("d = " + d);
System.out.println("e = " + e);
System.out.println("e2 = " + e2);
System.out.println("f = " + f);
}
//multiply with the energy of one photon hc/λ
//we don't have to multiply with the energy of one photon,
//as we have removed it from the initial calculation of the Beer Lambert Law
return numberOfAbsorbedPhotons;
}
238