def test_shape(self):
grid = grid_from_sphere(2)
space_const = function_space(grid,"DP",0)
space_lin = function_space(grid,"P",1)
op = laplace_slp(space_lin,space_lin,space_const).weak_form()
assert op.shape==(space_const.global_dof_count,space_lin.global_dof_count)
def test_shape(self):
grid = grid_from_sphere(2)
space_const = function_space(grid, "DP", 0)
space_lin = function_space(grid, "P", 1)
op = laplace_slp(space_lin, space_lin, space_const).weak_form()
assert op.shape == (space_const.global_dof_count,
space_lin.global_dof_count)
def grid(self):
return grid_from_sphere(3)
def complex_operator():
grid = grid_from_sphere(3)
space = function_space(grid, "DP", 0)
return helmholtz_slp(space, space, space, 1).weak_form()
def real_operator():
grid = grid_from_sphere(3)
space = function_space(grid, "DP", 0)
return laplace_slp(space, space, space).weak_form()
def space():
grid = grid_from_sphere(3)
space = function_space(grid,"DP",0)
return space
def grid(self):
return grid_from_sphere(3)
def space():
grid = grid_from_sphere(3)
space = function_space(grid, "DP", 0)
return space
def complex_operator():
grid = grid_from_sphere(3)
space = function_space(grid,"DP",0)
return helmholtz_slp(space,space,space,1).weak_form()
def real_operator():
grid = grid_from_sphere(3)
space = function_space(grid,"DP",0)
return laplace_slp(space,space,space).weak_form()
