Skip to content

Commit

Permalink
Use EvaluationFlags interface
Browse files Browse the repository at this point in the history
  • Loading branch information
bergbauer committed Jul 3, 2023
1 parent cc22d5b commit 6681f2b
Showing 1 changed file with 8 additions and 8 deletions.
16 changes: 8 additions & 8 deletions tests/matrix_free/ecl_fe_q.h
Original file line number Diff line number Diff line change
Expand Up @@ -150,10 +150,10 @@ test(const unsigned int geometry = 0,
{
phi.reinit(cell);
phi.read_dof_values(src);
phi.evaluate(false, true, false);
phi.evaluate(EvaluationFlags::gradients);
for (unsigned int q = 0; q < phi.n_q_points; ++q)
phi.submit_gradient(phi.get_gradient(q), q);
phi.integrate(false, true);
phi.integrate(EvaluationFlags::gradients);
phi.distribute_local_to_global(dst);
}
},
Expand All @@ -164,9 +164,9 @@ test(const unsigned int geometry = 0,
phi_p.reinit(face);

phi_m.read_dof_values(src);
phi_m.evaluate(true, true);
phi_m.evaluate(EvaluationFlags::values | EvaluationFlags::gradients);
phi_p.read_dof_values(src);
phi_p.evaluate(true, true);
phi_p.evaluate(EvaluationFlags::values | EvaluationFlags::gradients);
VectorizedArrayType sigmaF =
(std::abs(
(phi_m.normal_vector(0) * phi_m.inverse_jacobian(0))[dim - 1]) +
Expand All @@ -187,9 +187,9 @@ test(const unsigned int geometry = 0,
phi_m.submit_value(average_valgrad, q);
phi_p.submit_value(-average_valgrad, q);
}
phi_m.integrate(true, true);
phi_m.integrate(EvaluationFlags::values | EvaluationFlags::gradients);
phi_m.distribute_local_to_global(dst);
phi_p.integrate(true, true);
phi_p.integrate(EvaluationFlags::values | EvaluationFlags::gradients);
phi_p.distribute_local_to_global(dst);
}
},
Expand Down Expand Up @@ -241,10 +241,10 @@ test(const unsigned int geometry = 0,
{
phi.reinit(cell);
phi.read_dof_values(src);
phi.evaluate(false, true, false);
phi.evaluate(EvaluationFlags::gradients);
for (unsigned int q = 0; q < phi.n_q_points; ++q)
phi.submit_gradient(phi.get_gradient(q), q);
phi.integrate(false, true);
phi.integrate(EvaluationFlags::gradients);

for (unsigned int face = 0; face < GeometryInfo<dim>::faces_per_cell;
++face)
Expand Down

0 comments on commit 6681f2b

Please sign in to comment.