matlab - mex経由でC++からMatlab(vpa?)に多精度変数を渡す方法は?

Question

私はmexとMatlabが初めてです。実数変数と行列変数を mex 関数から Matlab に渡す方法を学習しました。

GMP (代替の mpir および mpfr ライブラリ) をインストールした後、私の質問は次のとおりです。C++ および GMP を介して mex 関数を構築している場合、Matlab のネイティブ関数のように「多精度」可変出力引数を渡す方法は?

簡単な実例はありますか？

以下は、複数の精度関連ヘッダーを追加した mex 関数コード サンプルです。「多倍精度」pi を 3 番目の出力引数として設定する方法は? またはprintf/mexPrintfを文字列として？

PS。正しくビルドするには、mpir、mpfr、および C++ ラッパー (mpreal.h) をインストールする必要があります。

#include <iostream>
#include <math.h>
#include "mex.h"

// The multiple precision libraries
// #include "mpir.h"
#include "mpreal.h"


using namespace mpfr;
using namespace std;

extern void _main();

/****************************/
class MyData {

public:
  void display();
  void set_data(double v1, double v2);
  MyData(double v1 = 0, double v2 = 0);
  ~MyData() { }
private:
  double val1, val2;
};

MyData::MyData(double v1, double v2)
{
  val1 = v1;
  val2 = v2;
}

void MyData::display()
{
#ifdef _WIN32
    mexPrintf("Value1 = %g\n", val1);
    mexPrintf("Value2 = %g\n\n", val2);
#else
  cout << "Value1 = " << val1 << "\n";
  cout << "Value2 = " << val2 << "\n\n";
#endif
}

void MyData::set_data(double v1, double v2) { val1 = v1; val2 = v2; }

/*********************/

static
void mexcpp(
        double num1,
        double num2
        )
{
#ifdef _WIN32
    mexPrintf("\nThe initialized data in object:\n");
#else
  cout << "\nThe initialized data in object:\n";
#endif
  MyData *d = new MyData; // Create a  MyData object
  d->display();           // It should be initialized to
                          // zeros
  d->set_data(num1,num2); // Set data members to incoming
                          // values
#ifdef _WIN32
  mexPrintf("After setting the object's data to your input:\n");
#else
  cout << "After setting the object's data to your input:\n";
#endif
  d->display();           // Make sure the set_data() worked
  delete(d);
  flush(cout);
  return;
}

void mexFunction(
         int          nlhs,
         mxArray      *[],
         int          nrhs,
         const mxArray *prhs[]
         )
{
  double      *vin1, *vin2;

  /* Check for proper number of arguments */

  if (nrhs != 2) {
    mexErrMsgIdAndTxt("MATLAB:mexcpp:nargin", 
            "MEXCPP requires two input arguments.");
  } else if (nlhs >= 1) {
    mexErrMsgIdAndTxt("MATLAB:mexcpp:nargout",
            "MEXCPP requires no output argument.");
  }

  vin1 = (double *) mxGetPr(prhs[0]);
  vin2 = (double *) mxGetPr(prhs[1]);


  ////////////// Here are the multiple precision settings
  const int digits = 64;
  mpreal::set_default_prec(mpfr::digits2bits(digits));
  const mpreal pi          =    mpfr::const_pi();
  ////////////// How can I set the multiple precision pi as a third output argument?
  ///// say vpa? or something else?
  mexPrintf(" The multiple precision pi is %g \n",pi);

  mexcpp(*vin1, *vin2);
  return;
}