1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
|
# Copyright (C) 2020 Matthew Scroggs
#
# This file is part of FFCx. (https://www.fenicsproject.org)
#
# SPDX-License-Identifier: LGPL-3.0-or-later
import numpy as np
import ffcx
import ffcx.codegeneration.jit
import ufl
def test_finite_element(compile_args):
ufl_element = ufl.FiniteElement("Lagrange", ufl.triangle, 1)
jit_compiled_elements, module, code = ffcx.codegeneration.jit.compile_elements(
[ufl_element], cffi_extra_compile_args=compile_args)
ufcx_element, ufcx_dofmap = jit_compiled_elements[0]
assert ufcx_element.topological_dimension == 2
assert ufcx_element.geometric_dimension == 2
assert ufcx_element.space_dimension == 3
assert ufcx_element.value_rank == 0
assert ufcx_element.value_size == 1
assert ufcx_element.reference_value_rank == 0
assert ufcx_element.reference_value_size == 1
assert ufcx_element.block_size == 1
assert ufcx_element.num_sub_elements == 0
assert ufcx_dofmap.block_size == 1
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_element_support_dofs == 3
assert ufcx_dofmap.num_entity_dofs[0] == 1
assert ufcx_dofmap.num_entity_dofs[1] == 0
assert ufcx_dofmap.num_entity_dofs[2] == 0
assert ufcx_dofmap.num_entity_dofs[3] == 0
for v in range(3):
vals = np.zeros(1, dtype=np.int32)
vals_ptr = module.ffi.cast("int *", module.ffi.from_buffer(vals))
ufcx_dofmap.tabulate_entity_dofs(vals_ptr, 0, v)
assert vals[0] == v
assert ufcx_dofmap.num_sub_dofmaps == 0
def test_vector_element(compile_args):
ufl_element = ufl.VectorElement("Lagrange", ufl.triangle, 1)
jit_compiled_elements, module, code = ffcx.codegeneration.jit.compile_elements(
[ufl_element], cffi_extra_compile_args=compile_args)
ufcx_element, ufcx_dofmap = jit_compiled_elements[0]
assert ufcx_element.topological_dimension == 2
assert ufcx_element.geometric_dimension == 2
assert ufcx_element.space_dimension == 6
assert ufcx_element.value_rank == 1
assert ufcx_element.value_shape[0] == 2
assert ufcx_element.value_size == 2
assert ufcx_element.reference_value_rank == 1
assert ufcx_element.reference_value_shape[0] == 2
assert ufcx_element.reference_value_size == 2
assert ufcx_element.block_size == 2
assert ufcx_element.num_sub_elements == 2
assert ufcx_dofmap.block_size == 2
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_element_support_dofs == 3
assert ufcx_dofmap.num_entity_dofs[0] == 1
assert ufcx_dofmap.num_entity_dofs[1] == 0
assert ufcx_dofmap.num_entity_dofs[2] == 0
assert ufcx_dofmap.num_entity_dofs[3] == 0
for v in range(3):
vals = np.zeros(1, dtype=np.int32)
vals_ptr = module.ffi.cast("int *", module.ffi.from_buffer(vals))
ufcx_dofmap.tabulate_entity_dofs(vals_ptr, 0, v)
assert vals[0] == v
assert ufcx_dofmap.num_sub_dofmaps == 2
def test_tensor_element(compile_args):
ufl_element = ufl.TensorElement("Lagrange", ufl.triangle, 1)
jit_compiled_elements, module, code = ffcx.codegeneration.jit.compile_elements(
[ufl_element], cffi_extra_compile_args=compile_args)
ufcx_element, ufcx_dofmap = jit_compiled_elements[0]
assert ufcx_element.topological_dimension == 2
assert ufcx_element.geometric_dimension == 2
assert ufcx_element.space_dimension == 12
assert ufcx_element.value_rank == 2
assert ufcx_element.value_shape[0] == 2
assert ufcx_element.value_shape[1] == 2
assert ufcx_element.value_size == 4
assert ufcx_element.reference_value_rank == 2
assert ufcx_element.reference_value_shape[0] == 2
assert ufcx_element.reference_value_shape[1] == 2
assert ufcx_element.reference_value_size == 4
assert ufcx_element.block_size == 4
assert ufcx_element.num_sub_elements == 4
assert ufcx_dofmap.block_size == 4
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_global_support_dofs == 0
assert ufcx_dofmap.num_element_support_dofs == 3
assert ufcx_dofmap.num_entity_dofs[0] == 1
assert ufcx_dofmap.num_entity_dofs[1] == 0
assert ufcx_dofmap.num_entity_dofs[2] == 0
assert ufcx_dofmap.num_entity_dofs[3] == 0
for v in range(3):
vals = np.zeros(1, dtype=np.int32)
vals_ptr = module.ffi.cast("int *", module.ffi.from_buffer(vals))
ufcx_dofmap.tabulate_entity_dofs(vals_ptr, 0, v)
assert vals[0] == v
assert ufcx_dofmap.num_sub_dofmaps == 4
|
