Example #1
0
 def QuadField1D(self,solution,quad):
     #diff = npl.norm(solution.q[1]-q_old[1],1)
     self._x = solution.x.centers.copy()
     quad[0,solution.frame] = solution.t
     quad[1,solution.frame] = sq(solution.q[0]*solution.q[1],solution.x.centers)
     quad[2,solution.frame] = sq(np.sqrt(solution.q[0]**2 + solution.q[1]**2),solution.x.centers)
     return quad
Example #2
0
def QuadField1D(solution, quad):
    #diff = npl.norm(solution.q[1]-q_old[1],1)
    quad[0, solution.frame] = solution.t
    quad[1, solution.frame] = sq(solution.q[0] * solution.q[1],
                                 solution.x.centers)
    quad[2, solution.frame] = sq(np.sqrt(solution.q[0]**2 + solution.q[1]**2),
                                 solution.x.centers)
    return quad
Example #3
0
 def QuadField2D(solution,quad):
     #diff = npl.norm(solution.q[1]-q_old[1],1)
     quad[0,solution.frame] = solution.t
     quad[1,solution.frame] = sq(np.sqrt(solution.q[1]*solution.q[2]**2 + (-solution.q[0]*solution.q[2])**2),solution.x.centers)
     quad[2,solution.frame] = sq(np.sqrt(solution.q[0]**2 + solution.q[1]**2 + solution.q[2]**2),solution.x.centers)    
     return quad