diff options
| author |  Axel Huebl <axel.huebl@plasma.ninja>
| 2021-11-22 18:04:18 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-11-22 18:04:18 -0800 |
| commit | 9e02b95c1f7310550d0ba3d89639b44b72a17de9 (patch) | |
| tree | 25c7bd7b477e9a502adf640fd6abaa1df70bb538 /Source/Particles | |
| parent | e55fef18dea27108247024fffba4fccd84206ac2 (diff) | |
| download | WarpX-9e02b95c1f7310550d0ba3d89639b44b72a17de9.tar.gz WarpX-9e02b95c1f7310550d0ba3d89639b44b72a17de9.tar.zst WarpX-9e02b95c1f7310550d0ba3d89639b44b72a17de9.zip | |
Python: Start 1D Support (#2582)
* Python: Start 1D Support
Start supporting 1D builds in Python.
* Fix 1D: PML FillZero unused
... since PMLs are not yet supported in 1D.
* BeamRelevant: Fix unused p_pos0
* FromTXYEFileLaserProfile: Not Impl in 1D
* QED Schwinger: 1D not Implemented
Fix unused warnings, add aborts.
* 1D RealVect/IntVect: Initialization
Use explicit scalar constructors, no braces.
Fix warning in clang 10.
* 1D NCI Filter: Fix unused members & Init
Unimplemented, but throws warnings.
* PSATD: 1D not Implemented
- remove compile warnings
- start porting some parts
* NCIGodfreyFilter: Clean up 2D & Else
Diffstat (limited to 'Source/Particles')
| -rw-r--r-- | Source/Particles/ElementaryProcess/QEDSchwingerProcess.H | 5 | ||||
| -rw-r--r-- | Source/Particles/MultiParticleContainer.cpp | 4 | ||||
| -rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 6 |
3 files changed, 10 insertions, 5 deletions
diff --git a/Source/Particles/ElementaryProcess/QEDSchwingerProcess.H b/Source/Particles/ElementaryProcess/QEDSchwingerProcess.H index c3582563b..6867544b8 100644 --- a/Source/Particles/ElementaryProcess/QEDSchwingerProcess.H +++ b/Source/Particles/ElementaryProcess/QEDSchwingerProcess.H @@ -78,7 +78,10 @@ struct SchwingerTransformFunc const amrex::Real total_weight) const noexcept { for (int n = 0; n < N; ++n){ -#if (AMREX_SPACEDIM == 3) +#if (AMREX_SPACEDIM == 1) + amrex::ignore_unused(dst1, dst2, i_dst1, i_dst2, N, total_weight); + amrex::Abort("SchwingerTransformFunc not implemented in 1D"); +#elif (AMREX_SPACEDIM == 3) dst1.m_rdata[m_weight_index][i_dst1+n] = total_weight/N; dst2.m_rdata[m_weight_index][i_dst2+n] = total_weight/N; #elif (AMREX_SPACEDIM == 2) diff --git a/Source/Particles/MultiParticleContainer.cpp b/Source/Particles/MultiParticleContainer.cpp index 70254aadd..fa0ffd2ee 100644 --- a/Source/Particles/MultiParticleContainer.cpp +++ b/Source/Particles/MultiParticleContainer.cpp @@ -1314,7 +1314,9 @@ MultiParticleContainer::doQEDSchwinger () // Get cell volume. In 2D the transverse size is // chosen by the user in the input file. amrex::Geometry const & geom = warpx.Geom(level_0); -#if (AMREX_SPACEDIM == 2) +#if (AMREX_SPACEDIM == 1) + const auto dV = geom.CellSize(0); // TODO: scale properly +#elif (AMREX_SPACEDIM == 2) const auto dV = geom.CellSize(0) * geom.CellSize(1) * m_qed_schwinger_y_size; #elif (AMREX_SPACEDIM == 3) diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index e34c4c28b..cdc1b0e02 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -1816,12 +1816,12 @@ PhysicalParticleContainer::applyNCIFilter ( const auto& nci_godfrey_filter_bxbyez = WarpX::GetInstance().nci_godfrey_filter_bxbyez; #if (AMREX_SPACEDIM == 1) - const Box& tbox = amrex::grow(box,{static_cast<int>(WarpX::noz)}); + const Box& tbox = amrex::grow(box, static_cast<int>(WarpX::noz)); #elif (AMREX_SPACEDIM == 2) - const Box& tbox = amrex::grow(box,{static_cast<int>(WarpX::nox), + const Box& tbox = amrex::grow(box, {static_cast<int>(WarpX::nox), static_cast<int>(WarpX::noz)}); #else - const Box& tbox = amrex::grow(box,{static_cast<int>(WarpX::nox), + const Box& tbox = amrex::grow(box, {static_cast<int>(WarpX::nox), static_cast<int>(WarpX::noy), static_cast<int>(WarpX::noz)}); #endif |