diff options
author | siddhu8990 | 2016-05-27 11:50:04 +0530 |
---|---|---|
committer | siddhu8990 | 2016-05-27 11:50:04 +0530 |
commit | 212c54043e454e41ca4b23a5b965d4dbd4b683fe (patch) | |
tree | c1b62c677facad021562292b823b5e21978c0e4b /src/c/differential_calculus/ode/doded.c | |
parent | a6d6a9c1e88b75668868af691c9731075a514ffb (diff) | |
download | Scilab2C_fossee_old-212c54043e454e41ca4b23a5b965d4dbd4b683fe.tar.gz Scilab2C_fossee_old-212c54043e454e41ca4b23a5b965d4dbd4b683fe.tar.bz2 Scilab2C_fossee_old-212c54043e454e41ca4b23a5b965d4dbd4b683fe.zip |
Support for basic ODE function added
Diffstat (limited to 'src/c/differential_calculus/ode/doded.c')
-rw-r--r-- | src/c/differential_calculus/ode/doded.c | 39 |
1 files changed, 0 insertions, 39 deletions
diff --git a/src/c/differential_calculus/ode/doded.c b/src/c/differential_calculus/ode/doded.c deleted file mode 100644 index b544c42..0000000 --- a/src/c/differential_calculus/ode/doded.c +++ /dev/null @@ -1,39 +0,0 @@ -// Copyright (C) 2016 - IIT Bombay - FOSSEE -// -// This file must be used under the terms of the CeCILL. -// This source file is licensed as described in the file COPYING, which -// you should have received as part of this distribution. The terms -// are also available at -// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt -// Author: Siddhesh Wani -// Organization: FOSSEE, IIT Bombay -// Email: toolbox@scilab.in - -#include "ode.h" -#include "types.h" -#include <gsl/gsl_errno.h> -#include <gsl/gsl_matrix.h> -#include <gsl/gsl_odeiv2.h> - - -void doded(double *initial_value, double start_time, double end_time, \ - char *ode_function, double nequs, double eps_abs, double eps_rel, \ - double step_size, double *out) -{ - double t = start_time; - //Setup ODE related parameters - gsl_odeiv2_system sys = {ode_function, NULL, 2, NULL}; - - gsl_odeiv2_step *s = gsl_odeiv2_step_alloc (gsl_odeiv2_step_rkf45, nequs); - gsl_odeiv2_control *c = gsl_odeiv2_control_y_new (eps_abs, eps_rel); - gsl_odeiv2_evolve *e = gsl_odeiv2_evolve_alloc (nequs); - - while(t < end_time) - { - gsl_odeiv2_evolve_apply_fixed_step (e, c, s, &sys, &t, step_size, &out); - } - - gsl_odeiv2_evolve_free (e); - gsl_odeiv2_control_free (c); - gsl_odeiv2_step_free (s); -} |