[Trilinos-Users] Epetra on TACC's Stampede

[Gyorgy Matyasfalvi]

Dave, Mike, John, thanks for the fast responses!
Yes John I am using the default module load trilinos/... and  module load
intel/13.0.2.146. In terms of Trilinos I have tried both 11.4.2 and
10.12.2, I get the same error. With going what Mike said, I am wondering if
the issue could be that I'm loading a different BLAS then with what
Trilinos has been installed?
I think for starters I'm going to try to find out what BLAS was specified
during installation. If that's not the issue I'll try to install Trilinos.
Thanks for the script!
Best,
Gyorgy


On Sun, May 11, 2014 at 11:35 PM, John T. Foster <john.foster at utsa.edu>wrote:

> Gyorgy,
>
> I assume your using the Trilinos on Stampede that is by installed and
> accessed with module load trilinos command. While Epetra shouldn’t be an
> issue, I’ve run into other problems because the TACC Trilinos is not the
> complete set of libraries, so I’ve just been building Trilinos myself and
> have no problems. If you’d like to try building Trilinos yourself, I’m
> attaching my CMake config script for Stampede for your reference.
>
> JTF
>
>
> On Sun, May 11, 2014 at 6:41 PM, Gyorgy Matyasfalvi <matyasfalvi at gmail.com
> > wrote:
>
>> Dear User Community:
>>
>> I have already created a TACC ticket regarding this issue but I'm curious
>> if anyone has successfully used Epetra on Stampede?
>>
>> The code below segfaults when computing the sup-norm:
>>
>>
>> **************************************************************************************
>> #include <iostream>
>> #include "mpi.h"
>> #include "Epetra_MpiComm.h"
>> #include "Epetra_LocalMap.h"
>> #include "Epetra_Vector.h"
>> #include "Epetra_Version.h"
>> #include "mkl_cblas.h"
>>
>> #define N 1000
>>
>> int main(int argc, char* argv[]) {
>>
>> MPI_Init(&argc, &argv);
>>
>> Epetra_MpiComm comm(MPI_COMM_WORLD);
>>
>> std::cout<<Epetra_Version()<<std::endl<<std::endl<<"### TEST
>> ###"<<std::endl<<std::endl;
>>
>> Epetra_Map map(N,0,comm);
>>
>> std::cout<<"Map created"<<std::endl;
>>
>> Epetra_Vector x(map);
>>
>> std::cout<<"Vector created"<<std::endl;
>>
>> double norm;
>> x.NormInf(&norm);
>>
>> std::cout<<"sup-norm of x = "<<norm<<std::endl;
>>
>> MPI_Finalize();
>>
>> return 0;
>> }
>>
>>
>> **************************************************************************************
>>
>>
>> I get the following error message when debugging with ddt:
>>
>>
>> **************************************************************************************
>> Process 0:
>>
>> Memory error detected in mkl_blas_avx_idamax from
>> /opt/apps/intel/13/composer_xe_2013.2.146/mkl/lib/intel64/libmkl_avx.so:
>>
>> null pointer dereference or unaligned memory access
>>
>> Note: the latter may sometimes occur spuriously if guard pages are
>> enabled
>>
>> Tip: Use the stack list and the local variables to explore your program's
>> current state and identify the source of the error.
>>
>>
>> **************************************************************************************
>>
>> Thanks a lot!
>> Best,
>> Gyorgy
>>
>> _______________________________________________
>> Trilinos-Users mailing list
>> Trilinos-Users at software.sandia.gov
>> https://software.sandia.gov/mailman/listinfo/trilinos-users
>>
>>
>
>
> --
> John T. Foster
> Assistant Professor of Mechanical Engineering
> The University of Texas at San Antonio
> AET 2.330
> One UTSA Circle
> San Antonio, TX 78249
> 210.458.5521
> john.foster at utsa.edu
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://software.sandia.gov/pipermail/trilinos-users/attachments/20140512/0406bf9f/attachment.html>