本文整理汇总了C++中Solver::LeapFrog方法的典型用法代码示例。如果您正苦于以下问题:C++ Solver::LeapFrog方法的具体用法?C++ Solver::LeapFrog怎么用?C++ Solver::LeapFrog使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Solver
的用法示例。
在下文中一共展示了Solver::LeapFrog方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: main
int main() {
vec A = zeros(values*N);
vec dAdt = zeros(values*N);
vec M = zeros(N);
double u, v, w;
double x, y, z;
double r, theta, phi;
long idum = -1;
// initial conditions
for (int i=0; i<N; i++) {
u = ran2(&idum);
v = ran2(&idum);
w = ran2(&idum);
phi = 2*pi*w;
r = R0*cbrt(u);
theta = acos(1-2*v);
x = r*sin(theta)*cos(phi);
y = r*sin(theta)*sin(phi);
z = r*cos(theta);
int j = values*i;
A(j) = x; // x - position
A(j+1) = y; // y - position
A(j+2) = z; // z - position
A(j+3) = 0.0; // vx - velocity
A(j+4) = 0.0; // vy - velocity
A(j+5) = 0.0; // vz - velocity
M(i) = sigma*gaussian_deviate(&idum) + mu;
}
Solver R = Solver(A, dAdt, M);
//R.RungeKutta4();
R.LeapFrog();
}