integrator.h

Go to the documentation of this file.

/*
 * integrator.h - integrator class definitions
 *
 * Copyright (C) 2004, 2005, 2006 Stefan Jahn <stefan@lkcc.org>
 *
 * This is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2, or (at your option)
 * any later version.
 * 
 * This software is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this package; see the file COPYING.  If not, write to
 * the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
 * Boston, MA 02110-1301, USA.  
 *
 * $Id: integrator.h 1825 2011-03-11 20:42:14Z ela $
 *
 */

#ifndef __INTEGRATOR_H__
#define __INTEGRATOR_H__

#include "states.h"

#define MODE_NONE 0
#define MODE_INIT 1

class integrator : public states<nr_double_t>
{
 public:
  // constructor and destructor set
  integrator ();
  integrator (const integrator &);
  ~integrator ();

  // integration specific
  typedef void (* integrate_func_t)
    (integrator *, int, nr_double_t, nr_double_t&, nr_double_t&);
  void setIntegration (integrate_func_t f) { integrate_func = f; }
  typedef void (* conductor_func_t)
    (integrator *, nr_double_t, nr_double_t&);
  void setConductance (conductor_func_t f) { conductor_func = f; }
  void integrate (int, nr_double_t, nr_double_t&, nr_double_t&);
  void conductor (nr_double_t, nr_double_t&);
  void setOrder (int o) { order = o; }
  int  getOrder (void) { return order; }
  void setMode (int s) { state = s; }
  int  getMode (void) { return state; }
  void setCoefficients (nr_double_t * c) { coefficients = c; }
  nr_double_t * getCoefficients (void) { return coefficients; }

 private:
  int order;
  int state;
  nr_double_t * coefficients;
  integrate_func_t integrate_func;
  conductor_func_t conductor_func;
};

#endif /* __INTEGRATOR_H__ */