@@ -527,37 +527,12 @@ namespace ntt {
527527 domain.mesh .extent (in::x2).first ,
528528 domain.mesh .extent (in::x3).first };
529529 const auto ext_current = m_pgen.ext_current ;
530- const auto dx = domain.mesh .metric .template sqrt_h_ <1 , 1 >({});
531- auto J = domain.fields .cur ;
532- real_t jx1 { ZERO };
533- real_t jx1_after { ZERO };
534- // clang-format off
535- Kokkos::parallel_reduce (
536- " Ampere" ,
537- domain.mesh .rangeActiveCells (),
538- KOKKOS_LAMBDA (index_t i1, real_t & val) {
539- val += J (i1, cur::jx1);
540- },
541- jx1);
542- jx1 /= domain.mesh .n_active (in::x1);
543- jx1 /= ppc0;
544530 Kokkos::parallel_for (
545531 " Ampere" ,
546532 domain.mesh .rangeActiveCells (),
547533 kernel::mink::CurrentsAmpere_kernel<M::Dim, decltype (ext_current)>(
548534 domain.fields .em , domain.fields .cur ,
549535 coeff, ppc0, ext_current, xmin, dx));
550- Kokkos::parallel_reduce (
551- " Ampere" ,
552- domain.mesh .rangeActiveCells (),
553- KOKKOS_LAMBDA (index_t i1, real_t & val) {
554- val += J (i1, cur::jx1);
555- },
556- jx1_after);
557- jx1_after /= domain.mesh.n_active(in::x1);
558- // clang-format on
559- std::cout << " jx1_before: " << jx1 << std::endl;
560- std::cout << " jx1_after: " << jx1_after << std::endl;
561536 } else {
562537 Kokkos::parallel_for (
563538 " Ampere" ,
0 commit comments