Reply to Runge-Kutta-Verfahren on Wed, 12 Oct 2005 16:36:56 GMT

fubar — Wed, 12 Oct 2005 16:36:56 GMT

Hallo,

ich wollte gerade meine Klassen zum Lösen gewöhnlicher DGL'en überarbeiten, dazu gehören auch Runge-Kutta-Verfahren:

\par u\_0 = y\_0 \par u_{j+1} = u\_j + h \sum\_{i=1}^m A\_i k\_i \par k\_1 = f(t\_j, y_j) \par k\_2 = f(t\_j + a\_2 h, u\_j + h b_{21}k_1) \par k\_3 = f(t\_j + a\_3 h, u\_j + h(b_{31}k\_1 + b\_{32}k_2) \par ... \par k\_m = f(t\_j + a\_m h, u\_j + h \sum_{i=1}^{m-1} b_{mi}k_i) \par

Allerdings will ich für die verschiedenen Verfahren nicht immer alles neu schreiben, deshalb habe ich das bis jetzt so gelöst:

template 
	class KuttaNystroem
	{
		protected:
			static const size_t order;
			static const T A[];
			static const T a[];
			static const T b[];
	};

	template 
	const size_t KuttaNystroem::order = 6;

	template 
	const T KuttaNystroem::A[] =
	{ ... }; 

	template 
	const T KuttaNystroem::a[] =
	{ ... };

	template 
	const T KuttaNystroem::b[] =
	{ ... };
	//------------------------------

	template 
	class England1
	{
		...
	}

Und die Implementierung des Verfahrens:

template
	<
		typename T,
		template class Method = ClassicRK
	>
	class RungeKutta : public Method
	{
		... 
	};

Ist das so sinnvoll, oder könnte man das irgendwie besser lösen? Wie würdet ihr das machen? Es geht mir vor allem die Speicherung der Koeffizienten.

Danke im voraus!