Commit 4267bada authored by Hsing-Yin Chang's avatar Hsing-Yin Chang
Browse files

change flambdas values to a range

parent 390428af
...@@ -348,7 +348,7 @@ def solve_ac(generation, ...@@ -348,7 +348,7 @@ def solve_ac(generation,
alambda = 1 alambda = 1
# rlambdas = Mtot/Ntot * 1e2**(np.arange(N_procs) - N_procs/2) # rlambdas = Mtot/Ntot * 1e2**(np.arange(N_procs) - N_procs/2)
rlambdas = Mtot/Ntot * 2**(np.arange(N_procs) - N_procs/2).astype(np.float) rlambdas = Mtot/Ntot * 2**(np.arange(N_procs) - N_procs/2).astype(np.float)
flambda = 1e5 * 10**(np.arange(N_procs) - N_procs//2).astype(np.float) flambdas = 1e5 * 10**(np.arange(N_procs) - N_procs//2).astype(np.float)
summary = Region((N_procs,), summary = Region((N_procs,),
{"rank": pygion.int32, "rlambda": pygion.float32, "v1": pygion.float32, "v2": pygion.float32}) {"rank": pygion.int32, "rlambda": pygion.float32, "v1": pygion.float32, "v2": pygion.float32})
...@@ -359,7 +359,7 @@ def solve_ac(generation, ...@@ -359,7 +359,7 @@ def solve_ac(generation,
solve( solve(
uregion, uregion_ups, ac, results_p[i], summary_p[i], uregion, uregion_ups, ac, results_p[i], summary_p[i],
weights, M, M_ups, Mtot, N, weights, M, M_ups, Mtot, N,
generation, i, alambda, rlambdas[i], flambda, generation, i, alambda, rlambdas[i], flambdas[i],
reciprocal_extent, use_reciprocal_symmetry, maxiter) reciprocal_extent, use_reciprocal_symmetry, maxiter)
iref = select_ac(generation, summary) iref = select_ac(generation, summary)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment