debug_script.m
0b5d00c4d901master

debug_script.m
View Options

			% system('cd ../results/dev/test_diag; ./helaz3_dbg; cd $HOME/HeLaZ/wk');
			% system('cd ../results/dev/test_diag; mpirun -np 2 ./helaz3_dbg 1 2 1; cd $HOME/HeLaZ/wk');
			system('cd ../results/dev/test_diag; mpirun -np 8 ./helaz3_dbg 2 2 2; cd $HOME/HeLaZ/wk');


			filename = '../results/dev/test_diag/outputs_00.h5';

			% test phi
			f_ = h5read(filename,'/data/var3d/phi/000004/');
			f_gatherv = h5read(filename,'/data/var3d/phi_gatherv/000004/');
			f_ = f_.real + 1i*f_.imaginary;
			f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
			err = sum(sum(sum(abs(f_-f_gatherv))))/sum(sum(sum(abs(f_))));
			disp(['error_phi = ',sprintf('%2.2e',err)]);

			% test Ni00
			f_ = h5read(filename,'/data/var3d/Ni00/000004/');
			f_gatherv = h5read(filename,'/data/var3d/Ni00_gatherv/000004/');
			f_ = f_.real + 1i*f_.imaginary;
			f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
			err = sum(sum(sum(abs(f_-f_gatherv))))/sum(sum(sum(abs(f_))));
			disp(['error_Ni00 = ',sprintf('%2.2e',err)]);

			% test Nipj
			f_ = h5read(filename,'/data/var5d/moments_i/000004/');
			f_gatherv = h5read(filename,'/data/var5d/moments_i_gatherv/000004/');
			f_ = f_.real + 1i*f_.imaginary;
			f_gatherv = f_gatherv.real + 1i*f_gatherv.imaginary;
			err = sum(sum(sum(sum(sum(abs(f_-f_gatherv))))))/sum(sum(sum(sum(sum(abs(f_))))));
			disp(['error_Nipj = ',sprintf('%2.2e',err)]);