SPAMS安装

SPAMS (SPArse Modeling Software)简介

是一个为解决各种稀疏估计问题的开源优化工具箱，其主页为http://spams-devel.gforge.inria.fr/index.html ，其可实现的功能如下：

• Dictionary learning and matrix factorization (NMF, sparse PCA, ...) 字典学习与矩阵分解 ，如（NMF非负矩阵分解，sparse PCA稀疏PCA，...）

• Solving sparse decomposition problems with LARS, coordinate descent, OMP, SOMP, proximal methods 使用LARS、快速下降法、OMP（Orthogonal Matching Pursuit，正交匹配追踪）、SOMP等算法解决稀疏分解问题

• Solving structured sparse decomposition problems (l1/l2, l1/linf, sparse group lasso, tree-structured regularization, structured sparsity with overlapping groups,...).

  解决结构化稀疏分解问题
  对于安装，下载安装包，解压到一个文件夹，如下图，解压后有不少文件、文件夹，其中的doc_spams.pdf是一个详尽的使用说明，介绍了各个函数接口如何使用（每个函数都有一个测试文件供测试，在test_release文件夹中）。
  [图片上传中...(image.png-c284ab-1511514252092-0)]

SPAMS安装

我的编译环境是window10, x64, vs2013, matlab2016a。安装主要分为两个步骤，设置编译器和编译compile.m文件
1.进入matlab，首先在 matlab命令窗口输入mex -setup，关联编译器，如果是Visual Studio的话，最好是2010版本以上（官网说明推荐）。

2.配置compile.m文件并编译

image.png

我自己安装的时候， 没有遇到找不到文件的问题， 在matlab把当前目录改为spams-matlab-**文件夹下就好了， compile.m用的是相对目录。

image.png image.png

如果报错的话， 就使用sublime或者notepad++批量修改文件路径为绝对路径
安装的时候，遇到isspace 不是 std 的函数问题， 打开prox/groups-graph.h文件， 把std::isspace 改成isspace就可以了。
当matlab命令行窗口没有红色报错的时候， 表示编译成功了， 截图

image.png

build文件下会新增很多文件

image.png

到此， 可以跑测试样例了
囧， 不能上传文件， 下面是我自己的compile.m文件

clear all;
get_architecture;

%%%%%%%%%%%%% COMPILER CONFIGURATION %%%%%%%%%%%%%%%%
% set up the compiler you want to use. Possible choices are
%    - 'mex' (default matlab compiler), this is the easy choice if your matlab
%              is correctly configured. Note that this choice might not compatible
%              with the option 'use_multithread=true'.
%    - 'icc' (intel compiler), usually produces the fastest code, but the
%              compiler is not free and not installed by default.
%    - 'gcc' (gnu compiler), good choice (for Mac, use gcc >= 4.6 for
%              the multi-threaded version, otherwise set use_multithread=false).
%              For windows, you need to have cygwin installed.
%    - 'clang'
%    - 'vs'  (visual studio compiler) for windows computers (10.0 or more is recommended)
%                for some unknown reason, the performance obtained with vs is poor compared to icc/gcc
compiler='mex';

 %%%%%%%%%%%% BLAS/LAPACK CONFIGURATION %%%%%%%%%%%%%%
% set up the blas/lapack library you want to use. Possible choices are
%    - builtin: blas/lapack shipped with Matlab,
%              same as mex: good choice if matlab is correctly configured.
%    - mkl: (intel math kernel library), usually the fastest, but not free.
%    - acml: (AMD Core math library), optimized for opteron cpus
%    - blas: (netlib version of blas/lapack), free
%    - atlas: (atlas version of blas/lapack), free,
% ==> you can also tweak this script to include your favorite blas/lapack library
blas='builtin';

%%%%%%%%%%%% MULTITHREADING CONFIGURATION %%%%%%%%%%%%%%
% set true if you want to use multi-threaded capabilities of the toolbox. You
% need an appropriate compiler for that (intel compiler, most recent gcc, or visual studio pro)
use_multithread=true; % (might not compatible with compiler=mex)
% if the compilation fails on Mac, try the single-threaded version.
% to run the toolbox on a cluster, it can be a good idea to deactivate this

use_64bits_integers=true;
% use this option if you have VERY large arrays/matrices
% this option allows such matrices, but may slightly reduce the speed of the computations.


use_mkl_threads=false;
% use this option is you use the mkl library and intends to use intensively BLAS3/lapack routines
% (for multiclass logistic regression, regularization with the trace norm for instance)
% this results in a loss of performance for many other functions

% if you use the options 'mex' and 'builtin', you can proceed with the compilation by
% typing 'compile' in the matlab shell. Otherwise, you need to set up a few path below.

path_matlab='/softs/bin/';
path_matlab='';

%%%%%%%%%%%% PATH CONFIGURATION %%%%%%%%%%%%%%%%%%%%
% only if you do not use the options 'mex' and 'builtin'
% set up the path to the compiler libraries that you intend to use below
add_flag='';
if strcmp(compiler,'gcc')
     if linux || mac
        % example when compiler='gcc' for Linux/Mac:    (path containing the files libgcc_s.*)
        path_to_compiler_libraries='/usr/lib/gcc/x86_64-redhat-linux/4.9.2/';
        path_to_compiler_libraries='/usr/lib/gcc/x86_64-redhat-linux/4.7.2/';
        path_to_compiler_libraries='/usr/lib/gcc/x86_64-linux-gnu/4.8.2/';
        path_to_compiler_libraries='/usr/lib/gcc/x86_64-linux-gnu/5/';
        path_to_compiler='/usr/bin/';
    else
        % example when compiler='gcc' for Windows+cygwin:    (the script does not
        % work at the moment in this configuration
        path_to_compiler='C:\cygwin\bin\';
        path_to_compiler_libraries='C:\cygwin\lib\gcc\i686-pc-cygwin\4.5.3\';
    end
elseif strcmp(compiler,'clang')
    path_to_compiler='/usr/bin/';
    path_to_compiler_libraries='/usr/lib/clang/3.5/lib/';
    path_to_libstd='/usr/lib/gcc/x86_64-linux-gnu/4.8/';
elseif strcmp(compiler,'open64')
    % example when compiler='gcc' for Linux/Mac:    (path containing libgcc_s.*)
    path_to_compiler_libraries='/opt/amdsdk/v1.0/x86_open64-4.2.4/lib/gcc-lib/x86_64-open64-linux/4.2.4/';
    path_to_compiler='/opt/amdsdk/v1.0/x86_open64-4.2.4/bin/';
elseif strcmp(compiler,'icc')
    if linux || mac
        % example when compiler='icc' for Linux/Mac
        path_to_gcccompiler_libraries='/usr/lib/gcc/x86_64-redhat-linux/4.7.2/';
        path_to_gcccompiler_libraries='/usr/lib/gcc/x86_64-linux-gnu/4.8/';
        path_to_gcccompiler_libraries='/usr/lib/gcc/x86_64-redhat-linux/4.9.2/';
        path_to_compiler_libraries='/opt/intel/composerxe/lib/intel64/';
        path_to_compiler='/opt/intel/composerxe/bin/';
        path_to_compiler='/scratch2/clear/mairal/intel/composerxe/bin/';
        path_to_compiler_libraries='/scratch2/clear/mairal/intel/composerxe/lib/intel64';
    else
        % example when compiler='icc' for Windows
        path_to_compiler_libraries='C:\Program Files (x86)\Intel\Composer XE\compiler\lib\intel64\';
        path_to_compiler='C:\Program Files (x86)\Intel\Composer XE\bin\intel64\';
        path_to_compiler_include='C:\Program Files (x86)\Intel\Composer XE\compiler\include\';
        path_to_vs='C:\Program Files (x86)\Microsoft Visual Studio 10.0\VC\bin\amd64\';
    end
elseif strcmp(compiler,'vs')
    % example when compiler='vs' for Windows
    path_to_compiler_libraries='C:\Program Files (x86)\Microsoft Visual Studio 10.0\VC\bin\amd64\';
    path_to_compiler=path_to_compiler_libraries;
elseif strcmp(compiler,'mex')
    % leave it blank when compiler='mex'
    path_to_compiler_libraries='';
    path_to_compiler='';
%    use_multithread=false;
end

% set up the path to the blas/lapack libraries.
if strcmp(blas,'mkl')
    if linux || mac
        path_to_blas='/opt/intel/composerxe/mkl/lib/intel64/';
        path_to_blas='/scratch2/clear/mairal/intel/composerxe/mkl/lib/intel64/';
    else
        path_to_blas='C:\Program Files (x86)\Intel\Composer XE\mkl\lib\intel64\';
    end
elseif strcmp(blas,'blas')
    if linux || mac
        path_to_blas='/usr/lib64/';
    else
        path_to_blas='?';
    end
elseif strcmp(blas,'atlas')
    if linux || mac
        path_to_blas='/usr/lib64/atlas/';
    else
        path_to_blas='?';
    end
elseif strcmp(blas,'acml')
    if linux || mac
        %path_to_blas=' /opt/amdsdk/v1.0/acml/open64_64/lib/';
        path_to_blas='/opt/acml5.0.0/ifort64/lib/';
    else
        path_to_blas='?';
    end
elseif strcmp(blas,'builtin')
    % leave it to /usr/lib/ for built-in:
    path_to_blas='/';
end

if mac
    add_flag=' -mmacosx-version-min=10.6';
end

debug=false;
if debug
    use_multithread=false;
end
%%%%%%%%%%%% END OF THE CONFIGURATION %%%%%%%%%%%%%%
% Do not touch what is below this line, unless you know what you are doing
out_dir='./build/';
mkdir(out_dir);

COMPILE = {
                '-I./linalg/ -I./prox/ prox/mex/mexSvmMisoOneVsRest.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexProjSplx.cpp',
                % compile dictLearn toolbox
                '-I./linalg/ -I./decomp/ -I./prox/ -I./dictLearn/ dictLearn/mex/mexArchetypalAnalysis.cpp',
                '-I./linalg/ -I./decomp/ -I./prox/ -I./dictLearn/ dictLearn/mex/mexTrainDL.cpp',
                '-I./linalg/ -I./decomp/ -I./prox/ -I./dictLearn/ dictLearn/mex/mexStructTrainDL.cpp',
                '-I./linalg/ -I./decomp/ -I./prox/ -I./dictLearn/ dictLearn/mex/mexTrainDL_Memory.cpp',
                % compile dag toolbox
                '-I./dags/ -I./linalg/ dags/mex/mexRemoveCyclesGraph.cpp',
                '-I./dags/ -I./linalg/ dags/mex/mexCountPathsDAG.cpp',
                '-I./dags/ -I./linalg/ dags/mex/mexCountConnexComponents.cpp',
                % compile proximal toolbox
                '-I./linalg/ -I./prox/ prox/mex/mexIncrementalProx.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexStochasticProx.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexFistaFlat.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexFistaTree.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexFistaGraph.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexFistaPathCoding.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexProximalFlat.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexProximalTree.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexProximalGraph.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexProximalPathCoding.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexEvalPathCoding.cpp',
                % compile linalg toolbox
                '-I./linalg/ linalg/mex/mexCalcAAt.cpp',
                '-I./linalg/ linalg/mex/mexCalcXAt.cpp',
                '-I./linalg/ linalg/mex/mexCalcXY.cpp',
                '-I./linalg/ linalg/mex/mexCalcXYt.cpp',
                '-I./linalg/ linalg/mex/mexCalcXtY.cpp',
                '-I./linalg/ linalg/mex/mexConjGrad.cpp',
                '-I./linalg/ linalg/mex/mexInvSym.cpp',
                '-I./linalg/ linalg/mex/mexSort.cpp',
                '-I./linalg/ linalg/mex/mexNormalize.cpp',
                % compile decomp toolbox
                '-I./linalg/ -I./dictLearn/ -I./decomp/ decomp/mex/mexDecompSimplex.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexOMP.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexLasso.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexLassoWeighted.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexRidgeRegression.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexCD.cpp'
                '-I./linalg/ -I./decomp/ decomp/mex/mexL1L2BCD.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexLassoMask.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexOMPMask.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexSOMP.cpp',
                '-I./linalg/ -I./decomp/ decomp/mex/mexSparseProject.cpp',
                % compile image toolbox
                '-I./image/ -I./linalg/ -I./prox/ image/mex/mexConvFista.cpp',
                '-I./image/ -I./linalg/ image/mex/mexExtractPatches.cpp',
                '-I./image/ -I./linalg/ image/mex/mexCombinePatches.cpp',
                % misc
                '-I./linalg/ linalg/mex/mexBayer.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexGraphOfGroupStruct.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexGroupStructOfString.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexReadGroupStruct.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexSimpleGroupTree.cpp',
                '-I./linalg/ -I./prox/ prox/mex/mexTreeOfGroupStruct.cpp'};

if linux || mac
    fid=fopen('run_matlab.sh','w+');
    fprintf(fid,'#!/bin/sh\n');
end

if sixtyfourbits
    if debug
        DEFCOMMON='-largeArrayDims -DDEBUG';
    else
        DEFCOMMON='-largeArrayDims -DNDEBUG';
    end
else
    if debug
        DEFCOMMON='-DDEBUG';
    else
        DEFCOMMON='-DNDEBUG';
    end
end
if windows
    DEFCOMMON=[DEFCOMMON ' -DWINDOWS -DREMOVE_'];
end


DEFBLAS='';
if strcmp(blas,'mkl')
    if use_mkl_threads
        name_mkl='mkl_intel_thread';
    else
        name_mkl='mkl_sequential';
    end
    DEFBLAS='-DUSE_BLAS_LIB -DAXPBY';
    if strcmp(arch,'GLNXA64')
        if use_64bits_integers
            blas_link = sprintf('-Wl,--start-group %slibmkl_intel_ilp64.a %slib%s.a %slibmkl_core.a -Wl,--end-group -ldl',path_to_blas,path_to_blas,name_mkl,path_to_blas);
        else
            blas_link = sprintf('-Wl,--start-group %slibmkl_intel_lp64.a %slib%s.a %slibmkl_core.a -Wl,--end-group -ldl',path_to_blas,path_to_blas,name_mkl,path_to_blas);
        end
    elseif strcmp(arch,'GLNX86')
        blas_link = sprintf('-Wl,--start-group %slibmkl_intel.a %slib%s.a %slibmkl_core.a -Wl,--end-group',path_to_blas,path_to_blas,name_mkl,path_to_blas);
    elseif strcmp(arch,'MACI64')
        if use_64bits_integers
            blas_link = sprintf('%slibmkl_intel_ilp64.a %slib%s.a %slibmkl_core.a',path_to_blas,path_to_blas,name_mkl,path_to_blas);
        else
            blas_link = sprintf('%slibmkl_intel_lp64.a %slib%s.a %slibmkl_core.a',path_to_blas,path_to_blas,name_mkl,path_to_blas);
        end
    elseif strcmp(arch,'MACI') || strcmp(arch,'MAC')
        blas_link = sprintf('%slibmkl_intel.a %slib%s.a %slibmkl_core.a',path_to_blas,path_to_blas,name_mkl,path_to_blas);
    elseif strcmp(arch,'PCWIN64')
        if use_64bits_integers
            blas_link = sprintf(' -L"%s" -lmkl_intel_ilp64 -l%s -lmkl_core',path_to_blas,name_mkl);
        else
            blas_link = sprintf(' -L"%s" -lmkl_intel_lp64 -l%s -lmkl_core',path_to_blas,name_mkl);
        end
    elseif strcmp(arch,'PCWIN')
        blas_link = sprintf(' -L"%s" -lmkl_intel -l%s -lmkl_core',path_to_blas,name_mkl);
    else
        'unsupported achitecture'
        return;
    end
elseif strcmp(blas,'blas')
    DEFBLAS='-DUSE_BLAS_LIB';
    blas_link='-lblas -l:liblapack.so.3';
elseif strcmp(blas,'atlas')
    DEFBLAS='-DUSE_BLAS_LIB -DADD_ATL -DAXPBY';
    blas_link='-l:libatlas.so.3 -l:liblapack.so.3';
elseif strcmp(blas,'builtin')
    blas_link='-lmwblas -lmwlapack';
    DEFBLAS='-DUSE_BLAS_LIB';
elseif strcmp(blas,'acml')
    DEFBLAS='-DUSE_BLAS_LIB -DNEW_MATLAB';
    blas_link=sprintf(' -L%s -L/usr/lib/ -lacml ',path_to_blas);
    %blas_link=sprintf(' -L%s -L/usr/lib/ -lacml -lg2c -L%s -lmv',path_to_blas,'/opt/amdsdk/v1.0/x86_open64-4.2.4/lib/gcc-lib/x86_64-open64-linux/4.2.4/');
else
    'please provide a correct blas library';
    return;
end


if strcmp(blas,'builtin')
    if ~verLessThan('matlab','7.9.0')
        DEFBLAS=[DEFBLAS ' -DNEW_MATLAB_BLAS'];
    else
        DEFBLAS=[DEFBLAS ' -DOLD_MATLAB_BLAS'];
    end
end
if use_64bits_integers
    DEFBLAS=[DEFBLAS ' -DINT_64BITS'];
end


links_lib=[blas_link];
link_flags=' -O ';

if strcmp(compiler,'icc')
    if windows
        DEFCOMP=['"PATH=%PATH%;' path_to_compiler '";"' path_to_vs '"  COMPILER=icl' ' -I"C:\Program Files (x86)\Intel\Composer XE\compiler\include"'];
        %compile_flags='/Qvc10 /fast /QaxSSE2,SSE3,SSE4.1,SSE4.2,AVX /Qansi-alias  /Qopenmp /fp:fast=2 /MD /Oy- /GR /EHs /Zp8';
        compile_flags='/Qvc10  /Qopenmp /MD /QaxSSE2,SSE3,SSE4.1,SSE4.2,AVX,CORE-AVX2,CORE-AVX-I /O2';
    else
        DEFCOMP=sprintf('CXX=%s/icpc LDCXX=%s/icpc',path_to_compiler,path_to_compiler);
        compile_flags='-fPIC -march=native -pipe -w -w0 -fast -fomit-frame-pointer -fno-alias -align -falign-functions';
        %link_flags=[link_flags sprintf(' -cxxlib=%s',path_to_gcc_libraries)];
        %link_flags=[link_flags ' -gcc-version=430'];
    end
    links_lib=[links_lib ' -L"' path_to_compiler_libraries '" '];
    %links_lib=[links_lib ' -L"' path_to_compiler_libraries '" -L"' path_to_gcc_libraries '"  -lstdc++'];
    if mac
        fprintf(fid,'export LIB_INTEL=%s\n',path_to_compiler_libraries);
        fprintf(fid,sprintf('export DYLD_LIBRARY_PATH=$DYLD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
        fprintf(fid,'export DYLD_INSERT_LIBRARIES=$LIB_INTEL/libimf.dylib:$LIB_INTEL/libintlc.dylib:$LIB_INTEL/libiomp5.dylib:$LIB_INTEL/libsvml.dylib\n');
    elseif linux
        fprintf(fid,'export LIB_GCC=%s\n',path_to_gcccompiler_libraries);
        fprintf(fid,'export LIB_INTEL=%s\n',path_to_compiler_libraries);
        fprintf(fid,sprintf('export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
        fprintf(fid,'export LD_PRELOAD=$LIB_INTEL/libimf.so:$LIB_INTEL/libintlc.so.5:$LIB_INTEL/libiomp5.so:$LIB_INTEL/libsvml.so:$LIB_GCC/libstdc++.so\n');
    end
    if use_multithread
        if windows
            compile_flags=[compile_flags ' /Qopenmp'];
            %links_lib=[links_lib ' /nodefaultlib:vcomp libiomp5md.lib'];
            links_lib=[links_lib sprintf('-L"%s" -liomp5md -lmmd -lirc -lsvml_dispmd -ldecimal',path_to_compiler_libraries)];
        else
            compile_flags=[compile_flags ' -openmp'];
            links_lib=[links_lib ' -liomp5'];
            if mac || linux
                fprintf(fid,'export KMP_DUPLICATE_LIB_OK=true\n');
            end
        end
    end
elseif strcmp(compiler,'open64')
    DEFCOMP=sprintf('CXX=%s/opencc',path_to_compiler);
%    fprintf(fid,'export LIB_INTEL=%s\n',path_to_compiler_libraries);
    fprintf(fid,sprintf('export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
    %fprintf(fid,'export LD_PRELOAD=$LIB_INTEL/libimf.so:$LIB_INTEL/libintlc.so.5:$LIB_INTEL/libiomp5.so:$LIB_INTEL/libsvml.so\n');
    compile_flags='-O2 -fomit-frame-pointer -march=opteron';
    links_lib=[links_lib ' -L"' path_to_compiler_libraries '" -L' path_to_blas ' -lstdc++'];
    if use_multithread
        compile_flags=[compile_flags ' -fopenmp'];
        links_lib=[links_lib ' -lopenmp'];
    end
elseif strcmp(compiler,'gcc')
    if windows
        DEFCOMP=['PATH=%PATH\%;' path_to_compiler ' COMPILER=g++-4' ' COMPFLAGS="-c -fexceptions" NAME_OBJECT=-o'];
    else
        DEFCOMP=sprintf('CXX=%s/g++',path_to_compiler);
    end
    if debug
        compile_flags='-O2 -g';
    else
        compile_flags='-O3 -mtune=native -fomit-frame-pointer -Wall';
    end
    links_lib=[links_lib ' -L"' path_to_compiler_libraries '" -L' path_to_blas];
    if mac
        fprintf(fid,'export LIB_GCC=%s\n',path_to_compiler_libraries);
        fprintf(fid,sprintf('export DYLD_LIBRARY_PATH=$DYLD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
        fprintf(fid,'export DYLD_INSERT_LIBRARIES=$LIB_GCC/libgfortran.so:$LIB_GCC/libgcc_s.so:$LIB_GCC/libstdc++.so:$LIB_GCC/libgomp.so\n');
    elseif linux
        fprintf(fid,'export LIB_GCC=%s\n',path_to_compiler_libraries);
        fprintf(fid,sprintf('export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
        fprintf(fid,'export LD_PRELOAD=$LIB_GCC/libgfortran.so:$LIB_GCC/libgcc_s.so:$LIB_GCC/libstdc++.so:$LIB_GCC/libgomp.so\n');
    end
    if use_multithread
        compile_flags=[compile_flags ' -fopenmp'];
        links_lib=[links_lib ' -lgomp'];
    end
elseif strcmp(compiler,'clang') ||strcmp(compiler,'clang++-libc++')
    DEFCOMP=sprintf('CXX=%s/clang++ -Dchar16_t=uint16_T',path_to_compiler);
    compile_flags='-stdlib=libstdc++ -O3 -mtune=native -fomit-frame-pointer -Wall';
    links_lib=[links_lib ' -L"' path_to_compiler_libraries '" -L' path_to_blas];
    if mac
        fprintf(fid,'export LIB_GCC=%s\n',path_to_compiler_libraries);
        fprintf(fid,sprintf('export DYLD_LIBRARY_PATH=$DYLD_LIBRARY_PATH:%s:%s\n',path_to_compiler_libraries,path_to_blas));
    elseif linux
        fprintf(fid,'export LIB_GCC=%s\n',path_to_compiler_libraries);
        fprintf(fid,'export LIB_STD=%s\n',path_to_libstd);
        fprintf(fid,sprintf('export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:%s:%s:\n',path_to_compiler_libraries,path_to_blas,path_to_libstd));
        fprintf(fid,'export LD_PRELOAD=$LIB_STD/libstdc++.so\n');
    end
    if use_multithread
        compile_flags=[compile_flags ' -fopenmp'];
        links_lib=[links_lib ' -lgomp'];
    end

elseif strcmp(compiler,'vs')
    DEFCOMP='COMPILER=cl';
    compile_flags='/c /02';
    if use_multithread
        compile_flags=[compile_flags ' /openmp'];
        links_lib=[links_lib ' -lvcomp'];
    end
elseif strcmp(compiler,'mex')
    DEFCOMP='';
    compile_flags=' -O';
    if mac || linux
        fprintf(fid,sprintf('export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:%s\n',path_to_blas));
    end
    if use_multithread
        if (linux || mac)
            compile_flags=[compile_flags ' -fopenmp']; % we assume gcc
            links_lib=[links_lib ' -lgomp'];
        end
    end
else
    'unknown compiler'
    return;
end

if ~windows
    if use_mkl_threads
        fprintf(fid,[path_matlab 'matlab $* -singleCompThread -r \"addpath(''./build/''); addpath(''./test_release''); setenv(''MKL_NUM_THREADS'',''4''); setenv(''MKL_SERIAL'',''NO'');setenv(''MKL_DYNAMIC'',''FALSE'');"\n']);
    else
        fprintf(fid,[path_matlab 'matlab $* -singleCompThread -r \"addpath(''./build/''); addpath(''./test_release''); setenv(''MKL_NUM_THREADS'',''1''); setenv(''MKL_SERIAL'',''YES'');"\n']);
    end
    fclose(fid);
    !chmod +x run_matlab.sh
end

DEFS=[DEFBLAS ' ' DEFCOMMON ' ' DEFCOMP];

if mac
    compile_flags=[compile_flags add_flag];
end

for k = 1:length(COMPILE),
    str = COMPILE{k};
    fprintf('compilation of: %s\n',str);
    if windows
        str = [str ' -outdir ' out_dir, ' ' DEFS ' ' links_lib ' OPTIMFLAGS="' compile_flags '" '];
    else
        if verLessThan('matlab','8.3.0')
            str = [str ' -outdir ' out_dir, ' ' DEFS ' CXXOPTIMFLAGS="' compile_flags '" LDOPTIMFLAGS="' link_flags '" ' links_lib];
        else
            str = [str ' -outdir ' out_dir, ' ' DEFS ' CXXOPTIMFLAGS="' compile_flags '" LDOPTIMFLAGS="' link_flags '" ' ' LINKLIBS="$LINKLIBS ' links_lib '" '];
        end
    end
    args = regexp(str, '\s+', 'split');
    args = args(find(~cellfun(@isempty, args)));
    mex(args{:});
end

copyfile src_release/*.m build/