/// <summary>
/// Method to determine if photon datapoint should be tallied or not
/// </summary>
/// <param name="dp"></param>
/// <param name="photonDatabaseWriter"></param>
/// <returns></returns>
public bool DPBelongsToSurfaceVirtualBoundary(PhotonDataPoint dp,
PhotonDatabaseWriter photonDatabaseWriter)
{
if ((dp.StateFlag.HasFlag(PhotonStateType.PseudoDiffuseReflectanceVirtualBoundary) && // pMC uses regular PST
photonDatabaseWriter.VirtualBoundaryType == VirtualBoundaryType.pMCDiffuseReflectance))
{
return(true);
}
return(false);
}
public void validate_PhotonDatabase_deserialized_class_is_correct_when_using_WriteToFile()
{
// test serialization
new SimulationInput().ToFile("SimulationInputTest.txt");
string databaseFilename = "testphotondatabase";
using (var dbWriter = new PhotonDatabaseWriter(
VirtualBoundaryType.DiffuseReflectance, databaseFilename))
{
dbWriter.Write(new PhotonDataPoint(
new Position(1, 2, 3),
new Direction(0, 0, 1),
1.0, // weight
10, // time
PhotonStateType.None));
dbWriter.Write(new PhotonDataPoint(
new Position(4, 5, 6),
new Direction(1, 0, 0),
0.50,
100,
PhotonStateType.None));
}
// read the database from file, and verify the correct number of photons were written
var dbCloned = PhotonDatabase.FromFile(databaseFilename);
Assert.AreEqual(dbCloned.NumberOfElements, 2);
// manually enumerate through the first two elements (same as foreach)
// PhotonDatabase is designed so you don't have to have the whole thing
// in memory, so .ToArray() loses the benefits of the lazy-load data points
var enumerator = dbCloned.DataPoints.GetEnumerator();
// advance to the first point and test that the point is valid
enumerator.MoveNext();
var dp1 = enumerator.Current;
Assert.AreEqual(dp1.Position, new Position(1, 2, 3));
Assert.AreEqual(dp1.Direction, new Direction(0, 0, 1));
Assert.AreEqual(dp1.Weight, 1.0);
Assert.AreEqual(dp1.TotalTime, 10);
Assert.IsTrue(dp1.StateFlag.HasFlag(PhotonStateType.None));
// advance to the second point and test that the point is valid
enumerator.MoveNext();
var dp2 = enumerator.Current;
Assert.AreEqual(dp2.Position, new Position(4, 5, 6));
Assert.AreEqual(dp2.Direction, new Direction(1, 0, 0));
Assert.AreEqual(dp2.Weight, 0.5);
Assert.AreEqual(dp2.TotalTime, 100);
Assert.IsTrue(dp2.StateFlag.HasFlag(PhotonStateType.None));
}
public void validate_pMCDatabase_deserialized_class_is_correct_when_using_WriteToFile()
{
// test serialization
new SimulationInput().ToFile("SimulationInputTest.txt");
string photonDbFilename = "testpmcdatabase";
string collisionDbFilename = "testpmccollisiondatabase";
int numberOfSubregions = 3;
using (var dbWriter = new PhotonDatabaseWriter(
VirtualBoundaryType.pMCDiffuseReflectance, photonDbFilename))
{
using (var collisionDbWriter = new CollisionInfoDatabaseWriter(
VirtualBoundaryType.pMCDiffuseReflectance, collisionDbFilename, numberOfSubregions))
{
// write data for first photon (exit DP and collision info)
dbWriter.Write(new PhotonDataPoint(
new Position(1, 2, 3),
new Direction(0, 0, 1),
1.0, // weight
10, // time
PhotonStateType.None));
collisionDbWriter.Write(new CollisionInfo(numberOfSubregions)
{
new SubRegionCollisionInfo(10.0, 1000),
new SubRegionCollisionInfo(20.0, 2000),
new SubRegionCollisionInfo(30.0, 3000)
});
// write data for second photon (exit DP and collision info)
dbWriter.Write(new PhotonDataPoint(
new Position(4, 5, 6),
new Direction(1, 0, 0),
0.50,
100,
PhotonStateType.None));
collisionDbWriter.Write(new CollisionInfo(numberOfSubregions)
{
new SubRegionCollisionInfo(40.0, 4000),
new SubRegionCollisionInfo(50.0, 5000),
new SubRegionCollisionInfo(60.0, 6000)
});
}
}
// read the database from file, and verify the correct number of photons were written
var dbCloned = pMCDatabase.FromFile(photonDbFilename, collisionDbFilename);
Assert.AreEqual(dbCloned.DataPoints.Count(), 2);
// manually enumerate through the first two elements (same as foreach)
// PhotonDatabase is designed so you don't have to have the whole thing
// in memory, so .ToArray() loses the benefits of the lazy-load data points
var enumerator = dbCloned.DataPoints.GetEnumerator();
// advance to the first point and test that the point is valid
enumerator.MoveNext();
var dp1 = enumerator.Current;
// verify photon database entries for first photon
Assert.AreEqual(dp1.PhotonDataPoint.Position, new Position(1, 2, 3));
Assert.AreEqual(dp1.PhotonDataPoint.Direction, new Direction(0, 0, 1));
Assert.AreEqual(dp1.PhotonDataPoint.Weight, 1.0);
Assert.AreEqual(dp1.PhotonDataPoint.TotalTime, 10);
Assert.IsTrue(dp1.PhotonDataPoint.StateFlag.HasFlag(PhotonStateType.None));
// verify collision info for first photon
Assert.AreEqual(dp1.CollisionInfo[0].PathLength, 10.0);
Assert.AreEqual(dp1.CollisionInfo[0].NumberOfCollisions, 1000);
Assert.AreEqual(dp1.CollisionInfo[1].PathLength, 20.0);
Assert.AreEqual(dp1.CollisionInfo[1].NumberOfCollisions, 2000);
Assert.AreEqual(dp1.CollisionInfo[2].PathLength, 30.0);
Assert.AreEqual(dp1.CollisionInfo[2].NumberOfCollisions, 3000);
// advance to the second point and test that the point is valid
enumerator.MoveNext();
// verify photon database entries for second photon
var dp2 = enumerator.Current;
Assert.AreEqual(dp2.PhotonDataPoint.Position, new Position(4, 5, 6));
Assert.AreEqual(dp2.PhotonDataPoint.Direction, new Direction(1, 0, 0));
Assert.AreEqual(dp2.PhotonDataPoint.Weight, 0.5);
Assert.AreEqual(dp2.PhotonDataPoint.TotalTime, 100);
Assert.IsTrue(dp2.PhotonDataPoint.StateFlag.HasFlag(PhotonStateType.None));
// verify collision info for second photon
Assert.AreEqual(dp2.CollisionInfo[0].PathLength, 40.0);
Assert.AreEqual(dp2.CollisionInfo[0].NumberOfCollisions, 4000);
Assert.AreEqual(dp2.CollisionInfo[1].PathLength, 50.0);
Assert.AreEqual(dp2.CollisionInfo[1].NumberOfCollisions, 5000);
Assert.AreEqual(dp2.CollisionInfo[2].PathLength, 60.0);
Assert.AreEqual(dp2.CollisionInfo[2].NumberOfCollisions, 6000);
}
