-
Notifications
You must be signed in to change notification settings - Fork 14
/
mexutils.h
140 lines (118 loc) · 2.81 KB
/
mexutils.h
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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
//
// Written by Petter Strandmark 2010
//
// Contains a function for flushing the output buffer to
// MATLAB
//
#ifndef MEX_UTILS_PETTER
#define MEX_UTILS_PETTER
#include <map>
#include <string>
#include <stdexcept>
#include <vector>
using std::map;
using std::string;
using std::vector;
using std::runtime_error;
#include "mex.h"
#include <sstream>
#define ASSERT(cond) if (!(cond)) { std::stringstream sout; \
sout << "Error (line " << __LINE__ << "): " << #cond; \
throw runtime_error(sout.str()); }
namespace {
void flush_output()
{
//Hack
mexEvalString("drawnow");
}
}
class MexParams
{
public:
MexParams(int nrhs, const mxArray *prhs[])
{
if (nrhs < 1) {
return;
}
int i = 0;
if (mxGetClassID(prhs[i])==mxSTRUCT_CLASS) {
ASSERT(mxGetNumberOfElements(prhs[i]) == 1);
int nfield = mxGetNumberOfFields(prhs[i]);
for (int field=0;field<nfield;++field) {
const char* fieldname = mxGetFieldNameByNumber(prhs[i],field);
params[fieldname] = mxGetFieldByNumber(prhs[i],0,field);
}
i++;
}
for (;i<nrhs-1; i+=2) {
const mxArray* mxKey = prhs[i];
const mxArray* mxVal = prhs[i+1];
char buffer[1024];
if (mxGetString(mxKey, buffer, 1024)) {
throw runtime_error("Expected string parameter name");
}
params[buffer] = mxVal;
}
}
~MexParams()
{
}
template<typename T> T getdef(const string& key, T def );
template<typename T> T get(const string& key, T def=T() )
{
return getdef<T>(key,def);
}
private:
map<string,const mxArray*> params;
};
//
// Default: return a scalar
//
template<typename T> T MexParams::getdef(const string& key, T def)
{
const mxArray* mxVal = params[key];
if (!mxVal) {
//throw runtime_error("No such key: " + key);
return def;
}
ASSERT( mxGetNumberOfElements(mxVal) == 1 );
return mxGetScalar(mxVal);
}
//
// Vector of doubles
//
template<> vector<double> MexParams::getdef(const string& key, vector<double> def)
{
const mxArray* mxVal = params[key];
if (!mxVal) {
//throw runtime_error("No such key: " + key);
return def;
}
ASSERT( mxIsDouble(mxVal) );
ASSERT( !mxIsComplex(mxVal) );
ASSERT( !mxIsSparse(mxVal) );
size_t n = mxGetNumberOfElements(mxVal);
vector<double> vec(n);
double* ptr = (double*)mxGetPr(mxVal);
for (size_t i=0;i<n;++i) {
vec[i] = ptr[i];
}
return vec;
}
//
// String
//
template<> string MexParams::getdef(const string& key, string def)
{
const mxArray* mxVal = params[key];
if (!mxVal) {
//throw runtime_error("No such key: " + key);
return def;
}
char buffer[1024];
if (mxGetString(mxVal, buffer, 1024)) {
throw runtime_error("Expected string parameter value");
}
return string(buffer);
}
#endif