size_t GramSchmidt(size_t N, size_t Mold, size_t Mnew, F *V, size_t LDV,
CQMemManager &mem, size_t NRe = 0, double eps = 1e-12) {
return
GramSchmidt(N,Mold,Mnew,V,LDV,
[&](F* Vc){ return std::sqrt(std::abs(InnerProd<F>(N,Vc,1,Vc,1))); },
[&](size_t i, size_t j, F* Vi, size_t LDVi, F* Vj, size_t LDVj,
F* inner){
Gemm('C','N',i,j,N,F(1.),Vi,LDVi,Vj,LDVj,F(0.),inner,i);
},
mem,NRe,eps);
};
container LyapunovSpectrum::LocalLyapunov(){
++(*fiducial);
labels_values result(fiducial->get_labels_values());
for(auto i : result)
jacobian_parameters[i.first] = i.second;
//Integrating the linear equations and saving the vector on ortogonal_space
std::unique_ptr<FunctionCapsule> function(jacobian->Create(jacobian_parameters));
for(std::vector<labels_values>::iterator it(base.begin()); it != base.end(); ++it){
RungeKutta4Th linearized_sistem(*function, *it, fiducial->get_dt());
++linearized_sistem;
*it = linearized_sistem.get_labels_values();
}
return GramSchmidt();
}
