Commit 735bf55d authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Zero-init AM matrix

parent db3242be

src/cluster/cluster.cpp

+1 −9
Original line number Diff line number Diff line
@@ -116,7 +116,7 @@ void cluster(string config_file, string data_file, string output_path) { 
    }

    C2 *c2 = new C2(nsph, max_ici, npnt, npntts);

    complex<double> **am = new complex<double>*[mxndm];

    for (int ai = 0; ai < mxndm; ai++) am[ai] = new complex<double>[mxndm];

    for (int ai = 0; ai < mxndm; ai++) am[ai] = new complex<double>[mxndm]();

    const int ndi = c4->nsph * c4->nlim;

    C9 *c9 = new C9(ndi, c4->nlem, 2 * ndi, 2 * c4->nlem);

    double *gaps = new double[nsph]();
@@ -322,17 +322,9 @@ void cluster(string config_file, string data_file, string output_path) { 
	  }

	  if (jer != 0) break;

	} // i132 loop

	//printf("INFO: initializing matrix...");

	cms(am, c1, c1ao, c4, c6);

	//printf(" done.\n");

	//ccsam = summat(am, 960, 960);

	//printf("DEBUG: after CMS CCSAM = (%lE,%lE)\n", ccsam.real(), ccsam.imag());

	int ndit = 2 * nsph * c4->nlim;

	//printf("INFO: inverting matrix...");

	lucin(am, mxndm, ndit, jer);

	//printf(" done.\n");

	//ccsam = summat(am, 960, 960);

	//printf("DEBUG: after LUCIN CCSAM = (%lE,%lE)\n", ccsam.real(), ccsam.imag());

	if (jer != 0) break; // jxi488 loop: goes to memory clean

	ztm(am, c1, c1ao, c4, c6, c9);

	if (sconf->idfc >= 0) {