vag.hpp

Go to the documentation of this file.

/*===========================================================================
//
// File: vag.hpp
//
// Created: 2012-06-08 15:46:23+0200
//
// Authors: Knut-Andreas Lie      <Knut-Andreas.Lie@sintef.no>
//          Halvor M. Nilsen      <HalvorMoll.Nilsen@sintef.no>
//          Atgeirr F. Rasmussen  <atgeirr@sintef.no>
//          Xavier Raynaud        <Xavier.Raynaud@sintef.no>
//          Bård Skaflestad       <Bard.Skaflestad@sintef.no>
//
//==========================================================================*/


/*
  Copyright 2012 SINTEF ICT, Applied Mathematics.
  Copyright 2012 Statoil ASA.

  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_VAG_HPP_HEADER
#define OPM_VAG_HPP_HEADER



#include <opm/core/grid.h>
#include <istream>
#include <ostream>
#include <string>
#include <vector>

namespace Opm
{
    struct PosStruct{
    std::vector<int> pos;
        std::vector<int> value;
    };
    struct VAG{
    int number_of_vertices;
    int number_of_volumes;
    int number_of_faces;
    int number_of_edges;
    std::vector<double> vertices;
    PosStruct volumes_to_faces;
    PosStruct volumes_to_vertices;
    PosStruct faces_to_edges;
    PosStruct faces_to_vertices;
    std::vector<int> edges;
    std::vector<int> faces_to_volumes;
    std::vector<double> material;
    };
    void readVagGrid(std::istream& is,Opm::VAG& vag_grid);
    void writeVagFormat(std::ostream& os,Opm::VAG& vag_grid);
    template <typename T>
    void readVector(std::istream& is,std::vector<T>& vec){
    using namespace std;
    for(int i=0;i< int(vec.size());++i){
        is >> vec[i];
    }
    }
    template <typename T>
    void writeVector(std::ostream& os,std::vector<T>& vec,int n){
        typedef typename std::vector<T>::size_type sz_t;

        const sz_t nn = n;

        for (sz_t i = 0; i < vec.size(); ++i) {
            os << vec[i] << (((i % nn) == 0) ? '\n' : ' ');
        }

        if ((vec.size() % nn) != 0) {
            os << '\n';
        }
    }

    void readPosStruct(std::istream& is,int n,PosStruct& pos_struct);
    void writePosStruct(std::ostream& os,PosStruct& pos_struct);

    void vagToUnstructuredGrid(Opm::VAG& vag_grid,UnstructuredGrid& grid);

    void unstructuredGridToVag(UnstructuredGrid& grid, Opm::VAG& vag_grid);
}
#endif  /* OPM_VAG_HPP_HEADER */