search

gustavdelius / mizer

This is my old mizer fork that is not actively used any more but still has lots of material in its many branches.
http://gustavdelius.github.io/mizer/
0 stars 2 forks source link

Use same discretisation scheme in set_scaling_model as in project #69

Closed gustavdelius closed 3 years ago

gustavdelius commented 6 years ago

Currently the discretised steady-state solution produced in set_scaling model and addSpecies interprets n[i] as the abundance between w[i] and w[i+1] and as a consequence sets n equal to zero at the maximal weight of a species. The upwind-difference scheme used in project() does not do this.

gustavdelius commented 3 years ago

The scaling model is now set up by 'newTraitParams()' with the option 'perfect_scaling = TRUE'.