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)