def evaluate(self, vol_id, edge_id):
"""Calculate reflections (reverse outward momentum).
vol_id
edge_id
"""
t = self.domain.get_time()
value = self.function(t)
try:
stage = float(value)
except:
stage = float(value[0])
q = self.conserved_quantities
#q[0] = self.stage[vol_id, edge_id]
q[0] = stage
q[1] = self.xmom[vol_id, edge_id]
q[2] = self.ymom[vol_id, edge_id]
normal = self.normals[vol_id, 2 * edge_id:2 * edge_id + 2]
r = rotate(q, normal, direction=1)
r[1] = -r[1]
q = rotate(r, normal, direction=-1)
return q
def evaluate(self, vol_id, edge_id):
"""Calculate reflections (reverse outward momentum).
vol_id
edge_id
"""
t = self.domain.get_time()
value = self.function(t)
try:
stage = float(value)
except:
stage = float(value[0])
q = self.conserved_quantities
#q[0] = self.stage[vol_id, edge_id]
q[0] = stage
q[1] = self.xmom[vol_id, edge_id]
q[2] = self.ymom[vol_id, edge_id]
normal = self.normals[vol_id, 2*edge_id:2*edge_id+2]
r = rotate(q, normal, direction = 1)
r[1] = -r[1]
q = rotate(r, normal, direction = -1)
return q
def evaluate(self, vol_id, edge_id):
"""Calculate reflections (reverse outward momentum).
vol_id
edge_id
"""
q = self.conserved_quantities
q[0] = self.stage[vol_id, edge_id]
q[1] = self.xmom[vol_id, edge_id]
q[2] = self.ymom[vol_id, edge_id]
normal = self.normals[vol_id, 2 * edge_id:2 * edge_id + 2]
r = rotate(q, normal, direction=1)
r[1] = -r[1]
q = rotate(r, normal, direction=-1)
return q
def evaluate(self, vol_id, edge_id):
"""Calculate reflections (reverse outward momentum).
vol_id
edge_id
"""
q = self.conserved_quantities
q[0] = self.stage[vol_id, edge_id]
q[1] = self.xmom[vol_id, edge_id]
q[2] = self.ymom[vol_id, edge_id]
normal = self.normals[vol_id, 2*edge_id:2*edge_id+2]
r = rotate(q, normal, direction = 1)
r[1] = -r[1]
q = rotate(r, normal, direction = -1)
return q
