Eigen3 Vector zu Standard-Array (C++) wandeln

[Freddy (Gast)]

Hallo,

ich versuche mich gerade in die Eigen3 C++ Bibliothek einzuarbeiten.
Zum Import von Standard C++ Arrays finde ich die Funktion MAP, aber wie 
kann ich einen Eigen3 Vector umcasten um ihn als Standard -Double-Array 
weiterzuverwenden?

Im folgenden Beispiel soll ein Testmuster berechnet werden und mit QWT 
geplottet werden:

so geht es

//#include "mainwindow.h"
//#include <QApplication>

#include <iostream>
#include <vector>

#include <qwt_plot_curve.h>
#include <qwt_plot.h>
#include <qapplication.h>

#include <Eigen/Dense>

using Eigen::MatrixXd;
using Eigen::VectorXd;
using Eigen::Vector3f;
using namespace std;

int main(int argc, char *argv[])
{
    QApplication a(argc, argv);

    double x[101];
    double y[101];

    for ( int i = 0; i < 101; i++ ) {
        x[i] =  i / 10.0;
        y[i] = sin(x[i]);
    }

    QwtPlot plot;
    QwtPlotCurve *curve = new QwtPlotCurve();
    curve->setRawSamples(x, y, 101);
    curve->attach(&plot);
    plot.show();

    return a.exec();
}*/

so soll es

int main(int argc, char *argv[])
{
    QApplication a(argc, argv);

    double x[101];
    double y[101];

    VectorXd rowX;
    VectorXd rowY;
    rowX.setLinSpaced(100, 0 , 10);
    rowY.setLinSpaced(100, 0 , 100);
    rowY = rowX.array().sin() ;

    std::cout << rowX << std::endl;
    std::cout << rowY << std::endl;


    //x = rowX.Map(1,100);
    //y = rowY.Map(1,100);

    QwtPlot plot;
    QwtPlotCurve *curve = new QwtPlotCurve();
    curve->setRawSamples(rowX, rowY, 101);
    curve->attach(&plot);
    plot.show();

    return a.exec();
}

Kennt jemand eine Lösung?

[Karl H. (kbuchegg)]

Ich fürchte mit einem simplen Umcasten wirst du da nicht weiter kommen.

Ein VectorXd ist laut Doku erst mal nichts anderes als eine Sonderform 
einer allgemeinen Matrix.
Für eine Matrix gibt es die Member-Funktion .data() welche dir einen 
Pointer auf die internen Daten liefert. Wie diese internen Daten 
allerdings organisiert sind, und ob dir das weiterhilft, kann ich nicht 
sagen. Im schlimmsten Fall wirst du eben eine Kopie der Daten machen 
müssen, in dem du mit dem Subscript-Operator auf die einzelnen Werte des 
Vektor in einer Schleife zugreift und im Prinzip ein Duplikat der Daten 
in deinen lokalen Arrays machst.
Ich denke aber, dass man mit .data weiterkommt. Zumnindest kligt es 
danach, als ob dort ein möglicher Ansatzpunkt zur Lösung des Problems 
liegt.

[Freddy (Gast)]

ok danke, dann bleibt mir wohl erstmal nichrts anderes übrig wie alles 
umzukopieren:

int main(int argc, char *argv[])
{
    QApplication a(argc, argv);

    double x[101];
    double y[101];

    VectorXd rowX;
    VectorXd rowY;
    rowX.setLinSpaced(100, 0 , 10);
    rowY.setLinSpaced(100, 0 , 100);
    //rowY = rowX.array().sin() ;
    //rowY = rowX.array().cos() ;
    rowY = rowX.array().tan();


    int i;
    for(i = 0; i != rowX.size(); i++) {
        //std::cout << row(i) << "  " << std::endl;
        x[i] = rowX(i);
        y[i] = rowY(i);
    }

    QwtPlot plot;
    QwtPlotCurve *curve = new QwtPlotCurve();
    curve->setRawSamples(x,y, rowX.size());
    curve->attach(&plot);
    plot.show();


    return a.exec();
}

[Freddy (Gast)]

und wenn man mag kann man den Plot noch ein wenig aufhübschen:

    QwtPlot plot;
    QwtPlotCurve *curve = new QwtPlotCurve();
    curve->setRawSamples(x,y, rowX.size());
    curve->setTitle("Sinus");
    curve->setPen( Qt::blue, 4 ),

    plot.setCanvasBackground(QColor(Qt::white));
    plot.setAxisTitle(QwtPlot::xBottom, "X-Axis");
    plot.setAxisTitle(QwtPlot::yLeft, "Y-Axis");
    plot.setAxisScale( QwtPlot::yLeft, -2, 2);


    curve->attach(&plot);
    plot.resize( 600, 400 );
    plot.show();

[Freddy (Gast)]

ein kleiner Screenshot