fsh.h

Go to the documentation of this file.

/*
  Copyright 2010 SINTEF ICT, Applied Mathematics.

  This file is part of the Open Porous Media project (OPM).

  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.

  OPM 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 General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef OPM_FSH_HEADER_INCLUDED
#define OPM_FSH_HEADER_INCLUDED

#include <opm/core/grid.h>
#include <opm/core/pressure/legacy_well.h>
#include <opm/core/pressure/flow_bc.h>

#ifdef __cplusplus
extern "C" {
#endif


struct CSRMatrix;
struct fsh_impl;

struct fsh_data {
    int               max_ngconn;

    size_t            sum_ngconn2;

    /* Linear system */
    struct CSRMatrix *A;        
    double           *b;        
    double           *x;        
    struct fsh_impl  *pimpl;
};



void
fsh_destroy(struct fsh_data *h);


struct fsh_data *
cfsh_construct(struct UnstructuredGrid *G, well_t *W);



void
cfsh_assemble(struct FlowBoundaryConditions *bc,
              const double    *src,
              const double    *Binv,
              const double    *Biv,
              const double    *P,
              const double    *gpress,
              well_control_t  *wctrl,
              const double    *WI,
              const double    *BivW,
              const double    *wdp,
              struct fsh_data *h);


struct fsh_data *
ifsh_construct(struct UnstructuredGrid *G, well_t *W);


void
ifsh_assemble(struct FlowBoundaryConditions *bc,
              const double    *src,
              const double    *Binv,
              const double    *gpress,
              well_control_t  *wctrl,
              const double    *WI,
              const double    *wdp,
              struct fsh_data *h);




void
fsh_press_flux(struct UnstructuredGrid *G,
               const double *Binv, const double *gpress,
               struct fsh_data *h,
               double *cpress, double *fflux,
               double *wpress, double *wflux);

#ifdef __cplusplus
}
#endif


#endif  /* OPM_FSH_HEADER_INCLUDED */