portelli/LatAnalyze
1
0
Fork 0
mirror of https://github.com/aportelli/LatAnalyze.git synced 2025-04-11 03:20:46 +01:00
Code Releases Activity

FitInterface: methods to set global correlations

Browse Source
This commit is contained in:
Antonin Portelli 2016-04-05 12:50:02 +01:00
parent d0ef32fb9e
commit 41ea5c984e
2 changed files with 36 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
lib/FitInterface.cpp
View File

@ -285,6 +285,16 @@ void FitInterface::assumeXXCorrelated(const bool isCorr, const Index r1,
scheduleFitVarMatInit(); scheduleFitVarMatInit();
} }
void FitInterface::assumeXXCorrelated(const bool isCorr, const Index i1,
const Index i2)
{
for (Index r1 = 0; r1 < getXSize(i1); ++r1)
for (Index r2 = 0; r2 < getXSize(i2); ++r2)
{
assumeXXCorrelated(isCorr, r1, i1, r2, i2);
}
}
void FitInterface::assumeYYCorrelated(const bool isCorr, const Index k1, void FitInterface::assumeYYCorrelated(const bool isCorr, const Index k1,
const Index j1, const Index k2, const Index j1, const Index k2,
const Index j2) const Index j2)
@ -300,6 +310,18 @@ void FitInterface::assumeYYCorrelated(const bool isCorr, const Index k1,
scheduleFitVarMatInit(); scheduleFitVarMatInit();
} }
void FitInterface::assumeYYCorrelated(const bool isCorr, const Index j1,
const Index j2)
{
checkYDim(j1);
checkYDim(j2);
for (auto &p1: yDataIndex_[j1])
for (auto &p2: yDataIndex_[j2])
{
assumeYYCorrelated(isCorr, p1.first, j1, p2.first, j2);
}
}
void FitInterface::assumeXYCorrelated(const bool isCorr, const Index r, void FitInterface::assumeXYCorrelated(const bool isCorr, const Index r,
const Index i, const Index k, const Index i, const Index k,
const Index j) const Index j)
@ -312,6 +334,17 @@ void FitInterface::assumeXYCorrelated(const bool isCorr, const Index r,
scheduleFitVarMatInit(); scheduleFitVarMatInit();
} }
void FitInterface::assumeXYCorrelated(const bool isCorr, const Index i,
const Index j)
{
checkYDim(j);
for (Index r = 0; r < getXSize(i); ++r)
for (auto &p: yDataIndex_[j])
{
assumeXYCorrelated(isCorr, r, i, p.first, j);
}
}
// tests /////////////////////////////////////////////////////////////////////// // tests ///////////////////////////////////////////////////////////////////////
bool FitInterface::pointExists(const Index k) const bool FitInterface::pointExists(const Index k) const
{ {

3
lib/FitInterface.hpp
View File

@ -92,10 +92,13 @@ public:
void assumeXExact(const bool isExact, const Index i); void assumeXExact(const bool isExact, const Index i);
void assumeXXCorrelated(const bool isCorr, const Index r1, const Index i1, void assumeXXCorrelated(const bool isCorr, const Index r1, const Index i1,
const Index r2, const Index i2); const Index r2, const Index i2);
void assumeXXCorrelated(const bool isCorr, const Index i1, const Index i2);
void assumeYYCorrelated(const bool isCorr, const Index k1, const Index j1, void assumeYYCorrelated(const bool isCorr, const Index k1, const Index j1,
const Index k2, const Index j2); const Index k2, const Index j2);
void assumeYYCorrelated(const bool isCorr, const Index j1, const Index j2);
void assumeXYCorrelated(const bool isCorr, const Index r, const Index i, void assumeXYCorrelated(const bool isCorr, const Index r, const Index i,
const Index k, const Index j); const Index k, const Index j);
void assumeXYCorrelated(const bool isCorr, const Index i, const Index j);
// tests // tests
bool pointExists(const Index k) const; bool pointExists(const Index k) const;
bool pointExists(const Index k, const Index j) const; bool pointExists(const Index k, const Index j) const;