toastpp Slab

Using TOAST++, the DOT problem is solved in a slab geometry. The forward solution at various detector positions is compared to the analytical solution to the diffusion equation. The analytical solution was calculated using different boundary conditions than those used by TOAST++, so the solutions are similar but not exactly the same.

The code is written in MATLAB, and the steps are split into .m files. The steps are:

1. Generate slab nodes
2a. Generate mesh with tetgen
2b. Check mesh and remove bad elements
3. Generate initial values
4a. Generate detector positions
4b. Generate source positions
5. Generate source-detector (QM) file
6a. Solve the forward solution in TOAST++ and compare to the analytical solution
6b. Generate the FEM file from synthetic data generated in 6a
7. Reconstruct mus and mua in the slab using TOAST++
8. Plot solution

tetgen: http://wias-berlin.de/software/tetgen/

This code is a starting point and can be expanded to complicated geometries.