You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
assertnp.shape(Q) ==np.shape(S), "For node type "+type[0]+", Qmat and Smat have different shape"
42
+
assertnp.shape(Q) ==np.shape(S), "For node type "+coll.__class__.__name__+", Qmat and Smat have different shape"
67
43
shape=np.shape(Q)
68
-
assertshape[0] ==shape[1], "For node type "+type[0]+", Qmat / Smat are not quadratic"
44
+
assertshape[0] ==shape[1], "For node type "+coll.__class__.__name__+", Qmat / Smat are not quadratic"
69
45
SSum=np.cumsum(S[:,:],axis=0)
70
46
foriinrange(0,M):
71
47
assertnp.linalg.norm( Q[i,:] -SSum[i,:] ) <1e-15, "For node type "+coll.__class__.__name__+", Qmat and Smat did not satisfy the expected summation property."
72
48
73
49
74
-
# TEST 3:
75
-
# Check that the partial quadrature rules from Qmat entries have order equal to number of nodes M
assertabs(int_ex-int_coll)<1e-12, "For node type "+coll.__class__.__name__+", partial quadrature from Qmat rule failed to integrate polynomial of degree M-1 exactly for M = "+str(M)
95
64
96
-
# TEST 3:
97
-
# Check that the partial quadrature rules from Smat entries have order equal to number of nodes M
0 commit comments