1702 101
format(
' #calls ',a9,
' = ',i20)
1703 102
format(
' #calls ',a9,
' (N = ',i1,
') = ',i20)
1704 111
format(
' #calls ',a9,
' = ',i20,
' or ',f10.5,
' %')
1705 112
format(
' #calls ',a9,
' (N = ',i1,
') = ',i20,
' or ',f10.5,
' %')
1707 if (.not.monitoring) then
1708 if (infoutlev_cll.ge.1) then
1709 write(ninfout_cll,*)
'COLLIER: CritPointsMonitor not initialized'
1710 write(ninfout_cll,*)
' no statistics available '
1713 else if(ncpout_cll.eq.closed_cll) then
1714 if (infoutlev_cll.ge.1) then
1715 write(ninfout_cll,*)
'COLLIER: Output for critical points switched off'
1720 write(ncpout_cll,90)
1721 90
format (//
' Collier: Summary of critical points:')
1723 write(ncpout_cll,100)
1724 100
format (/
' Total numbers of calls of COLLIER functions')
1726 if( pointscnta_cll.ne.0) then
1727 write(ncpout_cll,101)
'A_cll',pointscnta_cll
1729 if( pointscntb_cll.ne.0) then
1730 write(ncpout_cll,101)
'B_cll',pointscntb_cll
1732 if( pointscntc_cll.ne.0) then
1733 write(ncpout_cll,101)
'C_cll',pointscntc_cll
1735 if( pointscntd_cll.ne.0) then
1736 write(ncpout_cll,101)
'D_cll',pointscntd_cll
1738 if( pointscnte_cll.ne.0) then
1739 write(ncpout_cll,101)
'E_cll',pointscnte_cll
1741 if( pointscntf_cll.ne.0) then
1742 write(ncpout_cll,101)
'F_cll',pointscntf_cll
1744 if( pointscntg_cll.ne.0) then
1745 write(ncpout_cll,101)
'G_cll',pointscntg_cll
1748 if(pointscnttn_cll(i).ne.0) then
1749 write(ncpout_cll,102)
'TN_cll',i,pointscnttn_cll(i)
1753 if( pointscntaten_cll.ne.0) then
1754 write(ncpout_cll,101)
'Aten_cll',pointscntaten_cll
1756 if( pointscntbten_cll.ne.0) then
1757 write(ncpout_cll,101)
'Bten_cll',pointscntbten_cll
1759 if( pointscntcten_cll.ne.0) then
1760 write(ncpout_cll,101)
'Cten_cll',pointscntcten_cll
1762 if( pointscntdten_cll.ne.0) then
1763 write(ncpout_cll,101)
'Dten_cll',pointscntdten_cll
1765 if( pointscnteten_cll.ne.0) then
1766 write(ncpout_cll,101)
'Eten_cll',pointscnteten_cll
1768 if( pointscntften_cll.ne.0) then
1769 write(ncpout_cll,101)
'Ften_cll',pointscntften_cll
1771 if( pointscntgten_cll.ne.0) then
1772 write(ncpout_cll,101)
'Gten_cll',pointscntgten_cll
1775 if(pointscnttnten_cll(i).ne.0) then
1776 write(ncpout_cll,102)
'TNten_cll',i,pointscnttnten_cll(i)
1780 write(ncpout_cll,110) reqacc_coli
1781 110
format (/
' Numbers of calls of COLLIER functions'/ &
1782 ' with an estimated accuracy worse than reqacc_coli =',es11.4)
1783 if(pointscntdb_cll.ne.0.and.accpointscntdb_cll.ne.0) then
1784 write(ncpout_cll,111)
'DB_cll',accpointscntdb_cll,accpointscntdb_cll/real(pointscntdb_cll)*1d2
1786 if(pointscnta_cll.ne.0.and.accpointscnta_cll.ne.0) then
1787 write(ncpout_cll,111)
'A_cll',accpointscnta_cll,accpointscnta_cll/real(pointscnta_cll)*1d2
1789 if(pointscntb_cll.ne.0.and.accpointscntb_cll.ne.0) then
1790 write(ncpout_cll,111)
'B_cll',accpointscntb_cll,accpointscntb_cll/real(pointscntb_cll)*1d2
1792 if(pointscntc_cll.ne.0.and.accpointscntc_cll.ne.0) then
1793 write(ncpout_cll,111)
'C_cll',accpointscntc_cll,accpointscntc_cll/real(pointscntc_cll)*1d2
1795 if(pointscntd_cll.ne.0.and.accpointscntd_cll.ne.0) then
1796 write(ncpout_cll,111)
'D_cll',accpointscntd_cll,accpointscntd_cll/real(pointscntd_cll)*1d2
1798 if(pointscnte_cll.ne.0.and.accpointscnte_cll.ne.0) then
1799 write(ncpout_cll,111)
'E_cll',accpointscnte_cll,accpointscnte_cll/real(pointscnte_cll)*1d2
1801 if(pointscntf_cll.ne.0.and.accpointscntf_cll.ne.0) then
1802 write(ncpout_cll,111)
'F_cll',accpointscntf_cll,accpointscntf_cll*1d2/pointscntf_cll
1804 if(pointscntg_cll.ne.0.and.accpointscntg_cll.ne.0) then
1805 write(ncpout_cll,111)
'G_cll',accpointscntg_cll,accpointscntg_cll*1d2/pointscntg_cll
1808 if(pointscnttn_cll(i).ne.0.and.accpointscnttn_cll(i).ne.0) then
1809 write(ncpout_cll,112)
'TN_cll',i,accpointscnttn_cll(i),accpointscnttn_cll(i)*1d2/pointscnttn_cll(i)
1813 if(pointscntdbten_cll.ne.0.and.accpointscntdbten_cll.ne.0) then
1814 write(ncpout_cll,111)
'DBten_cll',accpointscntdbten_cll,accpointscntdbten_cll/real(pointscntdbten_cll)*1d2
1816 if(pointscntaten_cll.ne.0.and.accpointscntaten_cll.ne.0) then
1817 write(ncpout_cll,111)
'Aten_cll',accpointscntaten_cll,accpointscntaten_cll/real(pointscntaten_cll)*1d2
1819 if(pointscntbten_cll.ne.0.and.accpointscntbten_cll.ne.0) then
1820 write(ncpout_cll,111)
'Bten_cll',accpointscntbten_cll,accpointscntbten_cll/real(pointscntbten_cll)*1d2
1822 if(pointscntcten_cll.ne.0.and.accpointscntcten_cll.ne.0) then
1823 write(ncpout_cll,111)
'Cten_cll',accpointscntcten_cll,accpointscntcten_cll/real(pointscntcten_cll)*1d2
1825 if(pointscntdten_cll.ne.0.and.accpointscntdten_cll.ne.0) then
1826 write(ncpout_cll,111)
'Dten_cll',accpointscntdten_cll,accpointscntdten_cll/real(pointscntdten_cll)*1d2
1828 if(pointscnteten_cll.ne.0.and.accpointscnteten_cll.ne.0) then
1829 write(ncpout_cll,111)
'Eten_cll',accpointscnteten_cll,accpointscnteten_cll/real(pointscnteten_cll)*1d2
1831 if(pointscntften_cll.ne.0.and.accpointscntften_cll.ne.0) then
1832 write(ncpout_cll,111)
'Ften_cll',accpointscntften_cll,accpointscntften_cll*1d2/pointscntften_cll
1834 if(pointscntgten_cll.ne.0.and.accpointscntgten_cll.ne.0) then
1835 write(ncpout_cll,111)
'Gten_cll',accpointscntgten_cll,accpointscntgten_cll*1d2/pointscntgten_cll
1838 if(pointscnttnten_cll(i).ne.0.and.accpointscnttnten_cll(i).ne.0) then
1839 write(ncpout_cll,112)
'TNten_cll',i,accpointscnttnten_cll(i),accpointscnttnten_cll(i)*1d2/pointscnttnten_cll(i)
1844 130
format (/
' Numbers of calls of COLLIER functions'/ &
1845 ' with an estimated accuracy worse than '/ &
1846 ' sqrt(reqacc_coli) =',es11.4)
1849 write(ncpout_cll,120) critacc_cll
1850 120
format (/
' Numbers of calls of COLLIER functions'/ &
1851 ' with an estimated accuracy worse than critacc_coli =',es11.4)
1853 if(pointscntdb_cll.ne.0.and.critpointscntdb_cll.ne.0) then
1854 write(ncpout_cll,111)
'DB_cll',critpointscntdb_cll,critpointscntdb_cll/real(pointscntdb_cll)*1d2
1856 if(pointscnta_cll.ne.0.and.critpointscnta_cll.ne.0) then
1857 write(ncpout_cll,111)
'A_cll',critpointscnta_cll,critpointscnta_cll/real(pointscnta_cll)*1d2
1859 if(pointscntb_cll.ne.0.and.critpointscntb_cll.ne.0) then
1860 write(ncpout_cll,111)
'B_cll',critpointscntb_cll,critpointscntb_cll/real(pointscntb_cll)*1d2
1862 if(pointscntc_cll.ne.0.and.critpointscntc_cll.ne.0) then
1863 write(ncpout_cll,111)
'C_cll',critpointscntc_cll,critpointscntc_cll/real(pointscntc_cll)*1d2
1865 if(pointscntd_cll.ne.0.and.critpointscntd_cll.ne.0) then
1866 write(ncpout_cll,111)
'D_cll',critpointscntd_cll,critpointscntd_cll/real(pointscntd_cll)*1d2
1868 if(pointscnte_cll.ne.0.and.critpointscnte_cll.ne.0) then
1869 write(ncpout_cll,111)
'E_cll',critpointscnte_cll,critpointscnte_cll/real(pointscnte_cll)*1d2
1871 if(pointscntf_cll.ne.0.and.critpointscntf_cll.ne.0) then
1872 write(ncpout_cll,111)
'F_cll',critpointscntf_cll,critpointscntf_cll*1d2/pointscntf_cll
1874 if(pointscntg_cll.ne.0.and.critpointscntg_cll.ne.0) then
1875 write(ncpout_cll,111)
'G_cll',critpointscntg_cll,critpointscntg_cll*1d2/pointscntg_cll
1878 if(pointscnttn_cll(i).ne.0.and.critpointscnttn_cll(i).ne.0) then
1879 write(ncpout_cll,112)
'TN_cll',i,critpointscnttn_cll(i),critpointscnttn_cll(i)*1d2/pointscnttn_cll(i)
1883 if(pointscntdbten_cll.ne.0.and.critpointscntdbten_cll.ne.0) then
1884 write(ncpout_cll,111)
'DBten_cll',critpointscntdbten_cll,critpointscntdbten_cll/real(pointscntdbten_cll)*1d2
1886 if(pointscntaten_cll.ne.0.and.critpointscntaten_cll.ne.0) then
1887 write(ncpout_cll,111)
'Aten_cll',critpointscntaten_cll,critpointscntaten_cll/real(pointscntaten_cll)*1d2
1889 if(pointscntbten_cll.ne.0.and.critpointscntbten_cll.ne.0) then
1890 write(ncpout_cll,111)
'Bten_cll',critpointscntbten_cll,critpointscntbten_cll/real(pointscntbten_cll)*1d2
1892 if(pointscntcten_cll.ne.0.and.critpointscntcten_cll.ne.0) then
1893 write(ncpout_cll,111)
'Cten_cll',critpointscntcten_cll,critpointscntcten_cll/real(pointscntcten_cll)*1d2
1895 if(pointscntdten_cll.ne.0.and.critpointscntdten_cll.ne.0) then
1896 write(ncpout_cll,111)
'Dten_cll',critpointscntdten_cll,critpointscntdten_cll/real(pointscntdten_cll)*1d2
1898 if(pointscnteten_cll.ne.0.and.critpointscnteten_cll.ne.0) then
1899 write(ncpout_cll,111)
'Eten_cll',critpointscnteten_cll,critpointscnteten_cll/real(pointscnteten_cll)*1d2
1901 if(pointscntften_cll.ne.0.and.critpointscntften_cll.ne.0) then
1902 write(ncpout_cll,111)
'Ften_cll',critpointscntften_cll,critpointscntften_cll*1d2/pointscntften_cll
1904 if(pointscntgten_cll.ne.0.and.critpointscntgten_cll.ne.0) then
1905 write(ncpout_cll,111)
'Gten_cll',critpointscntgten_cll,critpointscntgten_cll*1d2/pointscntgten_cll
1908 if(pointscnttnten_cll(i).ne.0.and.critpointscnttnten_cll(i).ne.0) then
1909 write(ncpout_cll,112)
'TNten_cll',i,critpointscnttnten_cll(i),critpointscnttnten_cll(i)*1d2/pointscnttnten_cll(i)
1915 if (qopened_check) then
1916 write(ncheckout_cll,290)
1917 290
format (//
' Collier: Summary of COLI/DD use in mode 3:')
1919 write(ncheckout_cll,300)
1920 300
format (/
' Total numbers of uses of COLI functions')
1922 if( pointscntdb_coli.ne.0) then
1923 write(ncheckout_cll,111)
'DB_cll',pointscntdb_coli, &
1924 pointscntdb_coli*1d2/(pointscntdb_coli+pointscntdb_dd)
1926 if( pointscnta_coli.ne.0) then
1927 write(ncheckout_cll,111)
'A_cll',pointscnta_coli, &
1928 pointscnta_coli*1d2/(pointscnta_coli+pointscnta_dd)
1930 if( pointscntb_coli.ne.0) then
1931 write(ncheckout_cll,111)
'B_cll',pointscntb_coli, &
1932 pointscntb_coli*1d2/(pointscntb_coli+pointscntb_dd)
1934 if( pointscntc_coli.ne.0) then
1935 write(ncheckout_cll,111)
'C_cll',pointscntc_coli, &
1936 pointscntc_coli*1d2/(pointscntc_coli+pointscntc_dd)
1938 if( pointscntd_coli.ne.0) then
1939 write(ncheckout_cll,111)
'D_cll',pointscntd_coli, &
1940 pointscntd_coli*1d2/(pointscntd_coli+pointscntd_dd)
1942 if( pointscnte_coli.ne.0) then
1943 write(ncheckout_cll,111)
'E_cll',pointscnte_coli, &
1944 pointscnte_coli*1d2/(pointscnte_coli+pointscnte_dd)
1946 if( pointscntf_coli.ne.0) then
1947 write(ncheckout_cll,111)
'F_cll',pointscntf_coli, &
1948 pointscntf_coli*1d2/(pointscntf_coli+pointscntf_dd)
1950 if( pointscntg_coli.ne.0) then
1951 write(ncheckout_cll,111)
'G_cll',pointscntg_coli, &
1952 pointscntg_coli*1d2/(pointscntg_coli+pointscntg_dd)
1955 if(pointscnttn_coli(i).ne.0) then
1956 write(ncheckout_cll,112)
'TN_cll',i,pointscnttn_coli(i), &
1957 pointscnttn_coli(i)*1d2/(pointscnttn_coli(i)+pointscnttn_dd(i))
1961 if( pointscntaten_coli.ne.0) then
1962 write(ncheckout_cll,111)
'Aten_cll',pointscntaten_coli, &
1963 pointscntaten_coli*1d2/(pointscntaten_coli+pointscntaten_dd)
1965 if( pointscntbten_coli.ne.0) then
1966 write(ncheckout_cll,111)
'Bten_cll',pointscntbten_coli, &
1967 pointscntbten_coli*1d2/(pointscntbten_coli+pointscntbten_dd)
1969 if( pointscntcten_coli.ne.0) then
1970 write(ncheckout_cll,111)
'Cten_cll',pointscntcten_coli, &
1971 pointscntcten_coli*1d2/(pointscntcten_coli+pointscntcten_dd)
1973 if( pointscntdten_coli.ne.0) then
1974 write(ncheckout_cll,111)
'Dten_cll',pointscntdten_coli, &
1975 pointscntdten_coli*1d2/(pointscntdten_coli+pointscntdten_dd)
1977 if( pointscnteten_coli.ne.0) then
1978 write(ncheckout_cll,111)
'Eten_cll',pointscnteten_coli, &
1979 pointscnteten_coli*1d2/(pointscnteten_coli+pointscnteten_dd)
1981 if( pointscntften_coli.ne.0) then
1982 write(ncheckout_cll,111)
'Ften_cll',pointscntften_coli, &
1983 pointscntften_coli*1d2/(pointscntften_coli+pointscntften_dd)
1985 if( pointscntgten_coli.ne.0) then
1986 write(ncheckout_cll,111)
'Gten_cll',pointscntgten_coli, &
1987 pointscntgten_coli*1d2/(pointscntgten_coli+pointscntgten_dd)
1990 if(pointscnttnten_coli(i).ne.0) then
1991 write(ncheckout_cll,112)
'TNten_cll',i,pointscnttnten_coli(i), &
1992 pointscnttnten_coli(i)*1d2/(pointscnttnten_coli(i)+pointscnttnten_dd(i))
1996 write(ncheckout_cll,210)
1997 210
format (/
' Total numbers of uses of DD functions')
1999 if( pointscntdb_dd.ne.0) then
2000 write(ncheckout_cll,111)
'DB_cll',pointscntdb_dd, &
2001 pointscntdb_dd*1d2/(pointscntdb_coli+pointscntdb_dd)
2003 if( pointscnta_dd.ne.0) then
2004 write(ncheckout_cll,111)
'A_cll',pointscnta_dd, &
2005 pointscnta_dd*1d2/(pointscnta_coli+pointscnta_dd)
2007 if( pointscntb_dd.ne.0) then
2008 write(ncheckout_cll,111)
'B_cll',pointscntb_dd, &
2009 pointscntb_dd*1d2/(pointscntb_coli+pointscntb_dd)
2011 if( pointscntc_dd.ne.0) then
2012 write(ncheckout_cll,111)
'C_cll',pointscntc_dd, &
2013 pointscntc_dd*1d2/(pointscntc_coli+pointscntc_dd)
2015 if( pointscntd_dd.ne.0) then
2016 write(ncheckout_cll,111)
'D_cll',pointscntd_dd, &
2017 pointscntd_dd*1d2/(pointscntd_coli+pointscntd_dd)
2019 if( pointscnte_dd.ne.0) then
2020 write(ncheckout_cll,111)
'E_cll',pointscnte_dd, &
2021 pointscnte_dd*1d2/(pointscnte_coli+pointscnte_dd)
2023 if( pointscntf_dd.ne.0) then
2024 write(ncheckout_cll,111)
'F_cll',pointscntf_dd, &
2025 pointscntf_dd*1d2/(pointscntf_coli+pointscntf_dd)
2027 if( pointscntg_dd.ne.0) then
2028 write(ncheckout_cll,111)
'G_cll',pointscntg_dd, &
2029 pointscntg_dd*1d2/(pointscntg_coli+pointscntg_dd)
2032 if(pointscnttn_dd(i).ne.0) then
2033 write(ncheckout_cll,112)
'TN_cll',i,pointscnttn_dd(i), &
2034 pointscnttn_dd(i)*1d2/(pointscnttn_coli(i)+pointscnttn_dd(i))
2038 if( pointscntaten_dd.ne.0) then
2039 write(ncheckout_cll,111)
'Aten_cll',pointscntaten_dd, &
2040 pointscntaten_dd*1d2/(pointscntaten_coli+pointscntaten_dd)
2042 if( pointscntbten_dd.ne.0) then
2043 write(ncheckout_cll,111)
'Bten_cll',pointscntbten_dd, &
2044 pointscntbten_dd*1d2/(pointscntbten_coli+pointscntbten_dd)
2046 if( pointscntcten_dd.ne.0) then
2047 write(ncheckout_cll,111)
'Cten_cll',pointscntcten_dd, &
2048 pointscntcten_dd*1d2/(pointscntcten_coli+pointscntcten_dd)
2050 if( pointscntdten_dd.ne.0) then
2051 write(ncheckout_cll,111)
'Dten_cll',pointscntdten_dd, &
2052 pointscntdten_dd*1d2/(pointscntdten_coli+pointscntdten_dd)
2054 if( pointscnteten_dd.ne.0) then
2055 write(ncheckout_cll,111)
'Eten_cll',pointscnteten_dd, &
2056 pointscnteten_dd*1d2/(pointscnteten_coli+pointscnteten_dd)
2058 if( pointscntften_dd.ne.0) then
2059 write(ncheckout_cll,111)
'Ften_cll',pointscntften_dd, &
2060 pointscntften_dd*1d2/(pointscntften_coli+pointscntften_dd)
2062 if( pointscntgten_dd.ne.0) then
2063 write(ncheckout_cll,111)
'Gten_cll',pointscntgten_dd, &
2064 pointscntgten_dd*1d2/(pointscntgten_coli+pointscntgten_dd)
2067 if(pointscnttnten_dd(i).ne.0) then
2068 write(ncheckout_cll,112)
'TNten_cll',i,pointscnttnten_dd(i), &
2069 pointscnttnten_dd(i)*1d2/(pointscnttnten_coli(i)+pointscnttnten_dd(i))
2073 write(ncheckout_cll,220)
2074 220
format (/
' Total numbers of calls of COLI/DD functions')
2076 if( pointscntdb_coli+pointscntdb_dd.ne.0) then
2077 write(ncheckout_cll,101)
'DB_cll',pointscntdb_coli+pointscntdb_dd
2079 if( pointscnta_coli+pointscnta_dd.ne.0) then
2080 write(ncheckout_cll,101)
'A_cll',pointscnta_coli+pointscnta_dd
2082 if( pointscntb_coli+pointscntb_dd.ne.0) then
2083 write(ncheckout_cll,101)
'B_cll',pointscntb_coli+pointscntb_dd
2085 if( pointscntc_coli+pointscntc_dd.ne.0) then
2086 write(ncheckout_cll,101)
'C_cll',pointscntc_coli+pointscntc_dd
2088 if( pointscntd_coli+pointscntd_dd.ne.0) then
2089 write(ncheckout_cll,101)
'D_cll',pointscntd_coli+pointscntd_dd
2091 if( pointscnte_coli+pointscnte_dd.ne.0) then
2092 write(ncheckout_cll,101)
'E_cll',pointscnte_coli+pointscnte_dd
2094 if( pointscntf_coli+pointscntf_dd.ne.0) then
2095 write(ncheckout_cll,101)
'F_cll',pointscntf_coli+pointscntf_dd
2097 if( pointscntg_coli+pointscntg_dd.ne.0) then
2098 write(ncheckout_cll,101)
'G_cll',pointscntg_coli+pointscntg_dd
2101 if(pointscnttn_coli(i)+pointscnttn_dd(i).ne.0) then
2102 write(ncheckout_cll,102)
'TN_cll',i,pointscnttn_coli(i)+pointscnttn_dd(i)
2106 if( pointscntdbten_coli+pointscntdbten_dd.ne.0) then
2107 write(ncheckout_cll,101)
'DBten_cll',pointscntdbten_coli+pointscntdbten_dd
2109 if( pointscntaten_coli+pointscntaten_dd.ne.0) then
2110 write(ncheckout_cll,101)
'Aten_cll',pointscntaten_coli+pointscntaten_dd
2112 if( pointscntbten_coli+pointscntbten_dd.ne.0) then
2113 write(ncheckout_cll,101)
'Bten_cll',pointscntbten_coli+pointscntbten_dd
2115 if( pointscntcten_coli+pointscntcten_dd.ne.0) then
2116 write(ncheckout_cll,101)
'Cten_cll',pointscntcten_coli+pointscntcten_dd
2118 if( pointscntdten_coli+pointscntdten_dd.ne.0) then
2119 write(ncheckout_cll,101)
'Dten_cll',pointscntdten_coli+pointscntdten_dd
2121 if( pointscnteten_coli+pointscnteten_dd.ne.0) then
2122 write(ncheckout_cll,101)
'Eten_cll',pointscnteten_coli+pointscnteten_dd
2124 if( pointscntften_coli+pointscntften_dd.ne.0) then
2125 write(ncheckout_cll,101)
'Ften_cll',pointscntften_coli+pointscntften_dd
2127 if( pointscntgten_coli+pointscntgten_dd.ne.0) then
2128 write(ncheckout_cll,101)
'Gten_cll',pointscntgten_coli+pointscntgten_dd
2131 if(pointscnttnten_coli(i)+pointscnttnten_dd(i).ne.0) then
2132 write(ncheckout_cll,102)
'TNten_cll',i,pointscnttnten_coli(i)+pointscnttnten_dd(i)
2136 write(ncheckout_cll,320)
2137 320
format (/
' Numbers of comparisons of functions between COLI and DD')
2139 if(checkcntdb_cll.ne.0) then
2140 write(ncheckout_cll,101)
'DB_cll',checkcntdb_cll
2143 if(checkcnt_cll(i).ne.0) then
2144 write(ncheckout_cll,102)
'TN_cll',i,checkcnt_cll(i)
2148 if(checkcntten_cll(i).ne.0) then
2149 write(ncheckout_cll,102)
'TNten_cll',i,checkcntten_cll(i)
2153 write(ncheckout_cll,330) checkacc_cll
2154 330
format (/
' Numbers of calls of COLLIER functions'/ &
2155 ' with difference between COLI and DD > checkacc_cll =',es11.4)
2157 if(diffcntdb_cll.ne.0) then
2158 write(ncheckout_cll,111)
'DB_cll',diffcntdb_cll, &
2159 diffcntdb_cll*1d2/checkcntdb_cll
2162 if(diffcnt_cll(i).ne.0) then
2163 write(ncheckout_cll,112)
'TN_cll',i,diffcnt_cll(i), &
2164 diffcnt_cll(i)*1d2/checkcnt_cll(i)
2167 if(diffcntec_cll.ne.0) then
2168 write(ncheckout_cll,111)
'E_cll',diffcntec_cll, &
2169 diffcntec_cll*1d2/checkcnt_cll(5)
2172 if(diffcntten_cll(i).ne.0) then
2173 write(ncheckout_cll,112)
'TNten_cll',i,diffcntten_cll(i), &
2174 diffcntten_cll(i)*1d2/checkcntten_cll(i)
2180 erreventcnt(1) = erreventcnt(1) + 1
2181 erreventcnt(errflag_cll) = erreventcnt(errflag_cll) + 1
2182 acceventcnt(1) = acceventcnt(1) + 1
2183 acceventcnt(accflag_cll) = acceventcnt(accflag_cll) + 1
2185 401
format(
' #calls all ',
' = ',i20)
2186 411
format(
' #calls with errors of level ',i3,
' = ',i20,
' or ',f10.5,
' %')
2187 421
format(
' #events all ',
' = ',i20)
2188 431
format(
' #events with errors of level ',i3,
' = ',i20,
' or ',f10.5,
' %')
2190 write(nerrout_cll,400)
2191 400
format (/
' Numbers of errors in COLI functions')
2193 write(nerrout_cll,411) i,errcntcoli(i),errcntcoli(i)/real(errcnt(1))*1d2
2195 write(nerrout_cll,401) errcnt(1)
2197 write(nerrout_cll,405)
2198 405
format (/
' Numbers of errors in DD functions')
2200 write(nerrout_cll,411) i,errcntdd(i),errcntdd(i)/real(errcnt(1))*1d2
2202 write(nerrout_cll,401) errcnt(1)
2204 write(nerrout_cll,410)
2205 410
format (/
' Numbers of errors in COLLIER functions')
2207 write(nerrout_cll,411) i,errcnt(i),errcnt(i)/real(errcnt(1))*1d2
2209 write(nerrout_cll,401) errcnt(1)
2211 write(nerrout_cll,415)
2212 415
format (/
' Numbers of errors in Events')
2214 write(nerrout_cll,431) i,erreventcnt(i),erreventcnt(i)/real(erreventcnt(1))*1d2
2216 write(nerrout_cll,421) erreventcnt(1)
2219 501
format(
' #calls all ',
' = ',i20)
2220 511
format(
' #calls with accuracy of level ',i3,
' = ',i20,
' or ',f10.5,
' %')
2221 521
format(
' #events all ',
' = ',i20)
2222 531
format(
' #events with accuracy of level ',i3,
' = ',i20,
' or ',f10.5,
' %')
2224 write(ncpout_cll,510)
2225 510
format (/
' Numbers of COLLIER calls with accuracy levels')
2227 write(ncpout_cll,511) i,acccnt(i),acccnt(i)/real(acccnt(1))*1d2
2229 write(ncpout_cll,501) acccnt(1)
2231 write(ncpout_cll,500)
2232 500
format (/
' Numbers of Events with accuracy levels')
2234 write(ncpout_cll,531) i,acceventcnt(i),acceventcnt(i)/real(acceventcnt(1))*1d2
2236 write(ncpout_cll,521) acceventcnt(1)