Changeset a5b9608 in tol


Ignore:
Timestamp:
Jul 29, 2015, 8:03:06 PM (11 years ago)
Author:
Jorge Pérez <josp.jorge@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
b44a3b0
Parents:
4e57976
Message:

refs #1865, conversión de tipo erronea, ahora compila en linux tanto con clapack como lapacke

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@6795 be4846ae-f051-483a-831f-b64e7d08527f

Location:
tol
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • tol/bmath/blinalg/tol_lapack.cpp

    r4e57976 ra5b9608  
    6868  int n = S.Columns();
    6969  double* X_ = X.GetData().GetBuffer();
    70   int res = clapack_dpotrf(int(CblasRowMajor), LAPACK_UPLO(Uplo), n, X_, n);
     70  int res = clapack_dpotrf( CblasRowMajor, LAPACK_UPLO(Uplo), n, X_, n );
    7171  replicateTriang(X,Uplo,false);
    7272  return(res);
    … …  
    103103  int n = A.Columns();
    104104  double* X_ = X.GetData().GetBuffer();
    105   int res = clapack_dpotri( int(CblasRowMajor), LAPACK_UPLO(Uplo), n, X_, n);
     105  int res = clapack_dpotri( CblasRowMajor, LAPACK_UPLO(Uplo), n, X_, n );
    106106  replicateTriang(X,Uplo,true);
    107107  return(res);
  • tol/bmath/blinalg/tol_lapack.h

    r4e57976 ra5b9608  
    3232#undef I
    3333#define LAPACK_UPLO(UpLo) (UpLo==CblasLower?'L':'U')
     34#define LAPACK_ORDER(Order) (int)(order)
    3435#define clapack_dpotrf LAPACKE_dpotrf
    3536#define clapack_dpotri LAPACKE_dpotri
    … …  
    3839#include <clapack.h>
    3940#define LAPACK_UPLO(UpLo) UpLo
     41#define LAPACK_ORDER(Order) Order
    4042END_DECLS
    4143#endif
  • tol/btol/vmatrix_type/vmat_chol.cpp

    r4e57976 ra5b9608  
    3939  double* x = (double*)L.s_.blasRdense_->x;
    4040  memcpy(x, y, n*n*sizeof(double));
    41   res = clapack_dpotrf(CblasColMajor, LAPACK_UPLO(CblasLower), n, x, n);
     41  res = clapack_dpotrf( CblasColMajor, LAPACK_UPLO(CblasLower), n, x, n );
    4242  for(i=0; i<n; i++) for(j=i+1; j<n; j++) x[j*n+i] = 0;
    4343  isOk = (res==0);
Note: See TracChangeset for help on using the changeset viewer.