In an earlier post I did a review of different C/C++ libraries for scientific computing (among them IT++, GSL, and others) and concluded that the armadillo linear algebra library has a lot of functions, is very easy to use, and very fast. So, I thought it would be very nice to make it work within mex functions. This could make functions in C++ both very legible and easy to write at the same time fast.

I managed to do this and in this post I explain how.

It is kind of a hack, but works smoothly without problems, although the compilation was first a hassle that took some time. See the end of this post on how to compile.

The trick is to first create an armadillo structure and then let its pointer indicate the input structure from matlab (dimensionality of the structure has to be adapted, if not you get segfaulted). When the function ends, armadillo wants to destroy its structures, so before, you have to make it point to the original location and adapt dimensionality (or else...).

Here comes a toy example of a mex file that should illustrate how to use it with matlab.

## Compilation:

If you don't have it yet, take 5 minutes to install the armadillo library.For compilation you should include the additional compiler switches -lcblas -DARMA_NO_DEBUG. In an ideal world you would only have to type this:

>> mex toy_example.cpp -lcblas -DARMA_NO_DEBUG

Unfortunately this is not an ideal world. While this compiled for me without problems, on execution of the function within matlab, I got this error a lot:

invalid mex file.

GLIBCXX_3.4.9' not found

I first sym-linked matlab-provided shared libraries (.so files) to the ones of my system. I am not sure this is a good idea and anyways it didn't solve the problem. I edited heavily the mexopts.sh file which gives the parameters to mex compiler, however this did not help, it didn't seem to take these parameters. I finally compiled from bash prompt by a variation of the commands I found mex used (use -v switch with mex). The commands that work for me are these (you might have to adapt them slightly):

>> g++ -c -I/opt/matlab/extern/include -I/opt/matlab/simulink/include -DMATLAB_MEX_FILE -ansi -fPIC -fno-omit-frame-pointer -pthread -DMX_COMPAT_32 -DNDEBUG -DARMA_NO_DEBUG "toy_example.cpp"

>>g++ -c -I/opt/matlab/extern/include -I/opt/matlab/simulink/include -DMATLAB_MEX_FILE -ansi -D_GNU_SOURCE -fexceptions -fPIC -fno-omit-frame-pointer -pthread -DMX_COMPAT_32 -O3 -DNDEBUG "/opt/matlab/extern/src/mexversion.c"

>>g++ -O3 -pthread -shared -Wl,--version-script,/opt/matlab/extern/lib/glnxa64/mexFunction.map -Wl,--no-undefined -o "toy_example.mexa64" mahal_dist.o mexversion.o -Wl,-rpath-link,/opt/matlab/bin/glnxa64 -L/opt/matlab/bin/glnxa64 -lmx -lmex -lmat -lm -lcblas

Hope you enjoy using this. Please also see my mex implementation of the mahalanobis function.

Happy coding! Please leave comments below.

New bench IT++ vs Armadillo

ReplyDeletehttps://sourceforge.net/projects/itpp/forums/forum/115656/topic/3740250

Your comments are welcome!

@Anonymous : Thanks for the link! There's an interesting discussion on speed differences between armadillo and IT++.

ReplyDeleteI must say I really like it. Your imformation is usefull. Thanks for share.

ReplyDeleteNice article! I am very pleased to be visiting

ReplyDeleteyour site and read the information that you post. From your site, I can learn

more about the armadillo. Initially I found it difficult to do it but because

your site now I can understand about armadillos and do not feel any more

trouble.

that's very interesting application! I would like to introduce it in my research in the next step! thanks for sharing

ReplyDeleteI had difficulty compiling your example. What library does access:rw come from? Thanks in advance!

ReplyDeleteGreat article !

ReplyDeleteI had no problem compiling with MATLAB R2012a. Please see this link to configure your MATLAB installation in order to use the system libraries in place of matlab ones :

http://judsonsnotes.com/notes/index.php?option=com_content&view=article&id=659%3ainstalling-matlab-in-ubuntu-1110&catid=37%3atech-notes&Itemid=59

Simpson told her she would like Coach Outlet Online the ending. Over the years Lisa had an on-and-off relationship with Simpson, but it would be closer in some ways than the one she had with her father.

ReplyDelete

ReplyDeletethanks so much for that great blog and thanks also for accepting my links thanks

طريقة عمل الدونات طريقة عمل البان كيك طريقة عمل الكنافة طريقة عمل البسبوسة طريقة عمل الكيك طريقة عمل عجينة البيتزا فوائد القرفه

thanks so much i like very so much your post

فوائد الحلبة فوائد الزنجبيل فوائد الرمان فوائد زيت السمسم علاج البواسير فوائد البصل فوائد اليانسون فوائد الكركم فوائد الزعتر علامات الحمل

thanks a lot for that great article (y) (y)

ReplyDeleteGateau au chocolat

Recette Gateau au Yaourt

recette poulet basquaise

recette cake salé

soupe au choux

Soupe a l oignon

cake au jambon

recette cookies