Skip to content

Commit d1c5304

Browse files
ykempfursg
ykempf
authored and
ursg
committed
Turn off clang format in small part of Hall term.
one more space Y and Z clang format bypass in Hall term too.
1 parent cb97de2 commit d1c5304

File tree

1 file changed

+104
-74
lines changed

1 file changed

+104
-74
lines changed

‎fieldsolver/ldz_hall.cpp‎

Lines changed: 104 additions & 74 deletions
Original file line numberDiff line numberDiff line change
@@ -472,39 +472,49 @@ void calculateEdgeHallTermXComponents(
472472
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EXHALL_011_111) = EXHall;
473473

474474
break;
475-
case 2:
475+
case 2:
476+
// clang-format off
476477
hallRhoq = FOURTH * (
477-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
478-
momentsGrid.get(i ,j-1,k )->at(fsgrids::moments::RHOQ) +
479-
momentsGrid.get(i ,j ,k-1)->at(fsgrids::moments::RHOQ) +
480-
momentsGrid.get(i ,j-1,k-1)->at(fsgrids::moments::RHOQ)
478+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
479+
+ momentsGrid.get(i, j - 1, k)->at(fsgrids::moments::RHOQ)
480+
+ momentsGrid.get(i, j, k - 1)->at(fsgrids::moments::RHOQ)
481+
+ momentsGrid.get(i, j - 1, k - 1)->at(fsgrids::moments::RHOQ)
481482
);
482-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
483-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EXHALL_000_100) = JXBX_000_100(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
483+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
484+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EXHALL_000_100) =
485+
JXBX_000_100(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
486+
(physicalconstants::MU_0 * hallRhoq);
484487
hallRhoq = FOURTH * (
485-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
486-
momentsGrid.get(i ,j+1,k )->at(fsgrids::moments::RHOQ) +
487-
momentsGrid.get(i ,j ,k-1)->at(fsgrids::moments::RHOQ) +
488-
momentsGrid.get(i ,j+1,k-1)->at(fsgrids::moments::RHOQ)
488+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
489+
+ momentsGrid.get(i, j + 1, k)->at(fsgrids::moments::RHOQ)
490+
+ momentsGrid.get(i, j, k - 1)->at(fsgrids::moments::RHOQ)
491+
+ momentsGrid.get(i, j + 1, k - 1)->at(fsgrids::moments::RHOQ)
489492
);
490-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
491-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EXHALL_010_110) = JXBX_010_110(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
493+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
494+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EXHALL_010_110) =
495+
JXBX_010_110(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
496+
(physicalconstants::MU_0 * hallRhoq);
492497
hallRhoq = FOURTH * (
493-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
494-
momentsGrid.get(i ,j-1,k )->at(fsgrids::moments::RHOQ) +
495-
momentsGrid.get(i ,j ,k+1)->at(fsgrids::moments::RHOQ) +
496-
momentsGrid.get(i ,j-1,k+1)->at(fsgrids::moments::RHOQ)
498+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
499+
+ momentsGrid.get(i, j - 1, k)->at(fsgrids::moments::RHOQ)
500+
+ momentsGrid.get(i, j, k + 1)->at(fsgrids::moments::RHOQ)
501+
+ momentsGrid.get(i, j - 1, k + 1)->at(fsgrids::moments::RHOQ)
497502
);
498-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
499-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EXHALL_001_101) = JXBX_001_101(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
503+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
504+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EXHALL_001_101) =
505+
JXBX_001_101(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
506+
(physicalconstants::MU_0 * hallRhoq);
500507
hallRhoq = FOURTH * (
501-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
502-
momentsGrid.get(i ,j+1,k )->at(fsgrids::moments::RHOQ) +
503-
momentsGrid.get(i ,j ,k+1)->at(fsgrids::moments::RHOQ) +
504-
momentsGrid.get(i ,j+1,k+1)->at(fsgrids::moments::RHOQ)
508+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
509+
+ momentsGrid.get(i, j + 1, k)->at(fsgrids::moments::RHOQ)
510+
+ momentsGrid.get(i, j, k + 1)->at(fsgrids::moments::RHOQ)
511+
+ momentsGrid.get(i, j + 1, k + 1)->at(fsgrids::moments::RHOQ)
505512
);
506-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
507-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EXHALL_011_111) = JXBX_011_111(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
513+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
514+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EXHALL_011_111) =
515+
JXBX_011_111(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
516+
(physicalconstants::MU_0 * hallRhoq);
517+
// clang-format on
508518
break;
509519

510520
default:
@@ -570,39 +580,49 @@ void calculateEdgeHallTermYComponents(
570580
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EYHALL_001_011) = EYHall;
571581
break;
572582

573-
case 2:
583+
case 2:
584+
// clang-format off
574585
hallRhoq = FOURTH * (
575-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
576-
momentsGrid.get(i-1,j ,k )->at(fsgrids::moments::RHOQ) +
577-
momentsGrid.get(i ,j ,k-1)->at(fsgrids::moments::RHOQ) +
578-
momentsGrid.get(i-1,j ,k-1)->at(fsgrids::moments::RHOQ)
586+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
587+
+ momentsGrid.get(i - 1, j, k)->at(fsgrids::moments::RHOQ)
588+
+ momentsGrid.get(i, j, k - 1)->at(fsgrids::moments::RHOQ)
589+
+ momentsGrid.get(i - 1, j, k - 1)->at(fsgrids::moments::RHOQ)
579590
);
580-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
581-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EYHALL_000_010) = JXBY_000_010(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
591+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
592+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EYHALL_000_010) =
593+
JXBY_000_010(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
594+
(physicalconstants::MU_0 * hallRhoq);
582595
hallRhoq = FOURTH * (
583-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
584-
momentsGrid.get(i+1,j ,k )->at(fsgrids::moments::RHOQ) +
585-
momentsGrid.get(i ,j ,k-1)->at(fsgrids::moments::RHOQ) +
586-
momentsGrid.get(i+1,j ,k-1)->at(fsgrids::moments::RHOQ)
596+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
597+
+ momentsGrid.get(i + 1, j, k)->at(fsgrids::moments::RHOQ)
598+
+ momentsGrid.get(i, j, k - 1)->at(fsgrids::moments::RHOQ)
599+
+ momentsGrid.get(i + 1, j, k - 1)->at(fsgrids::moments::RHOQ)
587600
);
588-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
589-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EYHALL_100_110) = JXBY_100_110(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
601+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
602+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EYHALL_100_110) =
603+
JXBY_100_110(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
604+
(physicalconstants::MU_0 * hallRhoq);
590605
hallRhoq = FOURTH * (
591-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
592-
momentsGrid.get(i-1,j ,k )->at(fsgrids::moments::RHOQ) +
593-
momentsGrid.get(i ,j ,k+1)->at(fsgrids::moments::RHOQ) +
594-
momentsGrid.get(i-1,j ,k+1)->at(fsgrids::moments::RHOQ)
606+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
607+
+ momentsGrid.get(i - 1, j, k)->at(fsgrids::moments::RHOQ)
608+
+ momentsGrid.get(i, j, k + 1)->at(fsgrids::moments::RHOQ)
609+
+ momentsGrid.get(i - 1, j, k + 1)->at(fsgrids::moments::RHOQ)
595610
);
596-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
597-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EYHALL_001_011) = JXBY_001_011(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
611+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
612+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EYHALL_001_011) =
613+
JXBY_001_011(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
614+
(physicalconstants::MU_0 * hallRhoq);
598615
hallRhoq = FOURTH * (
599-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
600-
momentsGrid.get(i+1,j ,k )->at(fsgrids::moments::RHOQ) +
601-
momentsGrid.get(i ,j ,k+1)->at(fsgrids::moments::RHOQ) +
602-
momentsGrid.get(i+1,j ,k+1)->at(fsgrids::moments::RHOQ)
616+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
617+
+ momentsGrid.get(i + 1, j, k)->at(fsgrids::moments::RHOQ)
618+
+ momentsGrid.get(i, j, k + 1)->at(fsgrids::moments::RHOQ)
619+
+ momentsGrid.get(i + 1, j, k + 1)->at(fsgrids::moments::RHOQ)
603620
);
604-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
605-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EYHALL_101_111) = JXBY_101_111(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
621+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
622+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EYHALL_101_111) =
623+
JXBY_101_111(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBZ), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
624+
(physicalconstants::MU_0 * hallRhoq);
625+
// clang-format on
606626
break;
607627

608628
default:
@@ -669,38 +689,48 @@ void calculateEdgeHallTermZComponents(
669689
break;
670690

671691
case 2:
692+
// clang-format off
672693
hallRhoq = FOURTH * (
673-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
674-
momentsGrid.get(i-1,j ,k )->at(fsgrids::moments::RHOQ) +
675-
momentsGrid.get(i ,j-1,k )->at(fsgrids::moments::RHOQ) +
676-
momentsGrid.get(i-1,j-1,k )->at(fsgrids::moments::RHOQ)
694+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
695+
+ momentsGrid.get(i - 1, j, k)->at(fsgrids::moments::RHOQ)
696+
+ momentsGrid.get(i, j - 1, k)->at(fsgrids::moments::RHOQ)
697+
+ momentsGrid.get(i - 1, j - 1, k)->at(fsgrids::moments::RHOQ)
677698
);
678-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
679-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EZHALL_000_001) = JXBZ_000_001(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
699+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
700+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EZHALL_000_001) =
701+
JXBZ_000_001(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
702+
(physicalconstants::MU_0 * hallRhoq);
680703
hallRhoq = FOURTH * (
681-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
682-
momentsGrid.get(i+1,j ,k )->at(fsgrids::moments::RHOQ) +
683-
momentsGrid.get(i ,j-1,k )->at(fsgrids::moments::RHOQ) +
684-
momentsGrid.get(i+1,j-1,k )->at(fsgrids::moments::RHOQ)
704+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
705+
+ momentsGrid.get(i + 1, j, k)->at(fsgrids::moments::RHOQ)
706+
+ momentsGrid.get(i, j - 1, k)->at(fsgrids::moments::RHOQ)
707+
+ momentsGrid.get(i + 1, j - 1, k)->at(fsgrids::moments::RHOQ)
685708
);
686-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
687-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EZHALL_100_101) = JXBZ_100_101(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
709+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
710+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EZHALL_100_101) =
711+
JXBZ_100_101(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
712+
(physicalconstants::MU_0 * hallRhoq);
688713
hallRhoq = FOURTH * (
689-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
690-
momentsGrid.get(i-1,j ,k )->at(fsgrids::moments::RHOQ) +
691-
momentsGrid.get(i ,j+1,k )->at(fsgrids::moments::RHOQ) +
692-
momentsGrid.get(i-1,j+1,k )->at(fsgrids::moments::RHOQ)
714+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
715+
+ momentsGrid.get(i - 1, j, k)->at(fsgrids::moments::RHOQ)
716+
+ momentsGrid.get(i, j + 1, k)->at(fsgrids::moments::RHOQ)
717+
+ momentsGrid.get(i - 1, j + 1, k)->at(fsgrids::moments::RHOQ)
693718
);
694-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
695-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EZHALL_010_011) = JXBZ_010_011(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
719+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
720+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EZHALL_010_011) =
721+
JXBZ_010_011(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
722+
(physicalconstants::MU_0 * hallRhoq);
696723
hallRhoq = FOURTH * (
697-
momentsGrid.get(i ,j ,k )->at(fsgrids::moments::RHOQ) +
698-
momentsGrid.get(i+1,j ,k )->at(fsgrids::moments::RHOQ) +
699-
momentsGrid.get(i ,j+1,k )->at(fsgrids::moments::RHOQ) +
700-
momentsGrid.get(i+1,j+1,k )->at(fsgrids::moments::RHOQ)
724+
momentsGrid.get(i, j, k)->at(fsgrids::moments::RHOQ)
725+
+ momentsGrid.get(i + 1, j, k)->at(fsgrids::moments::RHOQ)
726+
+ momentsGrid.get(i, j + 1, k)->at(fsgrids::moments::RHOQ)
727+
+ momentsGrid.get(i + 1, j + 1, k)->at(fsgrids::moments::RHOQ)
701728
);
702-
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq ) ? Parameters::hallMinimumRhoq : hallRhoq ;
703-
EHallGrid.get(i,j,k)->at(fsgrids::ehall::EZHALL_110_111) = JXBZ_110_111(perturbedCoefficients, BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i,j,k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) / (physicalconstants::MU_0 * hallRhoq);
729+
hallRhoq = (hallRhoq <= Parameters::hallMinimumRhoq) ? Parameters::hallMinimumRhoq : hallRhoq;
730+
EHallGrid.get(i, j, k)->at(fsgrids::ehall::EZHALL_110_111) =
731+
JXBZ_110_111(perturbedCoefficients, BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBX), BgBGrid.get(i, j, k)->at(fsgrids::bgbfield::BGBY), technicalGrid.DX, technicalGrid.DY, technicalGrid.DZ) /
732+
(physicalconstants::MU_0 * hallRhoq);
733+
// clang-format on
704734
break;
705735

706736
default:

0 commit comments

Comments
 (0)