summaryrefslogtreecommitdiff
path: root/src/c/differential_calculus/ode/doded.c
diff options
context:
space:
mode:
authorsiddhu89902016-05-27 11:50:04 +0530
committersiddhu89902016-05-27 11:50:04 +0530
commit212c54043e454e41ca4b23a5b965d4dbd4b683fe (patch)
treec1b62c677facad021562292b823b5e21978c0e4b /src/c/differential_calculus/ode/doded.c
parenta6d6a9c1e88b75668868af691c9731075a514ffb (diff)
downloadScilab2C_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.c39
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);
-}