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
|
/////////////////////////////////////////////////////////////
// //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing //
// http://earth.uq.edu.au/centre-geoscience-computing //
// //
// Primary Business: Brisbane, Queensland, Australia //
// Licensed under the Open Software License version 3.0 //
// http://www.apache.org/licenses/LICENSE-2.0 //
// //
/////////////////////////////////////////////////////////////
#ifndef ESYS_LSMSTRESSTENSORCOLLECTION_H
#define ESYS_LSMSTRESSTENSORCOLLECTION_H
#include "Foundation/vec3.h"
#include "Foundation/Matrix3.h"
#include "Foundation/StlIterator.h"
#include "Tools/StressCalculator/StressTensor.h"
#include <vector>
#include <map>
namespace esys
{
namespace lsm
{
template <typename TmplStressTensorCalculator>
class StressTensorCollection
{
public:
typedef TmplStressTensorCalculator StressCalculator;
typedef std::vector<StressTensor> StressTensorVector;
typedef ForwardIterator<StressTensorVector> StressTensorIterator;
StressTensorCollection(StressCalculator &stressCalculator)
: m_stressTensorVector(),
m_pStressCalculator(&stressCalculator)
{
}
template <typename TmplContactIterator>
void addContactIterator(TmplContactIterator it)
{
m_stressTensorVector.push_back(m_pStressCalculator->calculate(it));
}
template <typename TmplContactIteratorIterator>
void addContactIterators(TmplContactIteratorIterator it)
{
while (it.hasNext())
{
addContactIterator(it.next());
}
}
StressTensorIterator getIterator()
{
return StressTensorIterator(m_stressTensorVector);
}
int size() const
{
return m_stressTensorVector.size();
}
private:
StressTensorVector m_stressTensorVector;
StressCalculator *m_pStressCalculator;
};
}
}
#endif
|
