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
|
// Copyright (C) 2017 EDF
// All Rights Reserved
// This code is published under the GNU Lesser General Public License (GNU LGPL)
#include "StOpt/regression/LocalRegression.h"
#include "StOpt/regression/meshCalculationLocalRegression.h"
using namespace std;
using namespace Eigen;
namespace StOpt
{
LocalRegression::LocalRegression(const ArrayXi &p_nbMesh,
const bool &p_bRotationAndRecale):
LocalAdaptCellRegression(p_nbMesh, p_bRotationAndRecale) {}
LocalRegression::LocalRegression(const bool &p_bZeroDate,
const ArrayXXd &p_particles,
const ArrayXi &p_nbMesh,
const bool &p_bRotationAndRecale):
LocalAdaptCellRegression(p_bZeroDate, p_particles, p_nbMesh, p_bRotationAndRecale)
{
if ((!m_bZeroDate) && (p_nbMesh.size() != 0))
{
meshCalculationLocalRegression(m_particles, m_nbMesh, m_simToCell, m_mesh, m_mesh1D);
}
else
{
m_simToCell.setConstant(0);
m_nbMeshTotal = 1;
}
}
LocalRegression::LocalRegression(const bool &p_bZeroDate,
const ArrayXi &p_nbMesh,
const Array< array< double, 2>, Dynamic, Dynamic > &p_mesh,
const vector< shared_ptr< ArrayXd > > &p_mesh1D, const ArrayXd &p_meanX,
const ArrayXd &p_etypX, const MatrixXd &p_svdMatrix, const bool &p_bRotationAndRecale) :
LocalAdaptCellRegression(p_bZeroDate, p_nbMesh, p_mesh, p_meanX, p_etypX, p_svdMatrix, p_bRotationAndRecale)
{
if ((!m_bZeroDate) && (p_mesh1D.size() > 0))
{
m_mesh1D = p_mesh1D;
}
}
LocalRegression::LocalRegression(const LocalRegression &p_object) : LocalAdaptCellRegression(p_object)
{
const vector< shared_ptr< ArrayXd > > &mesh1D = p_object.getMesh1D();
m_mesh1D.resize(mesh1D.size());
for (size_t i = 0 ; i < m_mesh1D.size(); ++i)
m_mesh1D[i] = make_shared< ArrayXd>(*mesh1D[i]);
}
int LocalRegression:: particleToMesh(const ArrayXd &p_oneParticle)const
{
int nBase = p_oneParticle.size() + 1;
int iCell = 0 ;
int idecCell = 1;
for (int id = 0 ; id < nBase - 1 ; id++)
{
int iMesh = 1 ;
while ((p_oneParticle(id) > (*m_mesh1D[id])(iMesh)) && (iMesh < m_mesh1D[id]->size() - 1)) iMesh++;
iCell += (iMesh - 1) * idecCell;
idecCell *= m_mesh1D[id]->size() - 1;
}
return iCell ;
}
}
|
