def m3(self):
m = self.gbasis_at(self.upts)
if 'surf-flux' in self.antialias:
fp = [_proj_l2(self._iqrules[kind], self.facebases[kind])
for kind, proj, norm in self.faces]
m = m @ block_diag(fp)
return m
def m3(self):
m = self.gbasis_at(self.upts)
if 'surf-flux' in self.antialias:
fp = [_proj_l2(self._iqrules[kind],
self.facebases[kind])
for kind, proj, norm, area in self.faces]
m = np.dot(m, block_diag(fp))
return m
def m10(self):
return block_diag([self.m9]*self.ndims)
def m9(self):
return block_diag([self.m8]*self.ndims)
def m10(self):
return block_diag([self.m9] * self.ndims)
