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
59
60
61
62
63
64
65
66
67
68
69
|
#include <memory>
#include <geners/IOException.hh>
#include <geners/arrayIO.hh>
#include "StOpt/regression/ContinuationValueDetGeners.h"
#include "StOpt/core/utils/eigenGeners.h"
using namespace StOpt;
using namespace std;
bool ContinuationValueDetGeners::write(ostream &p_of, const wrapped_base &p_base,
const bool p_dumpId) const
{
// If necessary, write out the class id
const bool status = p_dumpId ? wrappedClassId().write(p_of) : true;
// Write the object data out
if (status)
{
const wrapped_type &w = dynamic_cast<const wrapped_type &>(p_base);
int isizeRows = w.getValues().rows();
int isizeCols = w.getValues().cols();
gs::write_pod(p_of, isizeRows);
gs::write_pod(p_of, isizeCols);
gs::write_pod_array(p_of, w.getValues().data(), isizeRows * isizeCols);
std::shared_ptr< StOpt::SpaceGrid > ptrGrid = w.getGrid();
bool bSharedPtr = (ptrGrid ? true : false);
gs::write_pod(p_of, bSharedPtr);
if (bSharedPtr)
gs::write_item(p_of, *w.getGrid());
}
// Return "true" on success
return status && !p_of.fail();
}
ContinuationValueDet *ContinuationValueDetGeners::read(const gs::ClassId &p_id, istream &p_in) const
{
// Validate the class id. You might want to implement
// class versioning here.
wrappedClassId().ensureSameId(p_id);
// Read in the object data
int isizeRows = 0;
gs::read_pod(p_in, &isizeRows);
int isizeCols = 0;
gs::read_pod(p_in, &isizeCols);
int isizeLoc = isizeRows * isizeCols;
Eigen::ArrayXXd values(isizeRows, isizeCols);
gs::read_pod_array(p_in, values.data(), isizeLoc);
bool bSharedPtr ;
gs::read_pod(p_in, &bSharedPtr);
CPP11_auto_ptr<StOpt::SpaceGrid> pgrid ;
if (bSharedPtr)
pgrid = gs::read_item<StOpt::SpaceGrid>(p_in);
std::shared_ptr<StOpt::SpaceGrid > pgridShared(std::move(pgrid));
// Check that the stream is in a valid state
if (p_in.fail()) throw gs::IOReadFailure("In BIO::read: input stream failure");
// Return the object
return new StOpt::ContinuationValueDet(pgridShared, values);
}
const gs::ClassId &ContinuationValueDetGeners::wrappedClassId()
{
static const gs::ClassId wrapId(gs::ClassId::makeId<wrapped_type>());
return wrapId;
}
|
