Commit 6f480760 authored by Dominic Kempf's avatar Dominic Kempf
Browse files

Fix jamps and averages

parent 6169e333
......@@ -392,11 +392,11 @@ class EulerBernoulli2DLocalOperator
// 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 * dt2vn_n_0 * sk_dtun));
r_n.accumulate(lfsu_n.child(1), i, E*I*(- flip_factor * sk_dt2un*dtvn_n_1 - flip_factor * 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));
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_s.accumulate(lfsu_s.child(0), i, -E*I*(- flip_factor * sk_dt2un*dtvn_s_0 - flip_factor * dt2vn_s_0 * sk_dtun));
r_s.accumulate(lfsu_s.child(1), i, -E*I*(- flip_factor * sk_dt2un*dtvn_s_1 - flip_factor * 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));
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_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