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
|
"This module provides simple C++ code for verification of UFC code."
# Copyright (C) 2010 Kristian B. Oelgaard
#
# This file is part of FFC.
#
# FFC is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# FFC is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with FFC. If not, see <http://www.gnu.org/licenses/>.
#
# First added: 2010-01-18
# Last changed: 2010-01-18
# Common code for all integral types
evaluate_basis_code = """\
#include <iostream>
#include <ufc.h>
#include "test.h"
int main()
{
// Create element
%(element)s element;
// Size of dof_values
// FIXME: This will not work for TensorElements
int N = element.value_dimension(0);
// Create values
double* dof_values = new double[N];
for (int i = 0; i < N; i++)
dof_values[i] = 0.0;
// Create cell and fill with some arbitrary data
double cell_coordinates[8][3] = {{0.90, 0.34, 0.45},
{0.56, 0.76, 0.83},
{0.98, 0.78, 0.19},
{0.12, 0.56, 0.66},
{0.96, 0.78, 0.63},
{0.11, 0.35, 0.49},
{0.51, 0.88, 0.65},
{0.98, 0.45, 0.01}};
ufc::cell cell;
cell.coordinates = new double * [8];
for (int i = 0; i < 8; i++)
{
cell.coordinates[i] = new double[3];
for (int j = 0; j < 3; j++)
cell.coordinates[i][j] = cell_coordinates[i][j];
}
// Random coordinates where we want to evaluate the basis functions
double coordinates[3] = {0.32, 0.51, 0.05};
// Loop element space dimension and call evaluate_basis.
for (unsigned int i = 0; i < element.space_dimension(); i++)
{
element.evaluate_basis(i, dof_values, coordinates, cell);
// Print values
for (int j = 0; j < N; j++)
std::cout << dof_values[j] << " ";
}
std::cout << std::endl;
return 0;
}
"""
|
