Commit 89ae5390 authored by Dominic Kempf's avatar Dominic Kempf
Browse files

Fix the accumulation in alpha_skeleton

parent 6f480760
......@@ -386,17 +386,17 @@ class EulerBernoulli2DLocalOperator
auto dt2vn_s_1 = (t[1] * t[0] * (jit_s[1][1] * (jit_s[1][1] * basis_s.hessian(i, 1, 1) + jit_s[1][0] * basis_s.hessian(i, 1, 0)) + jit_s[1][0] * (jit_s[1][1] * basis_s.hessian(i, 0, 1) + jit_s[1][0] * basis_s.hessian(i, 0, 0))) + t[0] * t[0] * (jit_s[0][1] * (jit_s[1][1] * basis_s.hessian(i, 1, 1) + jit_s[1][0] * basis_s.hessian(i, 1, 0)) + jit_s[0][0] * (jit_s[1][1] * basis_s.hessian(i, 0, 1) + jit_s[1][0] * basis_s.hessian(i, 0, 0)))) * t[1] + (t[1] * t[0] * (jit_s[1][1] * (jit_s[0][1] * basis_s.hessian(i, 1, 1) + jit_s[0][0] * basis_s.hessian(i, 1, 0)) + jit_s[1][0] * (jit_s[0][1] * basis_s.hessian(i, 0, 1) + jit_s[0][0] * basis_s.hessian(i, 0, 0))) + t[0] * t[0] * (jit_s[0][1] * (jit_s[0][1] * basis_s.hessian(i, 1, 1) + jit_s[0][0] * basis_s.hessian(i, 1, 0)) + jit_s[0][0] * (jit_s[0][1] * basis_s.hessian(i, 0, 1) + jit_s[0][0] * basis_s.hessian(i, 0, 0)))) * t[0];
double flip_factor = flipped ? -1.0 : 1.0;
// r_n.accumulate(lfsu_n.child(0), i, E*I*( penalty * sk_dtun * dtvn_n_0));
// r_n.accumulate(lfsu_n.child(1), i, E*I*( penalty * sk_dtun * dtvn_n_1));
// r_n.accumulate(lfsu_n.child(0), i, -E*I*( penalty * sk_dtun * dtvn_n_0));
// r_n.accumulate(lfsu_n.child(1), i, -E*I*( penalty * sk_dtun * dtvn_n_1));
// r_s.accumulate(lfsu_s.child(0), i, -E*I*( penalty * sk_dtun * dtvn_s_0));
// r_s.accumulate(lfsu_s.child(1), i, -E*I*( penalty * sk_dtun * dtvn_s_1));
// r_s.accumulate(lfsu_s.child(0), i, E*I*( penalty * sk_dtun * dtvn_s_0));
// r_s.accumulate(lfsu_s.child(1), i, E*I*( penalty * sk_dtun * dtvn_s_1));
r_n.accumulate(lfsu_n.child(0), i, -E*I*(- flip_factor * sk_dt2un*dtvn_n_0 + flip_factor * 0.5 * dt2vn_n_0 * sk_dtun));
r_n.accumulate(lfsu_n.child(1), i, -E*I*(- flip_factor * sk_dt2un*dtvn_n_1 + flip_factor * 0.5 * dt2vn_n_1 * sk_dtun));
r_n.accumulate(lfsu_n.child(0), i, -E*I*(- flip_factor * sk_dt2un*dtvn_n_0 + flip_factor * 0.5 * dt2vn_n_0 * sk_dtun + penalty * sk_dtun * dtvn_n_0));
r_n.accumulate(lfsu_n.child(1), i, -E*I*(- flip_factor * sk_dt2un*dtvn_n_1 + flip_factor * 0.5 * dt2vn_n_1 * sk_dtun + penalty * sk_dtun * dtvn_n_1));
r_s.accumulate(lfsu_s.child(0), i, -E*I*(flip_factor * sk_dt2un*dtvn_s_0 + flip_factor * 0.5 * dt2vn_s_0 * sk_dtun));
r_s.accumulate(lfsu_s.child(1), i, -E*I*(flip_factor * sk_dt2un*dtvn_s_1 + flip_factor * 0.5 * dt2vn_s_1 * sk_dtun));
r_s.accumulate(lfsu_s.child(0), i, -E*I*(flip_factor * sk_dt2un*dtvn_s_0 + flip_factor * 0.5 * dt2vn_s_0 * sk_dtun - penalty * sk_dtun * dtvn_s_0));
r_s.accumulate(lfsu_s.child(1), i, -E*I*(flip_factor * sk_dt2un*dtvn_s_1 + flip_factor * 0.5 * dt2vn_s_1 * sk_dtun - penalty * sk_dtun * dtvn_s_1));
// r_n.accumulate(lfsu_n.child(0), i, E*I*(- flip_factor * sk_dt2un*dtvn_n_0 - flip_factor * dt2vn_n_0 * sk_dtun + penalty * sk_dtun * dtvn_n_0));
// r_n.accumulate(lfsu_n.child(1), i, E*I*(- flip_factor * sk_dt2un*dtvn_n_1 - flip_factor * dt2vn_n_1 * sk_dtun + penalty * sk_dtun * dtvn_n_1));
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment