laxwen1.hpp

面向对象的卡尔曼滤波器源码

// laxwen1.hpp		1-D Lax Wendroff (single step) method

//		u_t + f_x = rhs
//				where rhs is optional, assumed zero otherwise


#ifndef _LAX_WENDROFF1_HPP_
#define _LAX_WENDROFF1_HPP_

#include <dynamics.hpp>

#ifndef SPACE_DERIV
#define SPACE_DERIV
typedef float (*SpaceDeriv)(Vector&,int);
#endif

class LaxWendroff1 : public Dynamics
{
	private:
          float r;
	  Vector u_new;
	  SpaceDeriv f;
	  SpaceDeriv rhs;
	public:
	  LaxWendroff1(SpaceDeriv f, float dx_in, float dt_in, int n_in,
		      SpaceDeriv user_rhs = NULL);
	 ~LaxWendroff1() { }
	  void step(const int steps = 1);
};

#endif