 gpu_vector<real> marker_volume;
 gpu_vector<real> grid_vel, delta_v;
 gpu_vector<real> rhs;
-gpu_vector<Mat33> marker_Fe, marker_Fe_hat, marker_Fp, marker_delta_F;
+gpu_vector<Mat33> marker_Fe, marker_Fe_hat, marker_Fp, marker_delta_F, PolarR, PolarS;
 gpu_vector<real> old_vel_node_mpm;
 gpu_vector<real> ml, mg, mg_p, ml_p;
 gpu_vector<real> dot_g_proj_norm;

    marker_Fe_hat.resize(host_settings.num_mpm_markers);
    marker_Fp.resize(host_settings.num_mpm_markers);
    marker_delta_F.resize(host_settings.num_mpm_markers);

    kInitFeFp<<<CONFIG(host_settings.num_mpm_markers)>>>(marker_Fe.data_d,   // output
                                                         marker_Fp.data_d);  // output

    kMultiplyA<<<CONFIG(size)>>>(pos.data_d,    // input
                                 input.data_d,  //
                                 old_vel_node_mpm.data_d,
                                 PolarR.data_d,         // input
                                 PolarS.data_d,         // input
                                 marker_Fe.data_d,      // input
                                 marker_volume.data_d,  // input
                                 output.data_d);

