Using the following: auto y = solver.particles.sliceCurrentPosition(); auto xpos_func = KOKKOS_LAMBDA( const int p ) { return y( p, 0 ); }; within a createOutputTimeSeries seems to output the reference positions.
Using the following:
auto y = solver.particles.sliceCurrentPosition();
auto xpos_func = KOKKOS_LAMBDA( const int p ) { return y( p, 0 ); };
within a createOutputTimeSeries seems to output the reference positions.