Merge pull request ITHACA-FV#623 from giovastabile/update_muq

added new MUQ implementations

Author: giovastabile

src/ITHACA_FOMPROBLEMS/Make/files

@@ -20,5 +20,6 @@ inverseLaplacianProblem/inverseLaplacianProblem.C
 inverseLaplacianProblem_CG/inverseLaplacianProblem_CG.C
 inverseLaplacianProblem_paramBC/inverseLaplacianProblem_paramBC.C
 CompressibleSteadyNS/CompressibleSteadyNS.C
+sequentialIHTP/sequentialIHTP.C
 
 LIB = $(FOAM_USER_LIBBIN)/libITHACA_FOMPROBLEMS

src/ITHACA_FOMPROBLEMS/sequentialIHTP/createFields.H

@@ -0,0 +1,49 @@
+/*---------------------------------------------------------------------------*\
+     ██╗████████╗██╗  ██╗ █████╗  ██████╗ █████╗       ███████╗██╗   ██╗
+     ██║╚══██╔══╝██║  ██║██╔══██╗██╔════╝██╔══██╗      ██╔════╝██║   ██║
+     ██║   ██║   ███████║███████║██║     ███████║█████╗█████╗  ██║   ██║
+     ██║   ██║   ██╔══██║██╔══██║██║     ██╔══██║╚════╝██╔══╝  ╚██╗ ██╔╝
+     ██║   ██║   ██║  ██║██║  ██║╚██████╗██║  ██║      ██║      ╚████╔╝
+     ╚═╝   ╚═╝   ╚═╝  ╚═╝╚═╝  ╚═╝ ╚═════╝╚═╝  ╚═╝      ╚═╝       ╚═══╝
+
+ * In real Time Highly Advanced Computational Applications for Finite Volumes
+ * Copyright (C) 2017 by the ITHACA-FV authors
+-------------------------------------------------------------------------------
+  License
+  This file is part of ITHACA-FV
+  ITHACA-FV 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.
+  ITHACA-FV 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 ITHACA-FV. If not, see <http://www.gnu.org/licenses/>.
+ \*---------------------------------------------------------------------------*/
+Info << "Create field T" << endl;
+_T = autoPtr<volScalarField>
+     (
+         new volScalarField
+         (
+             IOobject
+             (
+                 "T",
+                 runTime.timeName(),
+                 mesh,
+                 IOobject::MUST_READ,
+                 IOobject::NO_WRITE
+             ),
+             mesh
+         )
+     );
+
+
+//_DT = autoPtr<dimensionedScalar>
+//(
+//    "DT",
+//    dimensionSet(1, 1, -3, -1, 0, 0, 0),
+//    1.0
+//);
+

src/ITHACA_FOMPROBLEMS/sequentialIHTP/createRegularization.H

@@ -0,0 +1,16 @@
+IOdictionary regularizationDict
+(
+    IOobject
+    (
+        "regularizationDict",
+        runTime.constant(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::NO_WRITE
+    )
+);
+
+linSys_solver = regularizationDict.lookupOrDefault<word>("regularizationTechnique", "fullPivLU");
+TSVD_filter = regularizationDict.lookupOrDefault<int>("TSVD_filter", 0);
+Tikhonov_filter = regularizationDict.lookupOrDefault<double>("Tikhonov_filter", 0.0);
+CG_Nsteps = regularizationDict.lookupOrDefault<int>("CG_Nsteps", 0);

src/ITHACA_FOMPROBLEMS/sequentialIHTP/createThermocouples.H

@@ -0,0 +1,34 @@
+IOdictionary thermocouplesPosition  // Kabir: "Right now, thermocouplesPosition has as many thermocouples as thermocouplesDict has. It means that thermocouplesNum is equal to the number of thermocouples in thermocouplesDict."
+(    
+    IOobject
+    (
+        "thermocouplesDict",        // inside 06enKFwDF_3dIHTP.H the function named int getObservationSize() reads measurementsDict
+        runTime.constant(),         // inside 06enKFwDF_3dIHTP.H the function named Eigen::VectorXd observe(volScalarField field) reads measurementsDict
+        mesh,                       // inside 06enKFwDF_3dIHTP.H the function named Eigen::VectorXd observe(Eigen::VectorXd field) reads measurementsDict
+        IOobject::MUST_READ,        // inside 06enKFwDF_3dIHTP.H the function named void readObservationTimes() reads observationDeltaTimesteps and observationStartTimestep
+        IOobject::NO_WRITE
+    )
+);
+
+List<vector> TCpos(thermocouplesPosition.lookup("positions"));
+thermocouplesNum = TCpos.size();
+thermocouplesPosition.lookup("timeSamplesT0") >> timeSamplesT0;                        // timeSamplesT0     1.0;
+thermocouplesPosition.lookup("timeSamplesDeltaT") >> timeSamplesDeltaT;                // timeSamplesDeltaT 1.0;
+timeSamplesNum = (runTime.endTime().value() - timeSamplesT0) / timeSamplesDeltaT + 1;  // 20, This is wrong. (It should be runTime.endTime().value()/Ntimes)/observationDeltaTimesteps
+Info << "First sampling time = " << timeSamplesT0 << endl;                             // 1
+Info << "Delta time between samples = " << timeSamplesDeltaT << endl;                  // 1
+Info << "Number of samplings = " << timeSamplesNum << endl;                            // 20, This is wrong. (It should be runTime.endTime().value()/Ntimes)/observationDeltaTimesteps
+Info << "Number of samplings should be calculated by (endtime / Ntimes) / observationDeltaTimesteps, so, the above calculated one is inncorect " << endl;
+
+
+
+// timeSamplesT0     1.0;      First sampling time
+// timeSamplesDeltaT 1.0;      Delta time between samples
+
+
+// observationDeltaTimesteps 2;       inside 06enKFwDF_3dIHTP.H
+// observationStartTimestep  1;       inside 06enKFwDF_3dIHTP.H
+
+// Number of samples in time = 50     it is equal to equal measurements.cols() inside 06enKFwDF_3dIHTP.H,   Ntimes/observationDeltaTimesteps
+// Number of sample in space = 100    it is equal to equal measurements.rows() inside 06enKFwDF_3dIHTP.H,   number of thermocouples in thermocouplesDict."
+