def create_tools(self):
tools = []
options = self.options
if options.remesh > 0:
if options.remesh_eq == 'm4':
equations = [M4(dest='interpolate', sources=['fluid'])]
else:
equations = None
from pysph.solver.tools import SimpleRemesher
if options.scheme == 'wcsph':
props = ['u', 'v', 'au', 'av', 'ax', 'ay', 'arho']
elif options.scheme == 'tvf':
props = ['u', 'v', 'uhat', 'vhat',
'au', 'av', 'auhat', 'avhat']
elif options.scheme == 'edac':
if 'uhat' in self.particles[0].properties:
props = ['u', 'v', 'uhat', 'vhat', 'p',
'au', 'av', 'auhat', 'avhat', 'ap']
else:
props = ['u', 'v', 'p', 'au', 'av', 'ax', 'ay', 'ap']
elif options.scheme == 'iisph':
# The accelerations are not really needed since the current
# stepper is a single stage stepper.
props = ['u', 'v', 'p']
remesher = SimpleRemesher(
self, 'fluid', props=props,
freq=self.options.remesh, equations=equations
)
tools.append(remesher)
return tools
def create_tools(self):
tools = []
if self.options.remesh > 0:
from pysph.solver.tools import SimpleRemesher
remesher = SimpleRemesher(self,
'fluid',
props=['u', 'v', 'uhat', 'vhat'],
freq=self.options.remesh)
tools.append(remesher)
return tools
def create_tools(self):
tools = []
options = self.options
if options.remesh > 0:
if options.remesh_eq == 'm4':
equations = [M4(dest='interpolate', sources=['fluid'])]
else:
equations = None
from pysph.solver.tools import SimpleRemesher
if options.scheme == 'wcsph' or options.scheme == 'crksph':
props = ['u', 'v', 'au', 'av', 'ax', 'ay', 'arho']
elif options.scheme == 'pcisph':
props = ['u', 'v', 'p']
elif options.scheme == 'tvf':
props = [
'u', 'v', 'uhat', 'vhat', 'au', 'av', 'auhat', 'avhat'
]
elif options.scheme == 'edac':
if 'uhat' in self.particles[0].properties:
props = [
'u', 'v', 'uhat', 'vhat', 'p', 'au', 'av', 'auhat',
'avhat', 'ap'
]
else:
props = ['u', 'v', 'p', 'au', 'av', 'ax', 'ay', 'ap']
elif options.scheme == 'iisph' or options.scheme == 'isph':
# The accelerations are not really needed since the current
# stepper is a single stage stepper.
props = ['u', 'v', 'p']
elif options.scheme == 'gtvf':
props = [
'uhat', 'vhat', 'what', 'rho0', 'rhodiv', 'p0', 'auhat',
'avhat', 'awhat', 'arho', 'arho0'
]
remesher = SimpleRemesher(self,
'fluid',
props=props,
freq=self.options.remesh,
equations=equations)
tools.append(remesher)
if options.shift_freq > 0:
shift = ShiftPositions(self,
'fluid',
freq=self.options.shift_freq,
shift_kind=self.options.shift_kind,
correct_velocity=self.options.correct_vel,
parameter=self.options.shift_parameter)
tools.append(shift)
return tools
def create_tools(self):
'''
Set-up remesher tool
'''
tools = []
if self.remesh > 0:
equations = [M4(dest='interpolate', sources=['fluid'])]
equations = None
props = ['u', 'v', 'au', 'av', 'ax', 'ay', 'arho']
remesher = SimpleRemesher(self,
'fluid',
props=props,
freq=self.remesh,
equations=equations)
tools.append(remesher)
return tools