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

Use static methods to write transition matrix to file

parent 281e5c6c

src/sphere/sphere.cpp

+8 −3
Original line number Diff line number Diff line
@@ -369,11 +369,16 @@ void sphere(string config_file, string data_file, string output_path) { 
	} // i132

	if (idfc >= 0 and nsph == 1 and jxi == jwtm) {

	  // This is the condition that writes the transition matrix to output.

	  TransitionMatrix ttms(l_max, vk, exri, c1->rmi, c1->rei, sconf->get_radius(0));

	  string ttms_name = output_path + "/c_TTMS.hd5";

	  ttms.write_binary(ttms_name, "HDF5");

	  TransitionMatrix::write_binary(

					 ttms_name, l_max, vk, exri, c1->rmi, c1->rei,

					 sconf->get_radius(0), "HDF5"

					 );

	  ttms_name = output_path + "/c_TTMS";

	  ttms.write_binary(ttms_name);

	  TransitionMatrix::write_binary(

					 ttms_name, l_max, vk, exri, c1->rmi, c1->rei,

					 sconf->get_radius(0)

					 );

	}

	double cs0 = 0.25 * vk * vk * vk / half_pi;

	sscr0(tfsas, nsph, l_max, vk, exri, c1);