// callMatlab.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <windows.h>
#include <matrix.h>
#include "../liba.h"

class mxMatDbl
{
    
public:
    mxArray *m_array;
    mxMatDbl()
    {
        m_array=0;
    }
    mxMatDbl(int h,int w,const double * pData)
    {
        m_array=mxCreateDoubleMatrix(h,w, mxREAL); 
        if(pData)
            memcpy(mxGetPr(m_array),pData,h*w*sizeof(double));
    }

    double * getPtr()
    {
        return mxGetPr(m_array);
    }
    operator double *()
    {
        return mxGetPr(m_array);
    }
    operator mxArray *()
    {
        return m_array;
    }
    double & operator ()(int y,int x)
    {
        assert(m_array);
        int W=w();
        int H=h();
        assert(y>=0 && y<H && x>=0 && x<W);
        return *(mxGetPr(m_array)+ y*W + x);
    }

    int h()
    {
        return mxGetM(m_array);
    }
    int w()
    {
        return mxGetN(m_array);
    }
    void Release()
    {
        if(m_array)
        {
            mxDestroyArray(m_array);
            m_array=0;
        }
    }

    ~mxMatDbl()
    {
        Release();
    }

    void Print()
    {
        for (int y=0;y<h();y++)
        {
            for (int x=0;x<w();x++)
            {
                printf("%g ",(*this)(y,x));
            }
            printf("\n");
        }
    }
protected:
private:
};

void TestMxMatDbl()
{
    double data[6]={1.0 ,2.0 ,3.0 ,4.0,5.0,6.0};
    mxMatDbl m(2,3,data);
    printf("w %d,h %d\n",m.w(),m.h());
    for (int y=0;y<m.h();y++)
    {
        for (int x=0;x<m.w();x++)
        {
            printf("%g ",m(y,x));
        }
        printf("\n");
    }
}

void callMatlab()
{
    libaInitialize();
    double _r[]={0,0,0};
    mxMatDbl r(3,1,_r);
    mxMatDbl R(3,3,0);
    if(0)
    {//这一局怎么不行？？
        mxArray *prhs[1];prhs[0]=r.m_array;
        mxArray *plhs[1];plhs[0]=R.m_array;
        mlxOula(1, plhs, 1, prhs);
    }
    else
        mlxOula(1, &R.m_array, 1, &r.m_array);
    r.Print();
    R.Print();
    libaTerminate();
}
int main(int argc, char* argv[])
{
    //TestMxMatDbl();
    callMatlab();
    return 0;
}