Difference between revisions of "ApCoCoA-1:MatlabInterface"

From ApCoCoAWiki
Line 133: Line 133:
  
 
<matlab>Pts = [-1,0;0,0;1,0;99,1;99,0;99,-1];
 
<matlab>Pts = [-1,0;0,0;1,0;99,1;99,0;99,-1];
PreProcessing([0.9,0.9], Pts, 2)
+
Preprocess([0.9,0.9], Pts, 2)
  
 
ans =
 
ans =

Revision as of 12:28, 25 February 2008

MatlabToolbox Interface

The ApCoCoALib is written in C++ and therefore the potential number of users is limited to people who are comfortable writing C++ code. But many users of Computer Algebra systems or other software might benefit from the functionality provided by the CoCoALib. One such piece of software is Matlab which among other things allows the user to integrate C, Fortran and C++ libraries. The integration is not too difficult for libraries written in C or Fortran but there are severe restrictions for libraries written in C++. Due to that fact we are currently working on the necessary bits to provide the functionality of the CoCoALib in MatLab.

We want to provide the MatLabToolbox as binary blob as well as in source. Since the CoCoALib runs on a wide number of 32 and 64 bit platforms we expect to provide binaries for

  • x86 & x86-64 Linux
  • x86 & x86-64 Windows
  • 32bit PPC MacOSX
  • Sparc Solaris

This obviously depends on access to appropriate hardware and MatLab binaries.

Downloads

Caution: The current source codes are prepared to compile with MSVC 2005 only. The use of any other compiler but MSVC 2005 will lead to spectacular link failures since one cannot link C++ binaries produced by MSVC 2005 with for example MSVC 2003. You also need to have the runtime of MSVC 2005 installed on your system to make this work. In case you experience problems you can always ask questions in the ApCoCoA forum.

The process of building the MatlabInterface is described at Compile on Windows. Links for downloads are also offered there.

Design

C++ with exceptions and RTTI vs. MatLab

The design is quite complex due to the fact that integrating C++ code into a mex file has severe restrictions, most significantly:

  • Limitation to g++ 3.2.2 on Linux. Using g++ 3.2.1 or 3.2.3 might lead to crashes or unexpected behaviour in general. While this might have been ok two years ago most current systems now ship with g++ 4.0 or 4.1.
  • The Microsoft compiler does disable exception handling or RTTI per default. Mex also explicitly disables exceptions in C++ code in the mexopts.bat file. While one obviously is free to change that we have been unable to create a mex plugin compiled with exceptions and RTTI. Since those two features are essential to the CoCoALib a different way had to be found.

It has also been reported that migrating C++ code from R13 to R14 has usually not been as easy as one could hope. This is obviously only anecdotal evidence, as each of the google groups on the other hands provides plenty of people seeking help.

Using C wrapper around the ApCoCoALib

The solution is to use extern "C" functions inside a second DLL. That second DLL can be compiled with exceptions and RTTI. For a demonstration look at the following two code snippets demonstrating the integration of the PreProcessing function:

ApCoCoALib.cpp: <cpp>#include "CoCoA/library.H"

  1. include "ApCoCoA/Library.H"
  2. include "mex.h"
  3. include <vector>

using std::vector;

extern "C" __declspec(dllexport) int LibPreprocess(mxArray *Tolerance, mxArray *OrigArray,

                                               int aType, mxArray *PreProcessedArray){
 
 double *tolerance = mxGetPr(Tolerance);
 double *OrigData = mxGetPr(OrigArray);
 int rows = (int) mxGetM(OrigArray);
 int cols = (int) mxGetN(OrigArray);
 CoCoA::GlobalManager CoCoAFoundations;
 
 std::vector<double> toler(cols);
 std::vector<CoCoA::ApproxPts::ApproxPt> OrigPts(rows, CoCoA::ApproxPts::ApproxPt(cols));
 std::vector<CoCoA::ApproxPts::ApproxPt> PreprocessedPts;
 // Build tolerance parameter
 for (int i=0; i < cols; i++)
   {toler[i]=tolerance[i];}
 // Copy original data to vector
 for (int i=0; i < rows; ++i) {
   for (int j=0; j < cols; ++j){
     OrigPts[i][j]=OrigData[j*rows+i];
   }
 }
 // Call Preprocess function in ApCoCoALib
 switch ( aType ) {
   case 1 : 

CoCoA::ApproxPts::PreprocessGridAlgm(PreprocessedPts, OrigPts, toler); break;

   case 2 : 

CoCoA::ApproxPts::PreprocessAggrAlgm(PreprocessedPts, OrigPts, toler); break;

   case 3:

CoCoA::ApproxPts::PreprocessSubdivAlgm(PreprocessedPts, OrigPts, toler); break;

 }
 // Rearrange preprocessed points
 int rowsRes = PreprocessedPts.size();
 double *content;
 size_t noBytes = cols * rowsRes * sizeof(double);
 if ((content = (double*)mxRealloc(mxGetPr(PreProcessedArray), noBytes))==NULL){
   return -1;
 }
 mxSetPr(PreProcessedArray, content);
 mxSetM(PreProcessedArray, rowsRes);
 mxSetN(PreProcessedArray, cols);
 
 for (int i=0; i < rowsRes; i++) {
   for (int j=0; j < cols; j++) {

content[j*rowsRes+i] = PreprocessedPts[i][j];

   }
 }
 return 0;

}</cpp>

and Preprocess.c <c> #include "mex.h"

 #include "matrix.h"
 __declspec(dllimport) int LibPreprocess(mxArray *Tolerance, mxArray *OrigArray, int aType, mxArray *PreProcessedArray);
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])

{

 mxArray *myRes;
 int algo, retCode;
 if (nrhs != 3) 
   {mexErrMsgTxt("3 input arguments required.");}
 algo = mxGetScalar(prhs[2]);
 myRes = mxCreateDoubleMatrix(0, 0, mxREAL);
 // Run calculation
 retCode = LibPreprocess(prhs[0], prhs[1], algo, myRes);
 if (retCode==-1){
   mexErrMsgTxt("Error in Preprocess.c: Reallocation failed. Out of memory.");
 }
 // Assign to output
 plhs[0] = myRes;

}</c>

We use mex.bat to compile the second file into a mex plugin (as C code) while we use cl.exe to compile the first file into a DLL.

Functions provided

Obviously the goal is to provide every available function in the ApCoCoALib to users of Matlab. Due to our own needs at the moment the first functions to work will be:

  • PreProcessing
  • Approximate Buchberger-Möller
  • Syzygy computation for ideals

PreProcessing

The third example given in the CoCoA documentation (see PreprocessPts5) can be derived in Matlab with the following commands:

<matlab>Pts = [-1,0;0,0;1,0;99,1;99,0;99,-1]; Preprocess([0.9,0.9], Pts, 2)

ans =

    0     0
   99     0

</matlab>