You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
27 lines
511 B
27 lines
511 B
9 years ago
|
#include <assert.h>
|
||
|
|
#include <stdio.h>
|
||
|
|
|
||
|
|
#include <matrix/integration.hpp>
|
||
|
|
|
||
|
|
using namespace matrix;
|
||
|
|
|
||
|
|
Vector<float, 6> f(float t, Vector<float, 6> y);
|
||
|
|
|
||
|
|
Vector<float, 6> f(float t, Vector<float, 6> y) {
|
||
|
|
return ones<float, 6, 1>();
|
||
|
|
}
|
||
|
|
|
||
|
|
int main()
|
||
|
|
{
|
||
|
|
Vector<float, 6> y = ones<float, 6, 1>();
|
||
|
|
float h = 1;
|
||
|
|
float t = 1;
|
||
|
|
y.T().print();
|
||
|
|
integrate_rk4(f, y, t, h);
|
||
|
|
y.T().print();
|
||
|
|
assert(y == (ones<float, 6, 1>()*2));
|
||
|
|
return 0;
|
||
|
|
}
|
||
|
|
|
||
|
|
/* vim: set et fenc=utf-8 ff=unix sts=0 sw=4 ts=4 : */
|