Code indexing in gitaly is broken and leads to code not being visible to the user. We work on the issue with highest priority.

Skip to content
Snippets Groups Projects

Resolve "mixed-precision"

Merged montan_v requested to merge 146-mixed-precision into master
Compare and Show latest version
14 files
+ 3
1575
Compare changes
  • Side-by-side
  • Inline
Files
14
+ 3
3
@@ -250,7 +250,7 @@ namespace ippl {
// for the P3M collision modelling method, it indicates
// the splitting between Particle-Particle interactions
// and the Particle-Mesh computations).
double alpha = 1e6;
Trhs alpha = 1e6;
// calculate square of the mesh spacing for each dimension
Vector_t hrsq(hr_m * hr_m);
@@ -264,7 +264,7 @@ namespace ippl {
const int nghost = grn_m.getNghost();
const auto& ldom = layout_mp->getLocalNDIndex();
constexpr double ke = 2.532638e8;
constexpr Trhs ke = 2.532638e8;
// Kokkos parallel for loop to find (0,0,0) point and regularize
Kokkos::parallel_for(
@@ -277,7 +277,7 @@ namespace ippl {
const bool isOrig = (ig == 0 && jg == 0 && kg == 0);
double r = Kokkos::real(Kokkos::sqrt(view(i, j, k)));
Trhs r = Kokkos::real(Kokkos::sqrt(view(i, j, k)));
view(i, j, k) = (!isOrig) * ke * (Kokkos::erf(alpha * r) / r);
});
Loading