[petsc-users] Using STRUMPACK within PETSc
Satish Balay
balay at mcs.anl.gov
Mon Nov 16 11:21:25 CST 2015
>From STRUMPACK-Dense-1.1.1/src/StrumpackDensePackage_C.cpp
>>>>>>>>>
/* This C++ file implements the functions of the C interface. */
#include "StrumpackDensePackage.hpp"
extern "C" {
#include "StrumpackDensePackage.h"
<<<<<<<<<<
So looks like you have to use StrumpackDensePackage.h with "extern C". i.e mainSDP.cc should have:
#include <petsc.h>
extern "C" {
#include "StrumpackDensePackage.h"
}
[usually such details should be handled in the public include file - and not exposed to the user..]
Satish
On Mon, 16 Nov 2015, Niels Aage wrote:
>
> Hi all,
>
> First of all, I hope this is the right place to pose the following question. We're trying to create a shell preconditioner based on low rank properties of a linear system. For this to be efficient it requires a fast (mpi-parallelized) method for compressing (low rank) dense matrices. A tool that offers just this is STRUMPACK (http://portal.nersc.gov/project/sparse/strumpack/).
>
> So we simply need to link our petsc program to strumpack, but we have a hard time performing the linking. We have cooked our code down to a very simple test problem (see the attached code) that tries to invoke a single method from strumpack. This results in an "undefined reference to SDP_C_double_init(....)" even though the method clearly is in the object file "StrumpackDensePackage_C.o".
>
> Note that my PETSc (3.6.2) is build with scalapack support which is needed by strumpack - and that I have used the PETSc libs to build and link the strumpack test example "STRUMPACK-Dense-1.1.1/examples/c_example.c". And that this works just fine:-)
>
> We hope you have suggestions on how we can proceed - and let me know if any additional information is needed.
>
> Thanks,
> Niels Aage
>
> Associate Professor, Ph.D.
> Department of Mechanical Engineering,
> Section for Solid Mechanics
> Centre for Acoustic-Mechanical Micro Systems
> Technical University of Denmark, Building 404, DK-2800 Lyngby, Denmark
>
> Phone: (+45) 4525 4253, Fax: (+45) 4593 1475,
> E-mail: naage at mek.dtu.dk, Office: b404-122
> Group homepage: www.topopt.dtu.dk<http://www.topopt.dtu.dk>
> Centre homepage: www.camm.elektro.dtu.dk<http://www.camm.elektro.dtu.dk/>
>
More information about the petsc-users
mailing list