def __init__(self, matrix, source_function=None):
self.matrix = matrix
# TODO: check that A is diagonalizable with real eigenvalues
flux_function = flux_functions.ConstantMatrix(matrix)
super().__init__(flux_function, source_function)
def __init__(self, q, matrix, source_function=None):
flux_function = flux_functions.ConstantMatrix(matrix)
app.ExactOperator.__init__(self, q, flux_function, source_function)
def __init__(self, matrix, source_function=None):
# matrix should be shape (num_eqns, num_eqns, 2)
self.matrix = matrix
flux_function = flux_functions.ConstantMatrix(matrix)
super().__init__(flux_function, source_function)