 
 int main(int argc, char* argv[])
 {
+    cudaFree(0);
+
     // initialize timer
     clock_t t0, t1;
     double  time_in_sec, time_total;
+    double  time_mat_J;
 
     t0 = clock();
     std::string time_txt ("");
     time_total = 0.0;
+    time_mat_J = 0.0;
 
     // use spherical harmonic d function?
     const int use_5d = 1;

            }
        }





    std::cout << time_txt;
    std::cout << "Total time used " << time_total << " sec\n";

    std::cout << "Mat_J time used " << time_mat_J << " sec\n";


    //====== the end of program ========

