blob: 5770e905dfbd9f1f03381f1ab56fee42508885e3 (
plain)
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
|
#include <vector>
template <int rank, int dim> class Tensor;
template <int dim>
class Tensor<1,dim>
{
public:
explicit Tensor (const bool initialize = true);
Tensor (const Tensor<1,dim> &);
double values[(dim!=0) ? (dim) : 1];
};
template <int dim>
Tensor<1,dim>::Tensor (const Tensor<1,dim> &p)
{
for (unsigned int i=0; i<dim; ++i)
values[i] = p.values[i];
}
template <int dim>
class KellyErrorEstimator
{
struct PerThreadData
{
std::vector<std::vector<std::vector<Tensor<1,dim> > > > psi;
PerThreadData (const unsigned int n_solution_vectors,
const unsigned int n_components,
const unsigned int n_q_points);
};
};
template <int dim>
KellyErrorEstimator<dim>::PerThreadData::
PerThreadData (const unsigned int n_solution_vectors,
const unsigned int n_components,
const unsigned int n_q_points)
{
for (unsigned int i=0; i<n_solution_vectors; ++i)
for (unsigned int qp=0;qp<n_q_points;++qp)
psi[i][qp].resize(n_components);
}
template class KellyErrorEstimator<3>;
|
