AnsweredAssumed Answered

ACML version 5.3.1: dsyev slows down dramatically after dpotrf

Question asked by mrip on Sep 20, 2013
Latest reply on Sep 27, 2013 by mrip

Running ACML version 5.3.1, libacml_mp using fma4 on Opteron 6348 processors, Ubuntu 12.04.


Calls to dsyev (eigen decomposition) slow down dramatically -- by a factor of 10 -- after calling dpotrf (cholesky decomposition).

Here is a simple C program that reproduces the problem:


#include <stdio.h>

#include <stdlib.h>

#include <acml.h>

#include <time.h>


int main(void) {

  double * x = malloc(1000000 * sizeof(double));

  double * y = malloc(1000000 * sizeof(double));

  double * eig0 = malloc(1000000 * sizeof(double));

  double * eig1 = malloc(1000000 * sizeof(double));

  double * eigw = malloc(1000 * sizeof(double));

  double * chol = malloc(1000000 * sizeof(double));


  clock_t t0,t1;

  int info;

  int i;


  // generate a random matrix

  for(i = 0; i<1000000; ++i){

    x[i] = rand() / (double) RAND_MAX;


  // compute y = xx^T so that y is symmetric positive definite



  // make a copy of y for cholesky and eigen decompositions

  for(i = 0; i<1000000; ++i){

    chol[i] = y[i];

    eig0[i] = y[i];

    eig1[i] = y[i];



  // first eigenvalue test

  t0 = clock();


  t1 = clock();

  printf("Eigen decomposition time: %d\n", (t1-t0)/1000);


  // cholesky



  // second eigenvalue test, after cholesky

  t0 = clock();


  t1 = clock();

  printf("Eigen decomposition time: %d\n", (t1-t0)/1000);



Here is the output:

Eigen decomposition time: 8120

Eigen decomposition time: 95140


If I comment out the dpotrf line, then it works fine:

Eigen decomposition time: 8150

Eigen decomposition time: 8210


This seems like some kind of a bug.  Am I missing something?  Is there some kind of cache that I need to clear?




On edit: I get the same behavior whether I link against the gfortran64_mp or the gfortran64_fma4_mp versions of libacml_mp.