connor/differential-equations
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

This should be the final

Browse Source
This commit is contained in:
Connor Johnstone
2023-10-16 14:49:15 -06:00
parent 3a15323a9c
commit 5d0a7d6e84
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/problem.rs
View File

@@ -1,5 +1,5 @@
use nalgebra::SVector; use nalgebra::SVector;
use roots::{find_root_brent, DebugConvergency}; use roots::{find_root_brent, SimpleConvergency};
use super::callback::Callback; use super::callback::Callback;
use super::controller::{Controller, PIController}; use super::controller::{Controller, PIController};
@@ -30,7 +30,7 @@ where
} }
} }
pub fn solve(&mut self) -> Solution<S, D> { pub fn solve(&mut self) -> Solution<S, D> {
let mut convergency = DebugConvergency::new(1e-12, 50); let mut convergency = SimpleConvergency { eps: 1e-10, max_iter: 1000 };
let mut times: Vec<f64> = vec![self.ode.t]; let mut times: Vec<f64> = vec![self.ode.t];
let mut states: Vec<SVector<f64, D>> = vec![self.ode.y]; let mut states: Vec<SVector<f64, D>> = vec![self.ode.y];
let mut dense_coefficients: Vec<Vec<SVector<f64, D>>> = Vec::new(); let mut dense_coefficients: Vec<Vec<SVector<f64, D>>> = Vec::new();