Good job indeed! Im looking forward on integration!!
As soon as the 0.12 is history we should go right after that!
Code: Select all
// check if solving works
if (rel_error < 1e-5) {
//new maxStep scaling
double scale = subsys->maxStep() / h.lpNorm<Eigen::Infinity>();
if( scale < 1.)
h *= scale;
// compute par's new estimate and ||d_par||^2
x_new = x + h;
double h_norm = h.squaredNorm();
Code: Select all
// see if we are already finished
if (stop)
break;
double scale = subsys->maxStep() / h_dl.lpNorm<Eigen::Infinity>();
if ( scale < 1.)
h_dl *= scale;
// get the new values
Code: Select all
solver inputfile outputfile
That's the best notice anyway!ickby wrote:So whats the sense of this post then? I don't know, maybe the notice that I'm still working on it