search

waternumbers / dynatop

Dynamic TOPMODEL
https://waternumbers.github.io/dynatop/
3 stars 6 forks source link

Improve implimentation of the eigen value method for ODEs #4

Closed waternumbers closed 4 years ago

waternumbers commented 5 years ago

Current the eigen value method computes the constants in the solution using the eigen matrix M and initial condition x_0 as

c <- solve(M) %*% x_0

since solve(M) which computes the inverse of M need only be evaluated once.

However this is less numerically stable then using

c <- solve(M,x_0)

i suggest changing this to improve the stability of the solution

waternumbers commented 4 years ago

This calculation is removed from the code in v0.0.4 and onwards