Commit acf7a405 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Create PPM matrix files when using DEBUG_AM

parent a9b5090b

src/cluster/cluster.cpp

+8 −2
Original line number Diff line number Diff line
@@ -921,6 +921,7 @@ int cluster_jxi488_cycle( 
#ifdef DEBUG_AM

    VirtualAsciiFile *outam1 = new VirtualAsciiFile();

    string outam1_name = output_path + "/c_AM1_JXI" + to_string(jxi488) + ".txt";

    string outam1_ppm_name = output_path + "/c_AM1_JXI" + to_string(jxi488) + ".ppm";

    sprintf(virtual_line, " AM matrix after CMS before LUCIN\n");

    outam1->append_line(virtual_line);

    sprintf(virtual_line, " %d\n", ndit);
@@ -930,6 +931,7 @@ int cluster_jxi488_cycle( 
    write_dcomplex_matrix(outam1, cid->am, ndit, ndit, " %5d %5d (%17.8lE,%17.8lE)\n", 1);

    outam1->write_to_disk(outam1_name);

    delete outam1;

    write_matrix_as_ppm(cid->am[0], ndit, ndit, outam1_ppm_name);

#endif // DEBUG_AM

#ifdef USE_NVTX

    nvtxRangePop();
@@ -946,6 +948,7 @@ int cluster_jxi488_cycle( 
#ifdef DEBUG_AM

    VirtualAsciiFile *outam2 = new VirtualAsciiFile();

    string outam2_name = output_path + "/c_AM2_JXI" + to_string(jxi488) + ".txt";

    string outam2_ppm_name = output_path + "/c_AM2_JXI" + to_string(jxi488) + ".ppm";

    sprintf(virtual_line, " AM matrix after LUCIN before ZTM\n");

    outam2->append_line(virtual_line);

    sprintf(virtual_line, " %d\n", ndit);
@@ -955,6 +958,7 @@ int cluster_jxi488_cycle( 
    write_dcomplex_matrix(outam2, cid->am, ndit, ndit);

    outam2->write_to_disk(outam2_name);

    delete outam2;

    write_matrix_as_ppm(cid->am[0], ndit, ndit, outam2_ppm_name);

#endif // DEBUG_AM

#ifdef USE_NVTX

    nvtxRangePop();
@@ -973,15 +977,17 @@ int cluster_jxi488_cycle( 
#ifdef DEBUG_AM

    VirtualAsciiFile *outam3 = new VirtualAsciiFile();

    string outam3_name = output_path + "/c_AM3_JXI" + to_string(jxi488) + ".txt";

    sprintf(virtual_line, " AM matrix after ZTM\n");

    string outam3_ppm_name = output_path + "/c_AM3_JXI" + to_string(jxi488) + ".txt";

    sprintf(virtual_line, " AM0M matrix after ZTM\n");

    outam3->append_line(virtual_line);

    sprintf(virtual_line, " %d\n", ndit);

    outam3->append_line(virtual_line);  

    sprintf(virtual_line, " I1+1   I2+1    Real    Imag\n");

    outam3->append_line(virtual_line);

    write_dcomplex_matrix(outam3, cid->am, ndit, ndit);

    write_dcomplex_matrix(outam3, cid->c1->am0m, 2 * cid->c1->nlem, 2 * cid->c1->nlem);

    outam3->write_to_disk(outam3_name);

    delete outam3;

    write_matrix_as_ppm(cid->c1->am0m[0], 2 * cid->c1->nlem, 2 * cid->c1->nlem, outam3_ppm_name);

#endif // DEBUG_AM

  }

  interval_start = chrono::high_resolution_clock::now();