FlowModelComponents.h

Go to the documentation of this file.

#ifndef FLOWMODELMICOMPONENTS_HPP_
#define FLOWMODELMICOMPONENTS_HPP_
 
#include "MUQ/Utilities/MultiIndices/MultiIndex.h"
 
#include "MUQ/Modeling/LinearAlgebra/EigenLinearOperator.h"
 
#include "MUQ/SamplingAlgorithms/SamplingProblem.h"
 
#include "FlowEquation.h"
 
struct Discretization {
  Discretization(std::size_t const numCells);
 
  const std::size_t numCells;
 
  const std::size_t numNodes;
 
  const Eigen::VectorXd nodeLocs;
 
  const Eigen::VectorXd cellLocs;
};
 
struct Data {
  Data(Eigen::VectorXd const& x, Eigen::VectorXd const& soln, Eigen::VectorXd const& obsLoc, Eigen::VectorXd const& obs);
 
  const Eigen::VectorXd x;
 
  const Eigen::VectorXd soln;
 
  const Eigen::VectorXd obsLoc;
 
  const Eigen::VectorXd obs;
};
 
 
std::shared_ptr<muq::Modeling::LinearOperator> InterpolateOntoObs(Discretization const& mesh, Eigen::VectorXd const& obsLoc);
 
Eigen::VectorXd GetRecharge(Discretization const& mesh);
 
Eigen::VectorXd GetTrueLogConductivity(Discretization const& mesh);
 
Data GenerateData(std::size_t const numCells, std::size_t const numObs, double const obsVar);
 
std::vector<std::shared_ptr<muq::Modeling::ModPiece> > ConstructDensities(std::size_t const numLevels, std::size_t const baseRefinement);
 
std::shared_ptr<muq::Modeling::ModPiece> ConstructDensity(std::size_t const numCells, Data const& data, double const likelihoodVar);
 
#endif