File: fileMatrix.cpp

package info (click to toggle)
groops 0%2Bgit20250907%2Bds-1
  • links: PTS, VCS
  • area: non-free
  • in suites: forky, sid
  • size: 11,140 kB
  • sloc: cpp: 135,607; fortran: 1,603; makefile: 20
file content (58 lines) | stat: -rw-r--r-- 1,459 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
/***********************************************/
/**
* @file fileMatrix.cpp
*
* @brief Read/write Matrix.
*
* @author Torsten Mayer-Guerr
* @date 2013-02-07
*
*/
/***********************************************/

#define DOCSTRING_FILEFORMAT_Matrix

#include "base/import.h"
#include "inputOutput/fileArchive.h"
#include "files/fileFormatRegister.h"
#include "files/fileInstrument.h"
#include "files/fileMatrix.h"

GROOPS_REGISTER_FILEFORMAT(Matrix, FILE_MATRIX_TYPE)

/***********************************************/

void writeFileMatrix(const FileName &fileName, const const_MatrixSlice &x)
{
  try
  {
    OutFileArchive file(fileName, FILE_MATRIX_TYPE, FILE_MATRIX_VERSION);
    file<<nameValue("matrix", x);
  }
  catch(std::exception &e)
  {
    GROOPS_RETHROW(e)
  }
}

/***********************************************/

void readFileMatrix(const FileName &fileName, Matrix &x)
{
  try
  {
    InFileArchive file(fileName, ""/*arbitrary type*/, std::max(FILE_MATRIX_VERSION, FILE_INSTRUMENT_VERSION));
    if(file.type().empty() || (file.type() == FILE_MATRIX_TYPE))
      file>>nameValue("matrix", x);
    else if(file.type() == FILE_INSTRUMENT_TYPE)
      x = InstrumentFile::read(fileName).matrix();
    else
      throw(Exception("file type is '"+file.type()+"' but must be '"+FILE_MATRIX_TYPE+"' or '"+FILE_INSTRUMENT_TYPE+"'"));
  }
  catch(std::exception &e)
  {
    GROOPS_RETHROW(e)
  }
}

/***********************************************/