RadCluster
Graph-based cluster dynamics for irradiated materials
Loading...
Searching...
No Matches
Namespaces | Variables
diagnose_stall.py File Reference

Go to the source code of this file.

Namespaces

namespace  diagnose_stall
 

Variables

 diagnose_stall.MODULE_ROOT = Path(__file__).resolve().parent.parent.parent
 
 diagnose_stall.REPO_ROOT = MODULE_ROOT.parent
 
 diagnose_stall.m = importlib.import_module(f'radcluster_code.py_utils.{mod_name}')
 
 diagnose_stall.I
 
 diagnose_stall.V
 
 diagnose_stall.i_mobile
 
 diagnose_stall.v_mobile
 
 diagnose_stall.i_discrete
 
 diagnose_stall.v_discrete
 
 diagnose_stall.saved = sys.stdout, sys.stderr
 
 diagnose_stall.stdout
 
 diagnose_stall.stderr
 
 diagnose_stall.sim
 
dict diagnose_stall.overrides
 
 diagnose_stall.inp = sim.input_data
 
bool diagnose_stall.placed = False
 
 diagnose_stall.re = sim.rate_equations
 
 diagnose_stall.PM = re.n_mom
 
 diagnose_stall.flush
 
dict diagnose_stall.config
 
 diagnose_stall.results = sim.run(solver_config=config, save_output=False, timeout_s=60)
 
 diagnose_stall.t = results['t']
 
 diagnose_stall.y = results['y']
 
 diagnose_stall.i_d = re.i_discrete
 
 diagnose_stall.i_VAC = i_d + PM * re.I_bin
 
list diagnose_stall.edges = []
 
int diagnose_stall.edge = i_d + 1
 
int diagnose_stall.n_lo = edge
 
 diagnose_stall.n_hi = max(int(edge * 1.3), edge + 1)
 
 diagnose_stall.mu0 = y[i_d + PM*k, -1]
 
int diagnose_stall.mu1 = 2 else 0
 
int diagnose_stall.n_bar = mu1 / max(abs(mu0), 1e-300)
 
int diagnose_stall.bw = n_hi - n_lo
 
 diagnose_stall.S1 = sum(range(n_lo, n_hi))
 
 diagnose_stall.S2 = sum(n*n for n in range(n_lo, n_hi))
 
int diagnose_stall.det = bw * S2 - S1 * S1
 
int diagnose_stall.c_lo = 0
 
tuple diagnose_stall.phi0_lo = (S2 - S1 * n_lo) / det
 
tuple diagnose_stall.phi1_lo = (bw * n_lo - S1) / det
 
tuple diagnose_stall.phi0_hi = (S2 - S1 * (n_hi-1)) / det
 
tuple diagnose_stall.phi1_hi = (bw * (n_hi-1) - S1) / det
 
tuple diagnose_stall.c_hi = phi0_hi * mu0 + phi1_hi * mu1
 
str diagnose_stall.flag = ""
 
 diagnose_stall.v_d = re.v_discrete
 
 diagnose_stall.vac_start = i_VAC + v_d
 
list diagnose_stall.v_edges = []
 
int diagnose_stall.m_lo = edge
 
 diagnose_stall.m_hi = max(int(edge * 1.3), edge + 1)
 
 diagnose_stall.sia_neg = np.sum(y[:i_d, -1] < -1e-30)
 
 diagnose_stall.vac_neg = np.sum(y[i_VAC:i_VAC+v_d, -1] < -1e-30)
 
 diagnose_stall.c_boundary = y[i_d - 1, -1]