@jmcragun, take a look at this. Two commits here, the first sets up the advection diffusion integrator to register platelets and the activating chemical. It also sets up the cohesion stress variable. The second commit sets up the beta function and registers it with the forcing function.
Note that the way I've implemented the beta function requires that the setter in CohesionStressRHS to take a function of the form std::function<double(double, void*)>. You could also rewrite it to use std::bind().
@jmcragun, take a look at this. Two commits here, the first sets up the advection diffusion integrator to register platelets and the activating chemical. It also sets up the cohesion stress variable. The second commit sets up the beta function and registers it with the forcing function.
Note that the way I've implemented the beta function requires that the setter in CohesionStressRHS to take a function of the form
std::function<double(double, void*)>
. You could also rewrite it to usestd::bind()
.