flow 0.4.0
Kokkos cut-cell IBM incompressible Navier-Stokes solver + pnm pore extraction
Loading...
Searching...
No Matches
gp_pressure_classification.py File Reference

Go to the source code of this file.

Namespaces

namespace  gp_pressure_classification
 

Functions

 gp_pressure_classification.comps (open_faces, active)
 

Variables

 gp_pressure_classification.pack
 
 gp_pressure_classification.R
 
 gp_pressure_classification.pk = np.load(pack)
 
 gp_pressure_classification.box = np.asarray(pk["box"], float)
 
 gp_pressure_classification.G = np.round(box * R).astype(int)
 
 gp_pressure_classification.c = np.asarray(pk["centers"]) * R
 
 gp_pressure_classification.r = np.asarray(pk["scales"]) * R
 
 gp_pressure_classification.S = np.full(tuple(G), 1e30, np.float32)
 
list gp_pressure_classification.ax = [np.arange(G[k]) + 0.5 for k in range(3)]
 
 gp_pressure_classification.cs = c + sh * G
 
 gp_pressure_classification.keep = np.all((cs + (r + 3)[:, None] > 0) & (cs - (r + 3)[:, None] < G), axis=1)
 
 gp_pressure_classification.i0
 
 gp_pressure_classification.i1
 
 gp_pressure_classification.j0
 
 gp_pressure_classification.j1
 
 gp_pressure_classification.k0
 
 gp_pressure_classification.k1
 
 gp_pressure_classification.d
 
 gp_pressure_classification.out
 
 gp_pressure_classification.n = int(np.prod(G))
 
int gp_pressure_classification.fluidC = 0
 
list gp_pressure_classification.Fm = [0.5 * (np.roll(S, 1, axis=a) + S) >= 0 for a in range(3)]
 
 gp_pressure_classification.nfaces = sum(Fm[a].astype(np.int8) + np.roll(Fm[a], -1, axis=a).astype(np.int8) for a in range(3))
 
int gp_pressure_classification.hasface = nfaces > 0
 
list gp_pressure_classification.cur_open = [np.roll(fluidC, 1, axis=a) & fluidC for a in range(3)]
 
 gp_pressure_classification.nc
 
 gp_pressure_classification.tot
 
 gp_pressure_classification.orph
 
 gp_pressure_classification.nc2
 
 gp_pressure_classification.tot2
 
 gp_pressure_classification.orph2