/usr/include/lobpcg.h is in libhypre-dev 2.8.0b-1build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 | /*BHEADER**********************************************************************
* Copyright (c) 2008, Lawrence Livermore National Security, LLC.
* Produced at the Lawrence Livermore National Laboratory.
* This file is part of HYPRE. See file COPYRIGHT for details.
*
* HYPRE is free software; you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License (as published by the Free
* Software Foundation) version 2.1 dated February 1999.
*
* $Revision: 2.7 $
***********************************************************************EHEADER*/
#include "multivector.h"
#ifndef LOCALLY_OPTIMAL_BLOCK_PRECONDITIONED_CONJUGATE_GRADIENTS
#define LOCALLY_OPTIMAL_BLOCK_PRECONDITIONED_CONJUGATE_GRADIENTS
#define PROBLEM_SIZE_TOO_SMALL 1
#define WRONG_BLOCK_SIZE 2
#define WRONG_CONSTRAINTS 3
#define REQUESTED_ACCURACY_NOT_ACHIEVED -1
typedef struct {
double absolute;
double relative;
} lobpcg_Tolerance;
typedef struct {
/* these pointers should point to 2 functions providing standard lapack functionality */
HYPRE_Int (*dpotrf) (char *uplo, HYPRE_Int *n, double *a, HYPRE_Int *
lda, HYPRE_Int *info);
HYPRE_Int (*dsygv) (HYPRE_Int *itype, char *jobz, char *uplo, HYPRE_Int *
n, double *a, HYPRE_Int *lda, double *b, HYPRE_Int *ldb,
double *w, double *work, HYPRE_Int *lwork, HYPRE_Int *info);
} lobpcg_BLASLAPACKFunctions;
#ifdef __cplusplus
extern "C" {
#endif
HYPRE_Int
lobpcg_solve( mv_MultiVectorPtr blockVectorX,
void* operatorAData,
void (*operatorA)( void*, void*, void* ),
void* operatorBData,
void (*operatorB)( void*, void*, void* ),
void* operatorTData,
void (*operatorT)( void*, void*, void* ),
mv_MultiVectorPtr blockVectorY,
lobpcg_BLASLAPACKFunctions blap_fn,
lobpcg_Tolerance tolerance,
HYPRE_Int maxIterations,
HYPRE_Int verbosityLevel,
HYPRE_Int* iterationNumber,
/* eigenvalues; "lambda_values" should point to array containing <blocksize> doubles where <blocksi
ze> is the width of multivector "blockVectorX" */
double * lambda_values,
/* eigenvalues history; a pointer to the entries of the <blocksize>-by-(<maxIterations>+1) matrix s
tored
in fortran-style. (i.e. column-wise) The matrix may be a submatrix of a larger matrix, see next
argument; If you don't need eigenvalues history, provide NULL in this entry */
double * lambdaHistory_values,
/* global height of the matrix (stored in fotran-style) specified by previous argument */
HYPRE_Int lambdaHistory_gh,
/* residual norms; argument should point to array of <blocksize> doubles */
double * residualNorms_values,
/* residual norms history; a pointer to the entries of the <blocksize>-by-(<maxIterations>+1) matri
x
stored in fortran-style. (i.e. column-wise) The matrix may be a submatrix of a larger matrix, see
next
argument If you don't need residual norms history, provide NULL in this entry */
double * residualNormsHistory_values ,
/* global height of the matrix (stored in fotran-style) specified by previous argument */
HYPRE_Int residualNormsHistory_gh
);
#ifdef __cplusplus
}
#endif
#endif /* LOCALLY_OPTIMAL_BLOCK_PRECONDITIONED_CONJUGATE_GRADIENTS */
|