defect_N=0;
for(i=N-k+1; i<=N; i++) {
defect_tmp=1;
for(n=1; n<=steps; n++)
for(j=0; j<=N; j++)
defect_tmp-=P_tau_LH[n*(N+1)*(N+1)+i*(N+1)+j];
if(defect_tmp>defect_N)
defect_N=defect_tmp;
printf("%d %lf\n",i,defect_tmp);
}
printf("defect of P_tau_LH=%lf=%e defect_0=%lf, defect_N=%lf\n", max_double (defect_0,defect_N), max_double (defect_0,defect_N), defect_0, defect_N);
module Main where
import Text.Printf
idef i steps nn =
1 - sum [ n*(n+1)**2 + i*(n+1) + j | j <- [0..nn], n <- [0..steps]]
calcAndPrint i steps nn = do
let res = idef i nn steps
printf "%d %f" i res
return res
defect k nn steps =
maximum [ calcAndPrint i steps nn | i <- [nn - k+1 .. nn]]
main = do
let defect_0 = ...
let defect_N = defect <k> <nn> <steps>
printf "defect of P_tau_LH=%f=%e\n" defect_0
</steps></nn></k>
N = start = k = defect0 = calcDefectTemp defectStart i = let getInnerPairs steps N = [(i, j) | i <- [1..steps], j <- [0..N]] in let f acc (n, j) = acc (P_tau_LH !! (n*(N+1)*..... )) in foldl f defectStart (getInnerPairs steps N) calcDefectN defectTemps = maximum (0 : defectTemps) main = do let lst = map (\i -> (i, calcDefectTemp 1.0 i)) [N-k+1..N] mapM (\(idx, def) -> print $ show idx ++ "" ++ show def) lst let defectN = calcDefectN $ map (\(idx, def) -> def) lst print $ "defect of P_tau_LH=" ++ (show $ maximum [defect0, defectN]) ++ ... return ()
Find more questions by tags Haskell