@@ -809,84 +809,124 @@ INITIALIZE
809809 fprintf (stderr, "Error opening file\n");
810810 exit (-1);
811811 } else {
812- printf ("%s\n", "File opened...");
812+ MPI_MASTER(
813+ printf ("%s\n", "File opened...");
814+ );
813815 }
814816 /* spectral parameters */
815817 while ((fgets (line, linelength, fp) != NULL) && (strchr (line, '#') != NULL)) {
816- printf ("%s", line);
818+ MPI_MASTER(
819+ printf ("%s", line);
820+ );
817821 }
818822 kk = sscanf (line, " %le %le %le %le %le %le %le %le %le %le %le %le %le %le", ¶_sp[0], ¶_sp[1], ¶_sp[2], ¶_sp[3], ¶_sp[4], ¶_sp[5],
819823 ¶_sp[6], ¶_sp[7], ¶_sp[8], ¶_sp[9], ¶_sp[10], ¶_sp[11], ¶_sp[12], ¶_sp[13]);
820824
821825 /* a parameter of emission time distribuation */
822826 while ((fgets (line, linelength, fp) != NULL) && (strchr (line, '#') != NULL)) {
823- printf ("%s", line);
827+ MPI_MASTER(
828+ printf ("%s", line);
829+ );
824830 }
825831 kk = sscanf (line, " %le %le %le %le %le %le %le %le %le %le", ¶_a[0], ¶_a[1], ¶_a[2], ¶_a[3], ¶_a[4], ¶_a[5], ¶_a[6], ¶_a[7],
826832 ¶_a[8], ¶_a[9]);
827833 if (kk != 10) {
828- printf ("para_a: insufficient number of data entries read\n");
834+ MPI_MASTER(
835+ printf ("para_a: insufficient number of data entries read\n");
836+ );
829837 _ABORT ();
830838 }
831839
832840 /* b parameter of emission time distribuation */
833841 while ((fgets (line, linelength, fp) != NULL) && (strchr (line, '#') != NULL)) {
834- printf ("%s", line);
842+ MPI_MASTER(
843+ printf ("%s", line);
844+ );
835845 }
836846 kk = sscanf (line, " %le %le %le %le %le %le %le %le %le %le", ¶_b[0], ¶_b[1], ¶_b[2], ¶_b[3], ¶_b[4], ¶_b[5], ¶_b[6], ¶_b[7],
837847 ¶_b[8], ¶_b[9]);
838848 if (kk != 10) {
839- printf ("para_b: insufficient number of data entries read\n");
849+ MPI_MASTER(
850+ printf ("para_b: insufficient number of data entries read\n");
851+ );
840852 _ABORT ();
841853 }
842854
843855 /* R parameter of emission time distribuation */
844856 while ((fgets (line, linelength, fp) != NULL) && (strchr (line, '#') != NULL)) {
845- printf ("%s", line);
857+ MPI_MASTER(
858+ printf ("%s", line);
859+ );
846860 }
847861 kk = sscanf (line, " %le %le %le %le %le %le %le %le %le %le", ¶_R[0], ¶_R[1], ¶_R[2], ¶_R[3], ¶_R[4], ¶_R[5], ¶_R[6], ¶_R[7],
848862 ¶_R[8], ¶_R[9]);
849863 if (kk != 10) {
850- printf ("para_R: kk=%d insufficient number of data entries read\n", kk);
864+ MPI_MASTER(
865+ printf ("para_R: kk=%d insufficient number of data entries read\n", kk);
866+ );
851867 _ABORT ();
852868 }
853869
854870 /* to parameter of emission time distribuation */
855871 while ((fgets (line, linelength, fp) != NULL) && (strchr (line, '#') != NULL)) {
856- printf ("%s", line);
872+ MPI_MASTER(
873+ printf ("%s", line);
874+ );
857875 }
858876 kk = sscanf (line, " %le %le %le %le %le %le %le %le %le %le", ¶_to[0], ¶_to[1], ¶_to[2], ¶_to[3], ¶_to[4], ¶_to[5], ¶_to[6],
859877 ¶_to[7], ¶_to[8], ¶_to[9]);
860878 if (kk != 10) {
861- printf ("para_R: insufficient number of data entries read\n");
879+ MPI_MASTER(
880+ printf ("para_R: insufficient number of data entries read\n");
881+ );
862882 _ABORT ();
863883 }
864884
865885 printf ("\n");
866886 for (k = 0; k < 14; k++) {
867- printf (" para_sp(%d)=%e\n", k, para_sp[k]);
887+ MPI_MASTER(
888+ printf (" para_sp(%d)=%e\n", k, para_sp[k]);
889+ );
868890 }
869- printf ("\n");
891+ MPI_MASTER(
892+ printf ("\n");
893+ );
870894 for (k = 0; k < 10; k++) {
871- printf (" para_a(%d)=%e\n", k, para_a[k]);
895+ MPI_MASTER(
896+ printf (" para_a(%d)=%e\n", k, para_a[k]);
897+ );
872898 }
873- printf ("\n");
899+ MPI_MASTER(
900+ printf ("\n");
901+ );
874902 for (k = 0; k < 10; k++) {
875- printf (" para_b(%d)=%e\n", k, para_b[k]);
903+ MPI_MASTER(
904+ printf (" para_b(%d)=%e\n", k, para_b[k]);
905+ );
876906 }
877- printf ("\n");
907+ MPI_MASTER(
908+ printf ("\n");
909+ );
878910 for (k = 0; k < 8; k++) {
879- printf (" para_R(%d)=%e\n", k, para_R[k]);
911+ MPI_MASTER(
912+ printf (" para_R(%d)=%e\n", k, para_R[k]);
913+ );
880914 }
881- printf ("\n");
915+ MPI_MASTER(
916+ printf ("\n");
917+ );
882918 for (k = 0; k < 10; k++) {
883- printf (" para_to(%d)=%e\n", k, para_to[k]);
919+ MPI_MASTER(
920+ printf (" para_to(%d)=%e\n", k, para_to[k]);
921+ );
884922 }
885923
886924 if (sample_E == 0) {
887925 CItot = prepare_cumulative_spectral_function (csfE, Emin, Emax, para_sp);
888926 if (CItot <= 0.0) {
889- printf (" choice of Emin and Emax gives zero neutron intensity\n");
927+ MPI_MASTER(
928+ printf (" choice of Emin and Emax gives zero neutron intensity\n");
929+ );
890930 _ABORT ();
891931 }
892932 // for(i=0;i<=1000;i++){
@@ -896,19 +936,29 @@ INITIALIZE
896936
897937 /* some checks */
898938 if (dist <= 0.0) {
899- printf (" dist must be greater zero\n");
939+ MPI_MASTER(
940+ printf (" dist must be greater zero\n");
941+ );
900942 _ABORT ();
901943 } else if (focus_xw <= 0.0) {
902- printf (" focus_xw must be greater zero\n");
944+ MPI_MASTER(
945+ printf (" focus_xw must be greater zero\n");
946+ );
903947 _ABORT ();
904948 } else if (focus_yh <= 0.0) {
905- printf (" focus_yh must be greater zero\n");
949+ MPI_MASTER(
950+ printf (" focus_yh must be greater zero\n");
951+ );
906952 _ABORT ();
907953 } else if (p_power <= 0) {
908- printf (" p_power must be greater zero\n");
954+ MPI_MASTER(
955+ printf (" p_power must be greater zero\n");
956+ );
909957 _ABORT ();
910958 } else if (n_pulses <= 0) {
911- printf (" n_pulses must be greater zero\n");
959+ MPI_MASTER(
960+ printf (" n_pulses must be greater zero\n");
961+ );
912962 _ABORT ();
913963 }
914964
0 commit comments