def create_scheme(self):
h0 = None
hdx = None
wcsph = WCSPHScheme(['fluid'], [],
dim=2,
rho0=rho0,
c0=c0,
h0=h0,
hdx=hdx,
nu=None,
gamma=7.0,
alpha=0.0,
beta=0.0)
tvf = TVFScheme(['fluid'], [],
dim=2,
rho0=rho0,
c0=c0,
nu=None,
p0=p0,
pb=None,
h0=h0)
edac = EDACScheme(['fluid'], [],
dim=2,
rho0=rho0,
c0=c0,
nu=None,
pb=p0,
h=h0)
iisph = IISPHScheme(fluids=['fluid'],
solids=[],
dim=2,
nu=None,
rho0=rho0,
has_ghosts=True)
crksph = CRKSPHScheme(fluids=['fluid'],
dim=2,
nu=None,
rho0=rho0,
h0=h0,
c0=c0,
p0=0.0)
gtvf = GTVFScheme(fluids=['fluid'],
solids=[],
dim=2,
rho0=rho0,
c0=c0,
nu=None,
h0=None,
pref=None)
pcisph = PCISPHScheme(fluids=['fluid'], dim=2, rho0=rho0, nu=None)
sisph = SISPHScheme(fluids=['fluid'],
solids=[],
dim=2,
nu=None,
rho0=rho0,
c0=c0,
alpha=0.0,
has_ghosts=True,
pref=p0,
rho_cutoff=0.2,
internal_flow=True,
gtvf=True)
isph = ISPHScheme(fluids=['fluid'],
solids=[],
dim=2,
nu=None,
rho0=rho0,
c0=c0,
alpha=0.0)
s = SchemeChooser(default='tvf',
wcsph=wcsph,
tvf=tvf,
edac=edac,
iisph=iisph,
crksph=crksph,
gtvf=gtvf,
pcisph=pcisph,
sisph=sisph,
isph=isph)
return s