SPILST.hpp

Go to the documentation of this file.

#ifndef SPILST_H
#define SPILST_H
#include <iostream>
#include <vector>
#include <petscksp.h>
#include <slepceps.h>
#include <slepcpep.h>
#include <string>
#include <tuple>
#include "SPIMat.hpp"
#include "SPIVec.hpp"
#include "SPIbaseflow.hpp"
#include "SPIgrid.hpp"
#include "SPIparams.hpp"
 
namespace SPI{
    std::tuple<PetscScalar, SPIVec> LST_temporal(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow,SPIVec q=SPIVec()); // temporal LST solution
    std::tuple<PetscScalar, SPIVec> LST_spatial(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow, SPIVec q=SPIVec()); // spatial LST solution
    std::tuple<PetscScalar, PetscScalar, SPIVec, SPIVec> LST_spatial_cg(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow); // spatial LST solution
    std::tuple<PetscScalar, PetscScalar, SPIVec, SPIVec> LSTNP_spatial(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow, SPIVec ql=SPIVec(),SPIVec qr=SPIVec()); // spatial non-parallel LST solution
    std::tuple<PetscScalar, SPIVec> LSTNP_spatial_right(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow, SPIVec qr=SPIVec()); // spatial non-parallel LST solution
    std::tuple<PetscScalar, SPIVec> LSTNP_spatial_right2(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow, SPIVec qr=SPIVec()); // spatial non-parallel LST solution
    std::tuple<std::vector<PetscScalar>, std::vector<SPIVec>> LSTNP_spatials_right(SPIparams &params, SPIgrid1D &grid, SPIbaseflow &baseflow, std::vector<PetscScalar> &alphas, std::vector<SPIVec> &qrs); // spatial non-parallel LST solution for multiple specific alphas
}
#endif // SPILST_H