configure

if test x$acx_fftw_ok = xyes; then

  LIBS="$FFTW_LIBS $LIBS"
	if test "$FFTW_WARN" != ""; then
	  { $as_echo "$as_me:${as_lineno-$LINENO}: WARNING: $FFTW_WARN" >&5
$as_echo "$as_me: WARNING: $FFTW_WARN" >&2;}
	fi
else

  as_fn_error $? "$FFTW_ERROR Exiting." "$LINENO" 5
fi


################## handle the ATLAS library(linear algebra) ##################




acx_atlas_ok=no
if test x$with_atlas_incdir = x; then
  if test x$with_atlas_libdir = x; then
    for ac_header in cblas.h clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

    if test x$acx_atlas_ok = xyes; then

$as_echo "#define ATLAS_BLAS_H \"cblas.h\"" >>confdefs.h


$as_echo "#define ATLAS_LAPACK_H \"clapack.h\"" >>confdefs.h

    else
      for ac_header in atlas/cblas.h atlas/clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

      if test x$acx_atlas_ok = xyes; then

$as_echo "#define ATLAS_BLAS_H \"atlas/cblas.h\"" >>confdefs.h


$as_echo "#define ATLAS_LAPACK_H \"atlas/clapack.h\"" >>confdefs.h

      else
        atlas_def=/usr/local/atlas
        for ac_header in $atlas_def/include/cblas.h $atlas_def/include/clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

        if test x$acx_atlas_ok = xyes; then

cat >>confdefs.h <<_ACEOF
#define ATLAS_BLAS_H "$atlas_def/include/cblas.h"
_ACEOF


cat >>confdefs.h <<_ACEOF
#define ATLAS_LAPACK_H "$atlas_def/include/clapack.h"
_ACEOF

        else
          atlas_def=/usr/atlas
          for ac_header in $atlas_def/include/cblas.h $atlas_def/include/clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

          if test x$acx_atlas_ok = xyes; then

cat >>confdefs.h <<_ACEOF
#define ATLAS_BLAS_H "$atlas_def/include/cblas.h"
_ACEOF


cat >>confdefs.h <<_ACEOF
#define ATLAS_LAPACK_H "$atlas_def/include/clapack.h"
_ACEOF

          else
            ATLAS_ERROR="CBLAS/LAPack include files not found!"
          fi
        fi
      fi
    fi
  else
    for ac_header in $with_atlas_libdir/include/cblas.h $with_atlas_libdir/include/clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

    if test x$acx_atlas_ok = xyes; then

cat >>confdefs.h <<_ACEOF
#define ATLAS_BLAS_H "$with_atlas_libdir/include/cblas.h"
_ACEOF


cat >>confdefs.h <<_ACEOF
#define ATLAS_LAPACK_H "$with_atlas_libdir/include/clapack.h"
_ACEOF

    else
      for ac_header in cblas.h clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

      if test x$acx_atlas_ok = xyes; then

cat >>confdefs.h <<_ACEOF
#define ATLAS_BLAS_H "cblas.h"
_ACEOF


cat >>confdefs.h <<_ACEOF
#define ATLAS_LAPACK_H "clapack.h"
_ACEOF

      else
        ATLAS_ERROR="CBLAS/LAPack include files not found in $with_atlas_libdir/include!"
      fi
    fi
  fi
else
  for ac_header in $with_atlas_incdir/cblas.h $with_atlas_incdir/clapack.h
do :
  as_ac_Header=`$as_echo "ac_cv_header_$ac_header" | $as_tr_sh`
ac_fn_c_check_header_mongrel "$LINENO" "$ac_header" "$as_ac_Header" "$ac_includes_default"
if eval test \"x\$"$as_ac_Header"\" = x"yes"; then :
  cat >>confdefs.h <<_ACEOF
#define `$as_echo "HAVE_$ac_header" | $as_tr_cpp` 1
_ACEOF
 acx_atlas_ok=yes
fi

done

  if test x$acx_atlas_ok = xyes; then

cat >>confdefs.h <<_ACEOF
#define ATLAS_BLAS_H "$with_atlas_incdir/cblas.h"
_ACEOF


cat >>confdefs.h <<_ACEOF
#define ATLAS_LAPACK_H "$with_atlas_incdir/clapack.h"
_ACEOF

  else
    ATLAS_ERROR="CBLAS/LAPack include files not found in $with_atlas_incdir!"
  fi
fi


if test x$acx_atlas_ok = xyes; then
  OLIBS="$LIBS"
  LIBS=""
  if test x$with_atlas_libdir = x; then
    { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

    { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

    if test x$acx_atlas_ok = xyes; then
      ATLAS_LIBPATH=""
    else
      atlas_def=/usr/local/atlas
      unset ac_cv_lib_lapack_clapack_dpotrf
      unset ac_cv_lib_cblas_cblas_dgemm
      acx_atlas_ok=yes
      { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$atlas_def/lib -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

      { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$atlas_def/lib -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

      if test x$acx_atlas_ok = xyes; then
        ATLAS_LIBPATH="-L$atlas_def/lib"
      else
        atlas_def=/usr/lib64/atlas
        unset ac_cv_lib_lapack_clapack_dpotrf
        unset ac_cv_lib_cblas_cblas_dgemm
        acx_atlas_ok=yes
        { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$atlas_def -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

        { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$atlas_def -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

        if test x$acx_atlas_ok = xyes; then
          ATLAS_LIBPATH="-L$atlas_def"
        else
          atlas_def=/usr/lib/atlas
          unset ac_cv_lib_lapack_clapack_dpotrf
          unset ac_cv_lib_cblas_cblas_dgemm
          acx_atlas_ok=yes
          { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$atlas_def -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

          { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$atlas_def -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

          if test x$acx_atlas_ok = xyes; then
            ATLAS_LIBPATH="-L$atlas_def"
          else
            atlas_def=/usr/atlas
            unset ac_cv_lib_lapack_clapack_dpotrf
            unset ac_cv_lib_cblas_cblas_dgemm
            acx_atlas_ok=yes
            { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$atlas_def/lib -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

            { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$atlas_def/lib -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

            if test x$acx_atlas_ok = xyes; then
              ATLAS_LIBPATH="-L$atlas_def/lib"
            else
              ATLAS_ERROR="CBLAS/LAPack library files not found at usual locations!"
            fi
          fi
        fi
      fi
    fi
  else
    { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$with_atlas_libdir -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

    { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$with_atlas_libdir -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

    if test x$acx_atlas_ok = xyes; then
      ATLAS_LIBPATH="-L$with_atlas_libdir"
    else
      unset ac_cv_lib_lapack_clapack_dpotrf
      unset ac_cv_lib_cblas_cblas_dgemm
      acx_atlas_ok=yes
      { $as_echo "$as_me:${as_lineno-$LINENO}: checking for clapack_dpotrf in -llapack" >&5
$as_echo_n "checking for clapack_dpotrf in -llapack... " >&6; }
if ${ac_cv_lib_lapack_clapack_dpotrf+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-llapack -L$with_atlas_libdir/lib -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char clapack_dpotrf ();
int
main ()
{
return clapack_dpotrf ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_lapack_clapack_dpotrf=yes
else
  ac_cv_lib_lapack_clapack_dpotrf=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_lapack_clapack_dpotrf" >&5
$as_echo "$ac_cv_lib_lapack_clapack_dpotrf" >&6; }
if test "x$ac_cv_lib_lapack_clapack_dpotrf" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBLAPACK 1
_ACEOF

  LIBS="-llapack $LIBS"

else
  acx_atlas_ok=no
fi

      { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lcblas... " >&6; }
if ${ac_cv_lib_cblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lcblas -L$with_atlas_libdir/lib -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_cblas_cblas_dgemm=yes
else
  ac_cv_lib_cblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_cblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_cblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_cblas_cblas_dgemm" = xyes; then :
  cat >>confdefs.h <<_ACEOF
#define HAVE_LIBCBLAS 1
_ACEOF

  LIBS="-lcblas $LIBS"

else
  acx_atlas_ok=no
fi

      if test x$acx_atlas_ok = xyes; then
        ATLAS_LIBPATH="-L$with_atlas_libdir/lib"
      else
        ATLAS_ERROR="CBLAS/LAPack library files not found in $with_atlas_libdir!"
      fi
    fi
  fi
  LIBS="$OLIBS"
fi




if test x"$acx_atlas_ok" = xyes; then

$as_echo "#define HAVE_ATLAS 1" >>confdefs.h

  if test x$use_pthreads = xyes; then
    { $as_echo "$as_me:${as_lineno-$LINENO}: checking for cblas_dgemm in -lptcblas" >&5
$as_echo_n "checking for cblas_dgemm in -lptcblas... " >&6; }
if ${ac_cv_lib_ptcblas_cblas_dgemm+:} false; then :
  $as_echo_n "(cached) " >&6
else
  ac_check_lib_save_LIBS=$LIBS
LIBS="-lptcblas $ATLAS_LIBPATH -lcblas -latlas -lm $LIBS"
cat confdefs.h - <<_ACEOF >conftest.$ac_ext
/* end confdefs.h.  */

/* Override any GCC internal prototype to avoid an error.
   Use char because int might match the return type of a GCC
   builtin and then its argument prototype would still apply.  */
#ifdef __cplusplus
extern "C"
#endif
char cblas_dgemm ();
int
main ()
{
return cblas_dgemm ();
  ;
  return 0;
}
_ACEOF
if ac_fn_c_try_link "$LINENO"; then :
  ac_cv_lib_ptcblas_cblas_dgemm=yes
else
  ac_cv_lib_ptcblas_cblas_dgemm=no
fi
rm -f core conftest.err conftest.$ac_objext \
    conftest$ac_exeext conftest.$ac_ext
LIBS=$ac_check_lib_save_LIBS
fi
{ $as_echo "$as_me:${as_lineno-$LINENO}: result: $ac_cv_lib_ptcblas_cblas_dgemm" >&5
$as_echo "$ac_cv_lib_ptcblas_cblas_dgemm" >&6; }
if test "x$ac_cv_lib_ptcblas_cblas_dgemm" = xyes; then :
  acx_atlast_ok=yes
else
  acx_atlast_ok=no
fi

    if test x$acx_atlast_ok = xyes; then
      ATLAS_LIBS="$ATLAS_LIBPATH -llapack -lptcblas -lcblas -latlas"
      LIBS="$OLIBS"


$as_echo "#define HAVE_ATLAS_MP 1" >>confdefs.h

      LIBS="$ATLAS_LIBS $LIBS"
	if test "$ATLAS_WARN" != ""; then
	  { $as_echo "$as_me:${as_lineno-$LINENO}: WARNING: $ATLAS_WARN" >&5
$as_echo "$as_me: WARNING: $ATLAS_WARN" >&2;}
	fi
    else
      ATLAS_LIBS="$ATLAS_LIBPATH -llapack -lcblas -latlas"
      LIBS="$OLIBS"

      ATLAS_WARN="CBLAS/LAPack was compiled without multithreading support!"

      LIBS="$ATLAS_LIBS $LIBS"
	if test "$ATLAS_WARN" != ""; then
	  { $as_echo "$as_me:${as_lineno-$LINENO}: WARNING: $ATLAS_WARN" >&5
$as_echo "$as_me: WARNING: $ATLAS_WARN" >&2;}
	fi
    fi
  else
    ATLAS_LIBS="$ATLAS_LIBPATH -llapack -lcblas -latlas"
    LIBS="$OLIBS"