Merge pull request #250 from elieberman3/ejl-fft-viz-fix-1-25

EVPFFT viz fix for non-cube voxels
lanl · Jan 15, 2025 · 83959e5 · 83959e5
2 parents 4e2c0c3 + cd3f1aa
commit 83959e5
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/src/EVPFFT/src/write_micro_state.cpp b/src/EVPFFT/src/write_micro_state.cpp
@@ -505,7 +505,7 @@ void EVPFFT::write_micro_state_pvtu()
           for (int jj = 1; jj <= 3; jj++) {
             dum += defgradavg(ii,jj)*xtmp[jj-1];
           }
-          xintp(ii,kx,ky,kz) = dum;
+          xintp(ii,kx,ky,kz) = dum*delt(ii);
         }
   }); // end FOR_ALL_CLASS
   Kokkos::fence();

diff --git a/src/LS-EVPFFT/src/evpfft.cpp b/src/LS-EVPFFT/src/evpfft.cpp
@@ -424,7 +424,7 @@ void EVPFFT::init_defgrad() {
       for (int ii = 1; ii <= 3; ii++) {
         defgradavg(jj,ii) = 0.0;
       }
-      defgradavg(jj,jj) = delt(jj);
+      defgradavg(jj,jj) = 1.0;
     } 
   }); // end FOR_ALL_CLASS
 

diff --git a/src/LS-EVPFFT/src/write_micro_state.cpp b/src/LS-EVPFFT/src/write_micro_state.cpp
@@ -599,7 +599,7 @@ void EVPFFT::write_micro_state_pvtu()
           for (int jj = 1; jj <= 3; jj++) {
             dum += defgradavg_dual(ii,jj)*xtmp[jj-1];
           }
-          xintp(ii,kx,ky,kz) = dum + ufintp(ii,kx,ky,kz);
+          xintp(ii,kx,ky,kz) = (dum + ufintp(ii,kx,ky,kz))*delt(ii);
         }
   }); // end FOR_ALL_CLASS
   Kokkos::fence();