Я пытаюсь построить два вектора, где я сохранил элементы двух mxArray
s из MATLAB (с использованием компилятора Visual C ++).
test1
а также test2
являются mxArray
s из MAT-файла, и когда я их печатаю, все в порядке. Но когда я запускаю программу, фигура пуста.
#include "mat.h"#include "NoneDynamic.h"#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <engine.h>
using namespace std;
mxArray *T = NULL, *D = NULL;
Engine *ep;
void matread(const char *file, std::vector<double>& v)
{
// open MAT-file
MATFile *pmat = matOpen("DataTesting.mat", "r");
if (pmat == NULL) return;
// extract the specified variable
mxArray *arr = matGetVariable(pmat, file);
if (arr != NULL && mxIsDouble(arr) && !mxIsEmpty(arr)) {
// copy data
mwSize num = mxGetNumberOfElements(arr);
double *pr = mxGetPr(arr);
if (pr != NULL) {
v.resize(num);
v.assign(pr, pr + num);
}
}
// cleanup
mxDestroyArray(arr);
matClose(pmat);
}
int main()
{
SeenData object;
vector<double> v;
matread("Test1", v);
vector <double> v1;
matread("test2", v1);
if (!(ep = engOpen(""))) {
fprintf(stderr, "\nCan't start MATLAB engine\n");
return EXIT_FAILURE;
}
T = mxCreateDoubleMatrix(1,100, mxREAL);
memcpy((vector<double> *)mxGetPr(T), (vector<double>*) &v, sizeof(v));
D = mxCreateDoubleMatrix(1,100, mxREAL);
memcpy((vector<double> * )mxGetPr(D), (vector<double>*) &v1, sizeof(v1));
engPutVariable(ep, "T", T);
engPutVariable(ep, "D", D);
engEvalString(ep, "plot (T,D)");
fgetc(stdin);// for pausing long enough to see the plot
engEvalString(ep, "close;");
mxDestroyArray(T);
mxDestroyArray(D);
engClose(ep);
}
Со следующими инструкциями
memcpy((vector<double> *)mxGetPr(T), (vector<double>*) &v, sizeof(v));
memcpy((vector<double> * )mxGetPr(D), (vector<double>*) &v1, sizeof(v1));
Я полагаю, ваше намерение было скопировано в mxGetPr(T)
а также mxGetPr(D)
данные, содержащиеся в v
а также v1
,
Но
1) std::vector<double>
класс с (я полагаю) указателем в нем; так sizeof(v)
а также sizeof(v1)
не размер содержащихся данных, а фиксированный размер (24 для моего clang ++)
2) мем копирует указатель v
(а также v1
) копирует структурный вектор, а не данные, содержащиеся в нем.
3) std::memcpy()
получать пустые указатели; так что ваш актерский состав бесполезен
если ты действительно хочу использовать std::memcpy()
и вы используете компилятор C ++ 11 / C ++ 14, я полагаю, вы можете написать
memcpy(mxGetPr(T), v.data(), sizeof(double)*v.size());
memcpy(mxGetPr(D), v1.data(), sizeof(double)*v1.size());
иначе старый добрый
double * tp = mxGetPr(T);
for ( std::size_t ui = 0 ; ui < v.size() ; ++ui )
tp[ui] = v[ui];double * td = mxGetPr(D);
for ( std::size_t ui = 0 ; ui < v1.size() ; ++ui )
td[ui] = v1[ui];
p.s .: извините за мой плохой английский
Других решений пока нет …