#! /bin/bash
#
# build up atlas_xxx.dll for use with MATLAB.
#
# 1) build atlas libraries (make config , make install)
# 2) run this script in the Atlas/lib/<arch> directory containing the created *.a files
#
# notes:
#  1) assumes cygwin install including mingw stuff
#  2) assumes you have the mingw libg2c.a library in /lib/mingw that is NOT PART OF CYGWIN distro
#     this is in the mingw distribution (www.mingw.org)

# base name of the DLL, .def, and import libraries (.lib) to be created
AtlasName=atlas-3.6.0-WinNT-P4SSE2
# set names for specific output files (see notes on MSVC lib.exe below)
defname=${AtlasName}.def
dllname=${AtlasName}.dll
gcclibname=${AtlasName}-gcc.lib
msvclibname=${AtlasName}-msvc.lib

# create the .def file
#minimal set of Fortran exports needed to produce an Atlas callable by Matlab 6.x
#MatlabExports=(\
#	dasum daxpy dcopy ddot dgbmv dgemm dgemv dger dgesv dgetrf \
#     dgetrs dnrm2 dposv dpotrf dpotrs drot drotg dsbmv dscal dspmv \
#	dspr dspr2 dswap dsymm dsymv dsyr dsyr2 dsyr2k dsyrk dtbmv \
#	dtbsv dtpmv dtpsv dtrmm dtrmv dtrsm dtrsv dzasum dznrm2 idamax \
#	izamax xerbla zaxpy zcopy zdotc zdotu zdscal zgbmv zgemm zgemv \
#	zgerc zgeru zgesv zgetrf zgetrs zhbmv zhemm zhemv zher zher2 \
#	zher2k zherk zhpmv zhpr zhpr2 zposv zpotrf zpotrs zrotg zscal \
#	zswap zsymm zsyr2k zsyrk ztbmv ztbsv ztpmv ztpsv ztrmm ztrmv \
#	ztrsm ztrsv)

# C BLAS interface
Cexports=(\
#single precision routines
#	cblas_sdsdot cblas_sdot cblas_snrm2 cblas_sasum cblas_scnrm2 \
#	cblas_scasum cblas_sscal cblas_isamax \
#	cblas_sswap cblas_scopy cblas_saxpy catlas_saxpby catlas_sset \
#	cblas_srotg cblas_srotmg cblas_srot cblas_srotm \
#	cblas_sgemv cblas_sgbmv cblas_strmv cblas_stbmv cblas_stpmv \
#	cblas_strsv cblas_stbsv cblas_stpsv \
#	cblas_ssymv	cblas_ssbmv cblas_sspmv cblas_sger cblas_ssyr \
#	cblas_sspr cblas_ssyr2 cblas_sspr2 \
#	cblas_sgemm cblas_ssymm cblas_ssyrk cblas_ssyr2k cblas_strmm cblas_strsm \
#double precision routines
#	cblas_dsdot cblas_ddot cblas_dnrm2 cblas_dasum cblas_dznrm2 \
#	cblas_dzasum cblas_dscal cblas_idamax \
#	cblas_dswap cblas_dcopy cblas_daxpy catlas_daxpby catlas_dset \
#	cblas_drotg cblas_drotmg cblas_drot cblas_drotm \
#	cblas_dgemv cblas_dgbmv cblas_dtrmv cblas_dtbmv cblas_dtpmv \
#	cblas_dtrsv cblas_dtbsv cblas_dtpsv \
#	cblas_dsymv cblas_dsbmv cblas_dspmv cblas_dger cblas_dsyr \
#	cblas_dspr cblas_dsyr2 cblas_dspr2 \
#	cblas_dgemm cblas_dsymm cblas_dsyrk cblas_dsyr2k cblas_dtrmm cblas_dtrsm \
#single precision complex routines
#	cblas_cdotu_sub cblas_cdotc_sub cblas_cscal cblas_csscal cblas_icamax \
#	cblas_cswap cblas_ccopy cblas_caxpy catlas_caxpby catlas_cset \
#	cblas_crotg cblas_csrot \
#	cblas_cgemv cblas_cgbmv	cblas_ctrmv cblas_ctbmv cblas_ctpmv \
#	cblas_ctrsv cblas_ctbsv cblas_ctpsv \
#	cblas_chemv cblas_chbmv cblas_chpmv cblas_cgeru cblas_cgerc \
#	cblas_cher cblas_chpr cblas_cher2 cblas_chpr2 \
#	cblas_cgemm cblas_csymm cblas_csyrk cblas_csyr2k cblas_ctrmm cblas_ctrsm \
        cblas_zgemm \
#	cblas_chemm cblas_cherk cblas_cher2k \
#double precision complex routines
#	cblas_zdotu_sub cblas_zdotc_sub cblas_zscal cblas_zdscal cblas_izamax \
#	cblas_zswap cblas_zcopy cblas_zaxpy catlas_zaxpby catlas_zset \
#	cblas_zrotg cblas_zdrot \
#	cblas_zgemv cblas_zgbmv cblas_ztrmv cblas_ztbmv cblas_ztpmv	\
#	cblas_ztrsv cblas_ztbsv cblas_ztpsv \
#	cblas_zhemv cblas_zhbmv cblas_zhpmv cblas_zgeru cblas_zgerc \
#	cblas_zher cblas_zhpr cblas_zher2 cblas_zhpr2 \
#	cblas_zgemm cblas_zsymm cblas_zsyrk cblas_zsyr2k cblas_ztrmm cblas_ztrsm \
        cblas_zgemm \
#	cblas_zhemm cblas_zherk cblas_zher2k \
#other stuff
	cblas_errprn \
# CLAPACK interface subset
#	clapack_sgesv clapack_sgetrf clapack_sgetrs clapack_sgetri \
#	clapack_sposv clapack_spotrf clapack_spotrs clapack_spotri \
#	clapack_slauum clapack_strtri \
	clapack_dgesv clapack_dgetrf clapack_dgetrs clapack_dgetri \
	clapack_dposv clapack_dpotrf clapack_dpotrs clapack_dpotri \
	clapack_dlauum clapack_dtrtri \
#	clapack_cgesv clapack_cgetrf clapack_cgetrs clapack_cgetri \
#	clapack_cposv clapack_cpotrf clapack_cpotrs clapack_cpotri \
#	clapack_clauum clapack_ctrtri \
	clapack_zgesv clapack_zgetrf clapack_zgetrs clapack_zgetri \
	clapack_zposv clapack_zpotrf clapack_zpotrs clapack_zpotri \
	clapack_zlauum clapack_ztrtri)
#Fortran BLAS interface
FortranExports=(\
#single precision routines
#	saxpy scasum scnrm2 scopy sdot sdsdot sgbmv sgemm \
#	sgemv sger snrm2 srot srotg srotm srotmg ssbmv sscal \
#	sspmv sspr sspr2 sswap ssymm ssymv ssyr ssyr2 ssyr2k \
#	ssyrk stbmv stbsv stpmv stpsv strmm strmv strsm strsv \
#double precision routines (dgemv and dtrsv are needed for complex SuperLU)
#	dasum daxpy dcabs1 dcopy ddot dgbmv dgemm \
#	dger dnrm2 drot drotg drotmg dsbmv dscal dsdot \
#	dspmv dspr dspr2 dswap dsymm dsymv dsyr dsyr2 dsyr2k \
#	dsyrk dtbmv dtbsv dtpmv dtpsv dtrmm dtrmv dtrsm \
#	dzasum dznrm2 icamax idamax isamax izamax lsame sasum \
	dgemv dtrsv \
#single precision complex routines
#	caxpy ccopy cdotc cdotu cgbmv cgemm cgerc cgeru chbmv \
#	chemm chemv cher cher2 cher2k cherk chpmv chpr chpr2 \
#	crotg cscal csrot csscal cswap csymm csyr2k csyrk \
#	ctbmv ctbsv ctpmv ctpsv ctrmm ctrmv ctrsm ctrsv \
#double precision complex routines (needed for complex SuperLU)
	zaxpy zcopy zdotc zdotu zdrot zdscal zgbmv zgemm \
	zgemv zgerc zgeru zhbmv zhemm zhemv zher zher2 zher2k \
	zherk zhpmv zhpr zhpr2 zrotg zscal zswap zsymm \
	zsyr2k zsyrk ztbmv ztbsv ztpmv ztpsv ztrmm ztrmv ztrsm ztrsv \
#other stuff
	xerbla \
# LAPACK interface subset
#	sgesv sgetrf sgetrs sgetri sposv spotrf spotrs spotri \
#	slauum strtri \
#	dgesv dgetrf dgetrs dgetri dposv dpotrf dpotrs dpotri \
#	dlauum dtrtri \
#	cgesv cgetrf cgetrs cgetri cposv cpotrf cpotrs cpotri \
#	clauum ctrtri \
	zgesv zgetrf zgetrs zgetri zposv zpotrf zpotrs zpotri \
	zlauum ztrtri)

#create exports file
echo "EXPORTS" > ${defname}
#write C exports (these are straight C, so do not need any name mangling or underscores)
for (( i = 0 ; i < ${#Cexports[*]} ; i++ )) ; do
	export=${Cexports[$i]}
	echo -e "${export}=${export}" >> ${defname}
done
#write Fortran exports (have to append trailing underscore)
for (( i = 0 ; i < ${#FortranExports[*]} ; i++ )) ; do
	export=${FortranExports[$i]}
	echo -e "${export}=${export}_" >> ${defname}
done

# libraries
CLIBPATH=/usr/lib/mingw
clib="$CLIBPATH/libg2c.a $CLIBPATH/libmoldname.a $CLIBPATH/libmsvcrt.a"

# Link Atlas DLL and create gcc import library
echo 'Linking DLL and creating gcc import library...'
gcc -mno-cygwin -shared -o ${dllname} ${defname} \
    liblapack.a libcblas.a libf77blas.a libatlas.a \
    -Wl,--out-implib=${gcclibname} \
    -Wl,--enable-auto-import \
    -Wl,--no-whole-archive ${clib}
#    -Wl,--export-all-symbols \
#-Wl,--whole-archive 

# Create import library for use with MS Visual C++
# VC's lib.exe tool produces a .lib which references a .dll with the same name as
# the .def that it's given, so those two filenames have to match
echo 'Creating MSVC import library...'
# paths for Visual Studio 6.0
#LIBPATH=`cygpath -u "C:\Program Files\Microsoft Visual Studio\VC98\bin"`
#VBPATH=`cygpath -u "C:\Program Files\Microsoft Visual Studio\Common\MSDev98\Bin"`
# paths for Visual Studio .NET
LIBPATH=`cygpath -u "C:\Program Files\Microsoft Visual Studio .NET 2003\VC7\bin"`
VBPATH=`cygpath -u "C:\Program Files\Microsoft Visual Studio .NET 2003\Common7\IDE"`
PATH="${PATH}:${VBPATH}:${LIBPATH}"
lib /DEF:${defname} /MACHINE:IX86 /OUT:${msvclibname}