From 84f481ef628aa5cf03460da4beddb7ccb970bbe4 Mon Sep 17 00:00:00 2001 From: Salvatore Filippone Date: Wed, 17 Jan 2007 15:08:50 +0000 Subject: [PATCH] Moving files from PSBLAS2 to MLD2P4. --- mld2p4/psb_umf_impl.c | 233 ------------------------------------------ 1 file changed, 233 deletions(-) delete mode 100644 mld2p4/psb_umf_impl.c diff --git a/mld2p4/psb_umf_impl.c b/mld2p4/psb_umf_impl.c deleted file mode 100644 index 902ce1a3..00000000 --- a/mld2p4/psb_umf_impl.c +++ /dev/null @@ -1,233 +0,0 @@ -/* - * MD2P4 - * Multilevel Domain Decomposition Parallel Preconditioner Package for PSBLAS - * for - * Parallel Sparse BLAS v2.0 - * - * (C) Copyright 2006 Salvatore Filippone University of Rome Tor Vergata - * Alfredo Buttari - * Daniela di Serafino Second University of Naples - * Pasqua D'Ambra ICAR-CNR - * - * Redistribution and use in source and binary forms, with or without - * modification, are permitted provided that the following conditions - * are met: - * 1. Redistributions of source code must retain the above copyright - * notice, this list of conditions and the following disclaimer. - * 2. Redistributions in binary form must reproduce the above copyright - * notice, this list of conditions, and the following disclaimer in the - * documentation and/or other materials provided with the distribution. - * 3. The name of the MD2P4 group or the names of its contributors may - * not be used to endorse or promote products derived from this - * software without specific written permission. - * - * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS - * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED - * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR - * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE MD2P4 GROUP OR ITS CONTRIBUTORS - * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR - * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF - * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS - * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN - * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) - * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE - * POSSIBILITY OF SUCH DAMAGE. - * - */ -/* This file is an interface to the UMFPACK routines for sparse - factorization. It was obtained by adapting umfpack_di_demo - under the original copyright terms reproduced below. - - PSBLAS v 2.0 */ - - -/* ===================== -UMFPACK Version 4.4 (Jan. 28, 2005), Copyright (c) 2005 by Timothy A. -Davis. All Rights Reserved. - -UMFPACK License: - - Your use or distribution of UMFPACK or any modified version of - UMFPACK implies that you agree to this License. - - THIS MATERIAL IS PROVIDED AS IS, WITH ABSOLUTELY NO WARRANTY - EXPRESSED OR IMPLIED. ANY USE IS AT YOUR OWN RISK. - - Permission is hereby granted to use or copy this program, provided - that the Copyright, this License, and the Availability of the original - version is retained on all copies. User documentation of any code that - uses UMFPACK or any modified version of UMFPACK code must cite the - Copyright, this License, the Availability note, and "Used by permission." - Permission to modify the code and to distribute modified code is granted, - provided the Copyright, this License, and the Availability note are - retained, and a notice that the code was modified is included. This - software was developed with support from the National Science Foundation, - and is provided to you free of charge. - -Availability: - - http://www.cise.ufl.edu/research/sparse/umfpack - -*/ - - - -#ifdef Add_ -#define psb_dumf_factor_ psb_dumf_factor_ -#define psb_dumf_solve_ psb_dumf_solve_ -#define psb_dumf_free_ psb_dumf_free_ -#endif -#ifdef AddDouble_ -#define psb_dumf_factor_ psb_dumf_factor__ -#define psb_dumf_solve_ psb_dumf_solve__ -#define psb_dumf_free_ psb_dumf_free__ -#endif -#ifdef NoChange -#define psb_dumf_factor_ psb_dumf_factor -#define psb_dumf_solve_ psb_dumf_solve -#define psb_dumf_free_ psb_dumf_free -#endif - - -#include -#ifdef Have_UMF_ -#include "umfpack.h" -#endif - -#ifdef LargeFptr -typedef long long fptr; /* 64-bit*/ -#else -typedef int fptr; /* 32-bit by default */ -#endif - -void -psb_dumf_factor_(int *n, int *nnz, - double *values, int *rowind, int *colptr, -#ifdef Have_UMF_ - fptr *symptr, - fptr *numptr, - -#else - void *symptr, - void *numptr, -#endif - int *info) - -{ - -#ifdef Have_UMF_ - double Info [UMFPACK_INFO], Control [UMFPACK_CONTROL]; - void *Symbolic, *Numeric ; - int i; - - - umfpack_di_defaults(Control); - - for (i = 0; i <= *n; ++i) --colptr[i]; - for (i = 0; i < *nnz; ++i) --rowind[i]; - *info = umfpack_di_symbolic (*n, *n, colptr, rowind, values, &Symbolic, - Control, Info); - - - if ( *info == UMFPACK_OK ) { - *info = 0; - } else { - printf("umfpack_di_symbolic() error returns INFO= %d\n", *info); - *info = -11; - *numptr = (fptr) NULL; - return; - } - - *symptr = (fptr) Symbolic; - - *info = umfpack_di_numeric (colptr, rowind, values, Symbolic, &Numeric, - Control, Info) ; - - - if ( *info == UMFPACK_OK ) { - *info = 0; - *numptr = (fptr) Numeric; - } else { - printf("umfpack_di_numeric() error returns INFO= %d\n", *info); - *info = -12; - *numptr = (fptr) NULL; - } - - for (i = 0; i <= *n; ++i) ++colptr[i]; - for (i = 0; i < *nnz; ++i) ++rowind[i]; -#else - fprintf(stderr," UMF Not Configured, fix make.inc and recompile\n"); - *info=-1; -#endif -} - - -void -psb_dumf_solve_(int *itrans, int *n, - double *x, double *b, int *ldb, -#ifdef Have_UMF_ - fptr *numptr, - -#else - void *numptr, -#endif - int *info) - -{ -#ifdef Have_UMF_ - double Info [UMFPACK_INFO], Control [UMFPACK_CONTROL]; - void *Symbolic, *Numeric ; - int i,trans; - - - umfpack_di_defaults(Control); - Control[UMFPACK_IRSTEP]=0; - - - if (*itrans == 0) { - trans = UMFPACK_A; - } else if (*itrans ==1) { - trans = UMFPACK_At; - } else { - trans = UMFPACK_A; - } - - *info = umfpack_di_solve(trans,NULL,NULL,NULL, - x,b,(void *) *numptr,Control,Info); - -#else - fprintf(stderr," UMF Not Configured, fix make.inc and recompile\n"); - *info=-1; -#endif - -} - - -void -psb_dumf_free_( -#ifdef Have_UMF_ - fptr *symptr, - fptr *numptr, - -#else - void *symptr, - void *numptr, -#endif - int *info) - -{ -#ifdef Have_UMF_ - void *Symbolic, *Numeric ; - Symbolic = (void *) *symptr; - Numeric = (void *) *numptr; - - umfpack_di_free_numeric(&Numeric); - umfpack_di_free_symbolic(&Symbolic); - *info=0; -#else - fprintf(stderr," UMF Not Configured, fix make.inc and recompile\n"); - *info=-1; -#endif -} - -