From a7c2084a14dd586cb396438bb0bcd711a6ca7074 Mon Sep 17 00:00:00 2001 From: Dave Goodwin Date: Mon, 17 May 2004 22:59:20 +0000 Subject: [PATCH] initial import --- Cantera/matlab/cantera/@ReactorNet/addReactor.m | 3 +++ Cantera/matlab/cantera/@ReactorNet/atol.m | 5 +++++ Cantera/matlab/cantera/@ReactorNet/rtol.m | 5 +++++ Cantera/matlab/cantera/@ReactorNet/setInitialTime.m | 9 +++++++++ Cantera/matlab/cantera/@ReactorNet/setMaxTimeStep.m | 12 ++++++++++++ Cantera/matlab/cantera/@ReactorNet/setTolerances.m | 7 +++++++ Cantera/matlab/cantera/@ReactorNet/time.m | 5 +++++ 7 files changed, 46 insertions(+) create mode 100644 Cantera/matlab/cantera/@ReactorNet/addReactor.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/atol.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/rtol.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/setInitialTime.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/setMaxTimeStep.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/setTolerances.m create mode 100644 Cantera/matlab/cantera/@ReactorNet/time.m diff --git a/Cantera/matlab/cantera/@ReactorNet/addReactor.m b/Cantera/matlab/cantera/@ReactorNet/addReactor.m new file mode 100644 index 000000000..b01a52147 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/addReactor.m @@ -0,0 +1,3 @@ +function addReactor(net, reactor) +% ADDREACTOR - Add a reactor to the network +reactornetmethods(4, reactornet_hndl(net), reactor_hndl(reactor)); diff --git a/Cantera/matlab/cantera/@ReactorNet/atol.m b/Cantera/matlab/cantera/@ReactorNet/atol.m new file mode 100644 index 000000000..388cd3145 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/atol.m @@ -0,0 +1,5 @@ +function t = atol(r) +% ATOL - absolute error tolerance +% +t = reactornetmethods(24, reactornet_hndl(r)); + diff --git a/Cantera/matlab/cantera/@ReactorNet/rtol.m b/Cantera/matlab/cantera/@ReactorNet/rtol.m new file mode 100644 index 000000000..6cf9e7690 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/rtol.m @@ -0,0 +1,5 @@ +function t = rtol(r) +% RTOL - relative error tolerance +% +t = reactornetmethods(23, reactornet_hndl(r)); + diff --git a/Cantera/matlab/cantera/@ReactorNet/setInitialTime.m b/Cantera/matlab/cantera/@ReactorNet/setInitialTime.m new file mode 100644 index 000000000..12cd32669 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/setInitialTime.m @@ -0,0 +1,9 @@ +function setInitialTime(r, t) +% SETINITIALTIME - Set the initial time at which to begin the time +% integration. +% +% If the time integration has already begun, this restarts the +% integrator using the current solution as the initial condition, +% setting the time to 't'. +% +reactornetmethods(5, reactornet_hndl(r), t); diff --git a/Cantera/matlab/cantera/@ReactorNet/setMaxTimeStep.m b/Cantera/matlab/cantera/@ReactorNet/setMaxTimeStep.m new file mode 100644 index 000000000..3ea28a02d --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/setMaxTimeStep.m @@ -0,0 +1,12 @@ +function setMaxTimeStep(r, maxstep) +% SETMAXTIMESTEP - Set the maximum time step. +% +% The integrator chooses a step size based on the desired error +% tolerances and the rate at which the solution is changing. In +% some cases, the solution changes very slowly at first, then very +% rapidly (ignition problems). In such cases, the integrator may +% choose a timestep that is too large, which leads to numerical +% problems later. Use this method to set an upper bound on the +% timestep. +% +reactornetmethods(6, reactornet_hndl(r), maxstep); diff --git a/Cantera/matlab/cantera/@ReactorNet/setTolerances.m b/Cantera/matlab/cantera/@ReactorNet/setTolerances.m new file mode 100644 index 000000000..978fbf167 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/setTolerances.m @@ -0,0 +1,7 @@ +function setTolerances(r, rtol, atol) +% SETTOLERANCES - Set error tolerances. +% +% rtol - scalar relative error tolerance +% atol - scalar absolute error tolerance +% +reactornetmethods(7, reactornet_hndl(r), rtol, atol); diff --git a/Cantera/matlab/cantera/@ReactorNet/time.m b/Cantera/matlab/cantera/@ReactorNet/time.m new file mode 100644 index 000000000..b93b94564 --- /dev/null +++ b/Cantera/matlab/cantera/@ReactorNet/time.m @@ -0,0 +1,5 @@ +function t = time(r) +% TIME - current value of the time. +% +t = reactornetmethods(22, reactornet_hndl(r)); +