patchRecovery.H

Go to the documentation of this file.

/*******************************************************************************
* Promesh                                                                      *
* Copyright (C) 2022, IllinoisRocstar LLC. All rights reserved.                *
*                                                                              *
* Promesh is the property of IllinoisRocstar LLC.                              *
*                                                                              *
* IllinoisRocstar LLC                                                          *
* Champaign, IL                                                                *
* www.illinoisrocstar.com                                                      *
* promesh@illinoisrocstar.com                                                  *
*******************************************************************************/
/*******************************************************************************
* This file is part of Promesh                                                 *
*                                                                              *
* This version of Promesh 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.                                       *
*                                                                              *
* Promesh 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 this program. If not, see <https://www.gnu.org/licenses/>.        *
*                                                                              *
*******************************************************************************/
#ifndef NEMOSYS_PATCHRECOVERY_H_
#define NEMOSYS_PATCHRECOVERY_H_

#include "nemosys_export.h"

#include "Integration/Cubature.H"

#include <Eigen/Core>

// pair type for coordinate and data there
//using pntDataPair = std::pair<std::vector<double>, std::vector<double>>;

// holds gauss points and data at these points as pairs
//using pntDataPairVec = std::vector<pntDataPair>;

class NEMOSYS_EXPORT PatchRecovery
{
  public:
    PatchRecovery(vtkDataSet *_dataSet, int _order,
                  const std::vector<int> &arrayIDs);

    ~PatchRecovery() = default;

    // disabled copy and copy-assignment constructors
    PatchRecovery(const PatchRecovery &that) = delete;
    PatchRecovery &operator=(const PatchRecovery &that) = delete;

    // approximates nodal values with polynomial approximant over patch of node
    // sets recovered nodal values as new point data
    void recoverNodalSolution(bool ortho);
    /* 1) approximates nodal values with polynomial approximant over patch of node
       2) sets error^2 in recovered values as new point data
       3) sets RMSE integrated over each cell, normalized by cell vol as new cell data
       4) sets patch averaged element size as new point data
       returns array of integrals of errors over entire mesh */
    std::vector<std::vector<double>> computeNodalError();

  private:
    std::unique_ptr<GaussCubature> cubature;
    int order;

    // get min max of each coord (for regularization)
    std::vector<double>
    getMinMaxCoords(const std::vector<std::vector<double>> &coords) const;
    // regularize coordinates as x_new = -1 + 2*(x_old - x_min)/(x_max - x_min)
    // for proper conditioning of least squares matrix
    void regularizeCoords(std::vector<std::vector<double>> &coords,
                          std::vector<double> &genNodeCoord) const;

    // extract coordinates and data from pntDataPair
    void extractAxesAndData(const pntDataPairVec &pntsAndData,
                            std::vector<std::vector<double>> &coords,
                            std::vector<Eigen::VectorXd> &data,
                            const std::vector<int> &numComponents,
                            int &pntNum) const;
};

#endif  // NEMOSYS_PATCHRECOVERY_H_