datamap.h

/******************************************************************************
 *   Copyright (C) 2005-2024 by the resistivity.net development team          *
 *   Carsten Rücker carsten@resistivity.net                                   *
 *                                                                            *
 *   Licensed under the Apache License, Version 2.0 (the "License");          *
 *   you may not use this file except in compliance with the License.         *
 *   You may obtain a copy of the License at                                  *
 *                                                                            *
 *       http://www.apache.org/licenses/LICENSE-2.0                           *
 *                                                                            *
 *   Unless required by applicable law or agreed to in writing, software      *
 *   distributed under the License is distributed on an "AS IS" BASIS,        *
 *   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. *
 *   See the License for the specific language governing permissions and      *
 *   limitations under the License.                                           *
 *                                                                            *
 ******************************************************************************/
 
#ifndef _BERT_DATAMAP__H
#define _BERT_DATAMAP__H
 
#include "bert.h"
 
#include <matrix.h>
#include <pos.h>
 
namespace GIMLI{

class DLLEXPORT DataMap {
public:
    DataMap();

    DataMap(const std::string & filename);

    DataMap(const DataMap & map);

    DataMap & operator = (const DataMap & map);

    bool isComplex() const { return complex_; }

    int save(const std::string & filename);

    int load(const std::string & filename);

    void collect(const std::vector < ElectrodeShape * > & electrodes,
                 const RMatrix & sol, bool isCEM=false);

    RVector data(const DataContainerERT & dat, bool reciprocity=false, bool imag=false);

    inline void setElectrodes(const std::vector < RVector3 > & elecs){ elecs_ = elecs; }

    inline const std::vector < RVector3 > & electrodes() const { return elecs_; }

    inline const RMatrix & map() const { return map_; }
 
protected:
    void copy_(const DataMap & map);

    std::vector < RVector3 > elecs_;

    RMatrix map_;
 
    bool complex_;
};
 
} //namespace BERT
 
#endif // _BERT_DATAMAP__H