17758a8cdSBarry Smith 27758a8cdSBarry Smith #if !defined(__TFS_H) 37758a8cdSBarry Smith #define __TFS_H 47758a8cdSBarry Smith 57758a8cdSBarry Smith /**********************************const.h************************************* 67758a8cdSBarry Smith 77758a8cdSBarry Smith Author: Henry M. Tufo III 87758a8cdSBarry Smith 97758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 107758a8cdSBarry Smith 117758a8cdSBarry Smith snail-mail: 127758a8cdSBarry Smith Division of Applied Mathematics 137758a8cdSBarry Smith Brown University 147758a8cdSBarry Smith Providence, RI 02912 157758a8cdSBarry Smith 167758a8cdSBarry Smith Last Modification: 177758a8cdSBarry Smith 6.21.97 187758a8cdSBarry Smith ***********************************const.h************************************/ 197758a8cdSBarry Smith 207758a8cdSBarry Smith /**********************************const.h************************************* 217758a8cdSBarry Smith File Description: 227758a8cdSBarry Smith ----------------- 237758a8cdSBarry Smith 247758a8cdSBarry Smith ***********************************const.h************************************/ 257758a8cdSBarry Smith #include "petsc.h" 26330ea6edSBarry Smith #include "petscsys.h" 27330ea6edSBarry Smith #if defined(PETSC_HAVE_STDLIB_H) 28330ea6edSBarry Smith #include <stdlib.h> 29330ea6edSBarry Smith #endif 309e8a9152SBarry Smith #if defined(PETSC_HAVE_MALLOC_H) 319e8a9152SBarry Smith #include <malloc.h> 329e8a9152SBarry Smith #endif 33330ea6edSBarry Smith #include "petscblaslapack.h" 34330ea6edSBarry Smith #include <limits.h> 35330ea6edSBarry Smith #include <float.h> 367758a8cdSBarry Smith 377758a8cdSBarry Smith #define X 0 387758a8cdSBarry Smith #define Y 1 397758a8cdSBarry Smith #define Z 2 407758a8cdSBarry Smith #define XY 3 417758a8cdSBarry Smith #define XZ 4 427758a8cdSBarry Smith #define YZ 5 437758a8cdSBarry Smith 447758a8cdSBarry Smith 457758a8cdSBarry Smith #define THRESH 0.2 467758a8cdSBarry Smith #define N_HALF 4096 477758a8cdSBarry Smith #define PRIV_BUF_SZ 45 487758a8cdSBarry Smith 497758a8cdSBarry Smith /*4096 8192 32768 65536 1048576 */ 507758a8cdSBarry Smith #define MAX_MSG_BUF 32768 517758a8cdSBarry Smith 527758a8cdSBarry Smith #define FULL 2 537758a8cdSBarry Smith #define PARTIAL 1 547758a8cdSBarry Smith #define NONE 0 557758a8cdSBarry Smith 567758a8cdSBarry Smith #define BYTE 8 577758a8cdSBarry Smith #define BIT_0 0x1 587758a8cdSBarry Smith #define BIT_1 0x2 597758a8cdSBarry Smith #define BIT_2 0x4 607758a8cdSBarry Smith #define BIT_3 0x8 617758a8cdSBarry Smith #define BIT_4 0x10 627758a8cdSBarry Smith #define BIT_5 0x20 637758a8cdSBarry Smith #define BIT_6 0x40 647758a8cdSBarry Smith #define BIT_7 0x80 657758a8cdSBarry Smith #define TOP_BIT INT_MIN 667758a8cdSBarry Smith #define ALL_ONES -1 677758a8cdSBarry Smith 687758a8cdSBarry Smith #define FALSE 0 697758a8cdSBarry Smith #define TRUE 1 707758a8cdSBarry Smith 717758a8cdSBarry Smith #define C 0 727758a8cdSBarry Smith 737758a8cdSBarry Smith 747758a8cdSBarry Smith #define MAX_VEC 1674 757758a8cdSBarry Smith #define FORMAT 30 767758a8cdSBarry Smith #define MAX_COL_LEN 100 777758a8cdSBarry Smith #define MAX_LINE FORMAT*MAX_COL_LEN 787758a8cdSBarry Smith #define DELIM " \n \t" 797758a8cdSBarry Smith #define LINE 12 807758a8cdSBarry Smith #define C_LINE 80 817758a8cdSBarry Smith 827758a8cdSBarry Smith #define REAL_MAX DBL_MAX 837758a8cdSBarry Smith #define REAL_MIN DBL_MIN 847758a8cdSBarry Smith 857758a8cdSBarry Smith #define UT 5 /* dump upper 1/2 */ 867758a8cdSBarry Smith #define LT 6 /* dump lower 1/2 */ 877758a8cdSBarry Smith #define SYMM 8 /* we assume symm and dump upper 1/2 */ 887758a8cdSBarry Smith #define NON_SYMM 9 897758a8cdSBarry Smith 907758a8cdSBarry Smith #define ROW 10 917758a8cdSBarry Smith #define COL 11 927758a8cdSBarry Smith 937758a8cdSBarry Smith #define EPS 1.0e-14 947758a8cdSBarry Smith #define EPS2 1.0e-07 957758a8cdSBarry Smith 967758a8cdSBarry Smith 977758a8cdSBarry Smith #define MPI 1 987758a8cdSBarry Smith #define NX 2 997758a8cdSBarry Smith 1007758a8cdSBarry Smith #define LOG2(x) (PetscScalar)log((double)x)/log(2) 1017758a8cdSBarry Smith #define SWAP(a,b) temp=(a); (a)=(b); (b)=temp; 1027758a8cdSBarry Smith #define P_SWAP(a,b) ptr=(a); (a)=(b); (b)=ptr; 1037758a8cdSBarry Smith 1047758a8cdSBarry Smith #define MAX_FABS(x,y) ((double)fabs(x)>(double)fabs(y)) ? ((PetscScalar)x) : ((PetscScalar)y) 1057758a8cdSBarry Smith #define MIN_FABS(x,y) ((double)fabs(x)<(double)fabs(y)) ? ((PetscScalar)x) : ((PetscScalar)y) 1067758a8cdSBarry Smith 1077758a8cdSBarry Smith /* specer's existence ... can be done w/MAX_ABS */ 1087758a8cdSBarry Smith #define EXISTS(x,y) ((x)==0.0) ? (y) : (x) 1097758a8cdSBarry Smith 1107758a8cdSBarry Smith #define MULT_NEG_ONE(a) (a) *= -1; 1117758a8cdSBarry Smith #define NEG(a) (a) |= BIT_31; 1127758a8cdSBarry Smith #define POS(a) (a) &= INT_MAX; 1137758a8cdSBarry Smith 1147758a8cdSBarry Smith 1157758a8cdSBarry Smith 1167758a8cdSBarry Smith 1177758a8cdSBarry Smith /**********************************types.h************************************* 1187758a8cdSBarry Smith 1197758a8cdSBarry Smith Author: Henry M. Tufo III 1207758a8cdSBarry Smith 1217758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 1227758a8cdSBarry Smith 1237758a8cdSBarry Smith snail-mail: 1247758a8cdSBarry Smith Division of Applied Mathematics 1257758a8cdSBarry Smith Brown University 1267758a8cdSBarry Smith Providence, RI 02912 1277758a8cdSBarry Smith 1287758a8cdSBarry Smith Last Modification: 1297758a8cdSBarry Smith 6.21.97 1307758a8cdSBarry Smith ***********************************types.h************************************/ 1317758a8cdSBarry Smith 1327758a8cdSBarry Smith /**********************************types.h************************************* 1337758a8cdSBarry Smith File Description: 1347758a8cdSBarry Smith ----------------- 1357758a8cdSBarry Smith 1367758a8cdSBarry Smith ***********************************types.h************************************/ 1373fdc5746SBarry Smith typedef PetscErrorCode (*vfp)(void*,void*,int,...); 138*52f87cdaSBarry Smith typedef PetscErrorCode (*rbfp)(PetscScalar *, PetscScalar *, PetscInt len); 1397758a8cdSBarry Smith #define vbfp MPI_User_function * 140*52f87cdaSBarry Smith typedef PetscInt (*bfp)(void*, void *, PetscInt *len, MPI_Datatype *dt); 1417758a8cdSBarry Smith 1427758a8cdSBarry Smith /***********************************comm.h************************************* 1437758a8cdSBarry Smith 1447758a8cdSBarry Smith Author: Henry M. Tufo III 1457758a8cdSBarry Smith 1467758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 1477758a8cdSBarry Smith 1487758a8cdSBarry Smith snail-mail: 1497758a8cdSBarry Smith Division of Applied Mathematics 1507758a8cdSBarry Smith Brown University 1517758a8cdSBarry Smith Providence, RI 02912 1527758a8cdSBarry Smith 1537758a8cdSBarry Smith Last Modification: 1547758a8cdSBarry Smith 6.21.97 1557758a8cdSBarry Smith ***********************************comm.h*************************************/ 1567758a8cdSBarry Smith 1577758a8cdSBarry Smith /***********************************comm.h************************************* 1587758a8cdSBarry Smith File Description: 1597758a8cdSBarry Smith ----------------- 1607758a8cdSBarry Smith 1617758a8cdSBarry Smith ***********************************comm.h*************************************/ 1627758a8cdSBarry Smith 1637758a8cdSBarry Smith /***********************************comm.h************************************* 1647758a8cdSBarry Smith Function: 1657758a8cdSBarry Smith 1667758a8cdSBarry Smith Input : 1677758a8cdSBarry Smith Output: 1687758a8cdSBarry Smith Return: 1697758a8cdSBarry Smith Description: 1707758a8cdSBarry Smith Usage: 1717758a8cdSBarry Smith ***********************************comm.h*************************************/ 1723fdc5746SBarry Smith extern PetscMPIInt my_id; 1733fdc5746SBarry Smith extern PetscMPIInt num_nodes; 1743fdc5746SBarry Smith extern PetscMPIInt floor_num_nodes; 1753fdc5746SBarry Smith extern PetscMPIInt i_log2_num_nodes; 1767758a8cdSBarry Smith 177*52f87cdaSBarry Smith extern PetscErrorCode giop(PetscInt *vals, PetscInt *work, PetscInt n, PetscInt *oprs); 178*52f87cdaSBarry Smith extern PetscErrorCode grop(PetscScalar *vals, PetscScalar *work, PetscInt n, PetscInt *oprs); 179*52f87cdaSBarry Smith extern PetscErrorCode gfop(void *vals, void *wk, PetscInt n, vbfp fp, MPI_Datatype dt, PetscInt comm_type); 1803fdc5746SBarry Smith extern PetscErrorCode comm_init(void); 181*52f87cdaSBarry Smith extern PetscErrorCode giop_hc(PetscInt *vals, PetscInt *work, PetscInt n, PetscInt *oprs, PetscInt dim); 182*52f87cdaSBarry Smith extern PetscErrorCode grop_hc(PetscScalar *vals, PetscScalar *work, PetscInt n, PetscInt *oprs, PetscInt dim); 183*52f87cdaSBarry Smith extern PetscErrorCode ssgl_radd(PetscScalar *vals, PetscScalar *work, PetscInt level, PetscInt *segs); 1847758a8cdSBarry Smith 1857758a8cdSBarry Smith #define MSGTAG0 101 1867758a8cdSBarry Smith #define MSGTAG1 1001 1877758a8cdSBarry Smith #define MSGTAG2 76207 1887758a8cdSBarry Smith #define MSGTAG3 100001 1897758a8cdSBarry Smith #define MSGTAG4 163841 1907758a8cdSBarry Smith #define MSGTAG5 249439 1917758a8cdSBarry Smith #define MSGTAG6 10000001 1927758a8cdSBarry Smith 1937758a8cdSBarry Smith 1947758a8cdSBarry Smith /**********************************error.h************************************* 1957758a8cdSBarry Smith 1967758a8cdSBarry Smith Author: Henry M. Tufo III 1977758a8cdSBarry Smith 1987758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 1997758a8cdSBarry Smith 2007758a8cdSBarry Smith snail-mail: 2017758a8cdSBarry Smith Division of Applied Mathematics 2027758a8cdSBarry Smith Brown University 2037758a8cdSBarry Smith Providence, RI 02912 2047758a8cdSBarry Smith 2057758a8cdSBarry Smith Last Modification: 2067758a8cdSBarry Smith 6.21.97 2077758a8cdSBarry Smith **********************************error.h*************************************/ 2087758a8cdSBarry Smith 2097758a8cdSBarry Smith /**********************************error.h************************************* 2107758a8cdSBarry Smith File Description: 2117758a8cdSBarry Smith ----------------- 2127758a8cdSBarry Smith 2137758a8cdSBarry Smith **********************************error.h*************************************/ 2147758a8cdSBarry Smith 2157758a8cdSBarry Smith /**********************************error.h************************************* 2167758a8cdSBarry Smith Function: error_msg_fatal() 2177758a8cdSBarry Smith 2187758a8cdSBarry Smith Input : formatted string and arguments. 2197758a8cdSBarry Smith Output: conversion printed to stdout. 2207758a8cdSBarry Smith Return: na. 2217758a8cdSBarry Smith Description: prints error message and terminates program. 2227758a8cdSBarry Smith Usage: error_msg_fatal("this is my %d'st test",test_num) 2237758a8cdSBarry Smith **********************************error.h*************************************/ 2243fdc5746SBarry Smith extern PetscErrorCode error_msg_fatal(const char msg[], ...); 2257758a8cdSBarry Smith 2267758a8cdSBarry Smith 2277758a8cdSBarry Smith 2287758a8cdSBarry Smith /**********************************error.h************************************* 2297758a8cdSBarry Smith Function: error_msg_warning() 2307758a8cdSBarry Smith 2317758a8cdSBarry Smith Input : formatted string and arguments. 2327758a8cdSBarry Smith Output: conversion printed to stdout. 2337758a8cdSBarry Smith Return: na. 2347758a8cdSBarry Smith Description: prints error message. 2357758a8cdSBarry Smith Usage: error_msg_warning("this is my %d'st test",test_num) 2367758a8cdSBarry Smith **********************************error.h*************************************/ 2373fdc5746SBarry Smith extern PetscErrorCode error_msg_warning(const char msg[], ...); 2387758a8cdSBarry Smith 2397758a8cdSBarry Smith /*$Id: vector.c,v 1.228 2001/03/23 23:21:22 balay Exp $*/ 2407758a8cdSBarry Smith /**********************************ivec.h************************************** 2417758a8cdSBarry Smith 2427758a8cdSBarry Smith Author: Henry M. Tufo III 2437758a8cdSBarry Smith 2447758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 2457758a8cdSBarry Smith 2467758a8cdSBarry Smith snail-mail: 2477758a8cdSBarry Smith Division of Applied Mathematics 2487758a8cdSBarry Smith Brown University 2497758a8cdSBarry Smith Providence, RI 02912 2507758a8cdSBarry Smith 2517758a8cdSBarry Smith Last Modification: 2527758a8cdSBarry Smith 6.21.97 2537758a8cdSBarry Smith ***********************************ivec.h*************************************/ 2547758a8cdSBarry Smith 2557758a8cdSBarry Smith /**********************************ivec.h************************************** 2567758a8cdSBarry Smith File Description: 2577758a8cdSBarry Smith ----------------- 2587758a8cdSBarry Smith 2597758a8cdSBarry Smith ***********************************ivec.h*************************************/ 2607758a8cdSBarry Smith 2617758a8cdSBarry Smith #define SORT_REAL 1 2627758a8cdSBarry Smith #define SORT_INTEGER 0 2637758a8cdSBarry Smith #define SORT_INT_PTR 2 2647758a8cdSBarry Smith 2657758a8cdSBarry Smith 2667758a8cdSBarry Smith #define NON_UNIFORM 0 2677758a8cdSBarry Smith #define GL_MAX 1 2687758a8cdSBarry Smith #define GL_MIN 2 2697758a8cdSBarry Smith #define GL_MULT 3 2707758a8cdSBarry Smith #define GL_ADD 4 2717758a8cdSBarry Smith #define GL_B_XOR 5 2727758a8cdSBarry Smith #define GL_B_OR 6 2737758a8cdSBarry Smith #define GL_B_AND 7 2747758a8cdSBarry Smith #define GL_L_XOR 8 2757758a8cdSBarry Smith #define GL_L_OR 9 2767758a8cdSBarry Smith #define GL_L_AND 10 2777758a8cdSBarry Smith #define GL_MAX_ABS 11 2787758a8cdSBarry Smith #define GL_MIN_ABS 12 2797758a8cdSBarry Smith #define GL_EXISTS 13 2807758a8cdSBarry Smith 2817758a8cdSBarry Smith 2827758a8cdSBarry Smith 2837758a8cdSBarry Smith /**********************************ivec.h************************************** 2847758a8cdSBarry Smith Function: 2857758a8cdSBarry Smith 2867758a8cdSBarry Smith Input : 2877758a8cdSBarry Smith Output: 2887758a8cdSBarry Smith Return: 2897758a8cdSBarry Smith Description: 2907758a8cdSBarry Smith Usage: 2917758a8cdSBarry Smith ***********************************ivec.h*************************************/ 292*52f87cdaSBarry Smith extern PetscInt *ivec_copy(PetscInt *arg1, PetscInt *arg2, PetscInt n); 2937758a8cdSBarry Smith 294*52f87cdaSBarry Smith extern PetscErrorCode ivec_zero(PetscInt *arg1, PetscInt n); 295*52f87cdaSBarry Smith extern PetscErrorCode ivec_set(PetscInt *arg1, PetscInt arg2, PetscInt n); 2967758a8cdSBarry Smith 297*52f87cdaSBarry Smith extern PetscInt ivec_lb(PetscInt *work, PetscInt n); 298*52f87cdaSBarry Smith extern PetscInt ivec_ub(PetscInt *work, PetscInt n); 299*52f87cdaSBarry Smith extern PetscInt ivec_sum(PetscInt *arg1, PetscInt n); 3007758a8cdSBarry Smith 301*52f87cdaSBarry Smith extern vfp ivec_fct_addr(PetscInt type); 3027758a8cdSBarry Smith 303*52f87cdaSBarry Smith extern PetscErrorCode ivec_non_uniform(PetscInt *arg1, PetscInt *arg2, PetscInt n, PetscInt *arg3); 304*52f87cdaSBarry Smith extern PetscErrorCode ivec_max(PetscInt *arg1, PetscInt *arg2, PetscInt n); 305*52f87cdaSBarry Smith extern PetscErrorCode ivec_min(PetscInt *arg1, PetscInt *arg2, PetscInt n); 306*52f87cdaSBarry Smith extern PetscErrorCode ivec_mult(PetscInt *arg1, PetscInt *arg2, PetscInt n); 307*52f87cdaSBarry Smith extern PetscErrorCode ivec_add(PetscInt *arg1, PetscInt *arg2, PetscInt n); 308*52f87cdaSBarry Smith extern PetscErrorCode ivec_xor(PetscInt *arg1, PetscInt *arg2, PetscInt n); 309*52f87cdaSBarry Smith extern PetscErrorCode ivec_or(PetscInt *arg1, PetscInt *arg2, PetscInt len); 310*52f87cdaSBarry Smith extern PetscErrorCode ivec_and(PetscInt *arg1, PetscInt *arg2, PetscInt len); 311*52f87cdaSBarry Smith extern PetscErrorCode ivec_lxor(PetscInt *arg1, PetscInt *arg2, PetscInt n); 312*52f87cdaSBarry Smith extern PetscErrorCode ivec_lor(PetscInt *arg1, PetscInt *arg2, PetscInt len); 313*52f87cdaSBarry Smith extern PetscErrorCode ivec_land(PetscInt *arg1, PetscInt *arg2, PetscInt len); 314*52f87cdaSBarry Smith extern PetscErrorCode ivec_and3( PetscInt *arg1, PetscInt *arg2, PetscInt *arg3, PetscInt n); 3157758a8cdSBarry Smith 316*52f87cdaSBarry Smith extern PetscErrorCode ivec_sort_companion(PetscInt *ar, PetscInt *ar2, PetscInt size); 317*52f87cdaSBarry Smith extern PetscErrorCode ivec_sort(PetscInt *ar, PetscInt size); 318*52f87cdaSBarry Smith extern PetscErrorCode SMI_sort(void *ar1, void *ar2, PetscInt size, PetscInt type); 319*52f87cdaSBarry Smith extern PetscInt ivec_binary_search(PetscInt item, PetscInt *list, PetscInt n); 320*52f87cdaSBarry Smith extern PetscInt ivec_linear_search(PetscInt item, PetscInt *list, PetscInt n); 3217758a8cdSBarry Smith 322*52f87cdaSBarry Smith extern PetscErrorCode ivec_sort_companion_hack(PetscInt *ar, PetscInt **ar2, PetscInt size); 3237758a8cdSBarry Smith 3247758a8cdSBarry Smith 325*52f87cdaSBarry Smith extern PetscErrorCode rvec_zero(PetscScalar *arg1, PetscInt n); 326*52f87cdaSBarry Smith extern PetscErrorCode rvec_one(PetscScalar *arg1, PetscInt n); 327*52f87cdaSBarry Smith extern PetscErrorCode rvec_set(PetscScalar *arg1, PetscScalar arg2, PetscInt n); 328*52f87cdaSBarry Smith extern PetscErrorCode rvec_copy(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 329*52f87cdaSBarry Smith extern PetscErrorCode rvec_scale(PetscScalar *arg1, PetscScalar arg2, PetscInt n); 3307758a8cdSBarry Smith 331*52f87cdaSBarry Smith extern vfp rvec_fct_addr(PetscInt type); 332*52f87cdaSBarry Smith extern PetscErrorCode rvec_add(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 333*52f87cdaSBarry Smith extern PetscErrorCode rvec_mult(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 334*52f87cdaSBarry Smith extern PetscErrorCode rvec_max(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 335*52f87cdaSBarry Smith extern PetscErrorCode rvec_max_abs(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 336*52f87cdaSBarry Smith extern PetscErrorCode rvec_min(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 337*52f87cdaSBarry Smith extern PetscErrorCode rvec_min_abs(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 338*52f87cdaSBarry Smith extern PetscErrorCode vec_exists(PetscScalar *arg1, PetscScalar *arg2, PetscInt n); 3397758a8cdSBarry Smith 3407758a8cdSBarry Smith 3417758a8cdSBarry Smith /*$Id: vector.c,v 1.228 2001/03/23 23:21:22 balay Exp $*/ 3427758a8cdSBarry Smith /***********************************gs.h*************************************** 3437758a8cdSBarry Smith 3447758a8cdSBarry Smith Author: Henry M. Tufo III 3457758a8cdSBarry Smith 3467758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 3477758a8cdSBarry Smith 3487758a8cdSBarry Smith snail-mail: 3497758a8cdSBarry Smith Division of Applied Mathematics 3507758a8cdSBarry Smith Brown University 3517758a8cdSBarry Smith Providence, RI 02912 3527758a8cdSBarry Smith 3537758a8cdSBarry Smith Last Modification: 3547758a8cdSBarry Smith 6.21.97 3557758a8cdSBarry Smith ************************************gs.h**************************************/ 3567758a8cdSBarry Smith 3577758a8cdSBarry Smith /***********************************gs.h*************************************** 3587758a8cdSBarry Smith File Description: 3597758a8cdSBarry Smith ----------------- 3607758a8cdSBarry Smith 3617758a8cdSBarry Smith ************************************gs.h**************************************/ 3627758a8cdSBarry Smith 3637758a8cdSBarry Smith /***********************************gs.h*************************************** 3647758a8cdSBarry Smith Type: gs_ADT 3657758a8cdSBarry Smith ------------ 3667758a8cdSBarry Smith 3677758a8cdSBarry Smith ************************************gs.h**************************************/ 3687758a8cdSBarry Smith 3697758a8cdSBarry Smith typedef struct gather_scatter_id *gs_ADT; 370*52f87cdaSBarry Smith typedef PetscErrorCode (*Rbfp)(PetscScalar *, PetscScalar *, PetscInt len); 3717758a8cdSBarry Smith 3727758a8cdSBarry Smith /***********************************gs.h*************************************** 3737758a8cdSBarry Smith Function: 3747758a8cdSBarry Smith 3757758a8cdSBarry Smith Input : 3767758a8cdSBarry Smith Output: 3777758a8cdSBarry Smith Return: 3787758a8cdSBarry Smith Description: 3797758a8cdSBarry Smith Usage: 3807758a8cdSBarry Smith ************************************gs.h**************************************/ 381*52f87cdaSBarry Smith extern gs_ADT gs_init(PetscInt *elms, PetscInt nel, PetscInt level); 3823fdc5746SBarry Smith extern PetscErrorCode gs_gop(gs_ADT gs_handle, PetscScalar *vals, const char *op); 383*52f87cdaSBarry Smith extern PetscErrorCode gs_gop_vec(gs_ADT gs_handle, PetscScalar *vals, const char *op, PetscInt step); 3843fdc5746SBarry Smith extern PetscErrorCode gs_gop_binary(gs_ADT gs, PetscScalar *vals, Rbfp fct); 385*52f87cdaSBarry Smith extern PetscErrorCode gs_gop_hc(gs_ADT gs_handle, PetscScalar *vals, const char *op, PetscInt dim); 3863fdc5746SBarry Smith extern PetscErrorCode gs_free(gs_ADT gs_handle); 387*52f87cdaSBarry Smith extern PetscErrorCode gs_init_msg_buf_sz(PetscInt buf_size); 388*52f87cdaSBarry Smith extern PetscErrorCode gs_init_vec_sz(PetscInt size); 3897758a8cdSBarry Smith 3907758a8cdSBarry Smith 3917758a8cdSBarry Smith 3927758a8cdSBarry Smith /*************************************xxt.h************************************ 3937758a8cdSBarry Smith Module Name: xxt 3947758a8cdSBarry Smith Module Info: need xxt.{c,h} gs.{c,h} comm.{c,h} ivec.{c,h} error.{c,h} 3957758a8cdSBarry Smith 3967758a8cdSBarry Smith author: Henry M. Tufo III 3977758a8cdSBarry Smith e-mail: hmt@asci.uchicago.edu 3987758a8cdSBarry Smith contact: 3997758a8cdSBarry Smith +--------------------------------+--------------------------------+ 4007758a8cdSBarry Smith |MCS Division - Building 221 |Department of Computer Science | 4017758a8cdSBarry Smith |Argonne National Laboratory |Ryerson 152 | 4027758a8cdSBarry Smith |9700 S. Cass Avenue |The University of Chicago | 4037758a8cdSBarry Smith |Argonne, IL 60439 |Chicago, IL 60637 | 4047758a8cdSBarry Smith |(630) 252-5354/5986 ph/fx |(773) 702-6019/8487 ph/fx | 4057758a8cdSBarry Smith +--------------------------------+--------------------------------+ 4067758a8cdSBarry Smith 4077758a8cdSBarry Smith Last Modification: 3.20.01 4087758a8cdSBarry Smith **************************************xxt.h***********************************/ 4097758a8cdSBarry Smith 4107758a8cdSBarry Smith /*************************************xxt.h************************************ 4117758a8cdSBarry Smith File Description: 4127758a8cdSBarry Smith **************************************xxt.h***********************************/ 4137758a8cdSBarry Smith 4147758a8cdSBarry Smith /*************************************xxt.h************************************ 4157758a8cdSBarry Smith Notes on Usage: 4167758a8cdSBarry Smith **************************************xxt.h***********************************/ 4177758a8cdSBarry Smith 4187758a8cdSBarry Smith 4197758a8cdSBarry Smith typedef struct xxt_CDT *xxt_ADT; 4207758a8cdSBarry Smith 4217758a8cdSBarry Smith 4227758a8cdSBarry Smith /*************************************xxt.h************************************ 4237758a8cdSBarry Smith Function: XXT_new() 4247758a8cdSBarry Smith 4257758a8cdSBarry Smith Input : 4267758a8cdSBarry Smith Output: 4277758a8cdSBarry Smith Return: ADT ptr or NULL upon failure. 4287758a8cdSBarry Smith Description: This function allocates and returns an xxt handle 4297758a8cdSBarry Smith Usage: xxt_handle = xxt_new(); 4307758a8cdSBarry Smith **************************************xxt.h***********************************/ 4317758a8cdSBarry Smith extern xxt_ADT XXT_new(void); 4327758a8cdSBarry Smith 4337758a8cdSBarry Smith 4347758a8cdSBarry Smith /*************************************xxt.h************************************ 4357758a8cdSBarry Smith Function: XXT_free() 4367758a8cdSBarry Smith 4377758a8cdSBarry Smith Input : pointer to ADT. 4387758a8cdSBarry Smith Output: 4397758a8cdSBarry Smith Return: 4407758a8cdSBarry Smith Description: This function frees the storage associated with an xxt handle 4417758a8cdSBarry Smith Usage: XXT_free(xxt_handle); 4427758a8cdSBarry Smith **************************************xxt.h***********************************/ 443*52f87cdaSBarry Smith EXTERN PetscInt XXT_free(xxt_ADT xxt_handle); 4447758a8cdSBarry Smith 4457758a8cdSBarry Smith 4467758a8cdSBarry Smith /*************************************xxt.h************************************ 4477758a8cdSBarry Smith Function: XXT_factor 4487758a8cdSBarry Smith 4497758a8cdSBarry Smith Input : ADT ptr, and pointer to object 4507758a8cdSBarry Smith Output: 4517758a8cdSBarry Smith Return: 0 on failure, 1 on success 4527758a8cdSBarry Smith Description: This function sets the xxt solver 4537758a8cdSBarry Smith 4547758a8cdSBarry Smith xxt assumptions: given n rows of global coarse matrix (E_loc) where 4557758a8cdSBarry Smith o global dofs N = sum_p(n), p=0,P-1 4567758a8cdSBarry Smith (i.e. row dist. with no dof replication) 4577758a8cdSBarry Smith (5.21.00 will handle dif replication case) 4587758a8cdSBarry Smith o m is the number of columns in E_loc (m>=n) 4597758a8cdSBarry Smith o local2global holds global number of column i (i=0,...,m-1) 4607758a8cdSBarry Smith o local2global holds global number of row i (i=0,...,n-1) 4617758a8cdSBarry Smith o mylocmatvec performs E_loc . x_loc where x_loc is an vector of 4627758a8cdSBarry Smith length m in 1-1 correspondence with local2global 4637758a8cdSBarry Smith (note that gs package takes care of communication). 4647758a8cdSBarry Smith (note do not zero out upper m-n entries!) 4657758a8cdSBarry Smith o mylocmatvec(void *grid_data, double *in, double *out) 4667758a8cdSBarry Smith 4677758a8cdSBarry Smith ML beliefs/usage: move this to to ML_XXT_factor routine 4687758a8cdSBarry Smith o my_ml holds address of ML struct associated w/E_loc, grid_data, grid_tag 4697758a8cdSBarry Smith o grid_tag, grid_data, my_ml used in 4707758a8cdSBarry Smith ML_Set_CSolve(my_ml, grid_tag, grid_data, ML_Do_CoarseDirect); 4717758a8cdSBarry Smith o grid_data used in 4727758a8cdSBarry Smith A_matvec(grid_data,v,u); 4737758a8cdSBarry Smith 4747758a8cdSBarry Smith Usage: 4757758a8cdSBarry Smith **************************************xxt.h***********************************/ 476*52f87cdaSBarry Smith extern PetscInt XXT_factor(xxt_ADT xxt_handle, /* prev. allocated xxt handle */ 477*52f87cdaSBarry Smith PetscInt *local2global, /* global column mapping */ 478*52f87cdaSBarry Smith PetscInt n, /* local num rows */ 479*52f87cdaSBarry Smith PetscInt m, /* local num cols */ 4807758a8cdSBarry Smith void *mylocmatvec, /* b_loc=A_local.x_loc */ 4817758a8cdSBarry Smith void *grid_data /* grid data for matvec */ 4827758a8cdSBarry Smith ); 4837758a8cdSBarry Smith 4847758a8cdSBarry Smith 4857758a8cdSBarry Smith /*************************************xxt.h************************************ 4867758a8cdSBarry Smith Function: XXT_solve 4877758a8cdSBarry Smith 4887758a8cdSBarry Smith Input : ADT ptr, b (rhs) 4897758a8cdSBarry Smith Output: x (soln) 4907758a8cdSBarry Smith Return: 4917758a8cdSBarry Smith Description: This function performs x = E^-1.b 4927758a8cdSBarry Smith Usage: 4937758a8cdSBarry Smith XXT_solve(xxt_handle, double *x, double *b) 4947758a8cdSBarry Smith XXT_solve(xxt_handle, double *x, NULL) 4957758a8cdSBarry Smith assumes x has been initialized to be b 4967758a8cdSBarry Smith **************************************xxt.h***********************************/ 497*52f87cdaSBarry Smith extern PetscInt XXT_solve(xxt_ADT xxt_handle, double *x, double *b); 4987758a8cdSBarry Smith 4997758a8cdSBarry Smith 5007758a8cdSBarry Smith /*************************************xxt.h************************************ 5017758a8cdSBarry Smith Function: XXT_stats 5027758a8cdSBarry Smith 5037758a8cdSBarry Smith Input : handle 5047758a8cdSBarry Smith Output: 5057758a8cdSBarry Smith Return: 5067758a8cdSBarry Smith Description: 5077758a8cdSBarry Smith factor stats 5087758a8cdSBarry Smith **************************************xxt.h***********************************/ 509*52f87cdaSBarry Smith extern PetscInt XXT_stats(xxt_ADT xxt_handle); 5107758a8cdSBarry Smith 5117758a8cdSBarry Smith 5127758a8cdSBarry Smith /*************************************xxt.h************************************ 5137758a8cdSBarry Smith Function: XXT_sp_1() 5147758a8cdSBarry Smith 5157758a8cdSBarry Smith Input : pointer to ADT 5167758a8cdSBarry Smith Output: 5177758a8cdSBarry Smith Return: 5187758a8cdSBarry Smith Description: sets xxt parameter 1 in xxt_handle 5197758a8cdSBarry Smith Usage: implement later 5207758a8cdSBarry Smith 5217758a8cdSBarry Smith void XXT_sp_1(xxt_handle,parameter 1 value) 5227758a8cdSBarry Smith **************************************xxt.h***********************************/ 5237758a8cdSBarry Smith 5247758a8cdSBarry Smith 5257758a8cdSBarry Smith /*************************************xyt.h************************************ 5267758a8cdSBarry Smith Module Name: xyt 5277758a8cdSBarry Smith Module Info: need xyt.{c,h} gs.{c,h} comm.{c,h} ivec.{c,h} error.{c,h} 5287758a8cdSBarry Smith 5297758a8cdSBarry Smith author: Henry M. Tufo III 5307758a8cdSBarry Smith e-mail: hmt@asci.uchicago.edu 5317758a8cdSBarry Smith contact: 5327758a8cdSBarry Smith +--------------------------------+--------------------------------+ 5337758a8cdSBarry Smith |MCS Division - Building 221 |Department of Computer Science | 5347758a8cdSBarry Smith |Argonne National Laboratory |Ryerson 152 | 5357758a8cdSBarry Smith |9700 S. Cass Avenue |The University of Chicago | 5367758a8cdSBarry Smith |Argonne, IL 60439 |Chicago, IL 60637 | 5377758a8cdSBarry Smith |(630) 252-5354/5986 ph/fx |(773) 702-6019/8487 ph/fx | 5387758a8cdSBarry Smith +--------------------------------+--------------------------------+ 5397758a8cdSBarry Smith 5407758a8cdSBarry Smith Last Modification: 3.20.01 5417758a8cdSBarry Smith **************************************xyt.h***********************************/ 5427758a8cdSBarry Smith 5437758a8cdSBarry Smith /*************************************xyt.h************************************ 5447758a8cdSBarry Smith File Description: 5457758a8cdSBarry Smith **************************************xyt.h***********************************/ 5467758a8cdSBarry Smith 5477758a8cdSBarry Smith /*************************************xyt.h************************************ 5487758a8cdSBarry Smith Notes on Usage: 5497758a8cdSBarry Smith **************************************xyt.h***********************************/ 5507758a8cdSBarry Smith 5517758a8cdSBarry Smith 5527758a8cdSBarry Smith 5537758a8cdSBarry Smith typedef struct xyt_CDT *xyt_ADT; 5547758a8cdSBarry Smith 5557758a8cdSBarry Smith 5567758a8cdSBarry Smith /*************************************xyt.h************************************ 5577758a8cdSBarry Smith Function: XYT_new() 5587758a8cdSBarry Smith 5597758a8cdSBarry Smith Input : 5607758a8cdSBarry Smith Output: 5617758a8cdSBarry Smith Return: ADT ptr or NULL upon failure. 5627758a8cdSBarry Smith Description: This function allocates and returns an xyt handle 5637758a8cdSBarry Smith Usage: xyt_handle = xyt_new(); 5647758a8cdSBarry Smith **************************************xyt.h***********************************/ 5657758a8cdSBarry Smith extern xyt_ADT XYT_new(void); 5667758a8cdSBarry Smith 5677758a8cdSBarry Smith 5687758a8cdSBarry Smith /*************************************xyt.h************************************ 5697758a8cdSBarry Smith Function: XYT_free() 5707758a8cdSBarry Smith 5717758a8cdSBarry Smith Input : pointer to ADT. 5727758a8cdSBarry Smith Output: 5737758a8cdSBarry Smith Return: 5747758a8cdSBarry Smith Description: This function frees the storage associated with an xyt handle 5757758a8cdSBarry Smith Usage: XYT_free(xyt_handle); 5767758a8cdSBarry Smith **************************************xyt.h***********************************/ 577*52f87cdaSBarry Smith EXTERN PetscInt XYT_free(xyt_ADT xyt_handle); 5787758a8cdSBarry Smith 5797758a8cdSBarry Smith 5807758a8cdSBarry Smith /*************************************xyt.h************************************ 5817758a8cdSBarry Smith Function: XYT_factor 5827758a8cdSBarry Smith 5837758a8cdSBarry Smith Input : ADT ptr, and pointer to object 5847758a8cdSBarry Smith Output: 5857758a8cdSBarry Smith Return: 0 on failure, 1 on success 5867758a8cdSBarry Smith Description: This function sets the xyt solver 5877758a8cdSBarry Smith 5887758a8cdSBarry Smith xyt assumptions: given n rows of global coarse matrix (E_loc) where 5897758a8cdSBarry Smith o global dofs N = sum_p(n), p=0,P-1 5907758a8cdSBarry Smith (i.e. row dist. with no dof replication) 5917758a8cdSBarry Smith (5.21.00 will handle dif replication case) 5927758a8cdSBarry Smith o m is the number of columns in E_loc (m>=n) 5937758a8cdSBarry Smith o local2global holds global number of column i (i=0,...,m-1) 5947758a8cdSBarry Smith o local2global holds global number of row i (i=0,...,n-1) 5957758a8cdSBarry Smith o mylocmatvec performs E_loc . x_loc where x_loc is an vector of 5967758a8cdSBarry Smith length m in 1-1 correspondence with local2global 5977758a8cdSBarry Smith (note that gs package takes care of communication). 5987758a8cdSBarry Smith (note do not zero out upper m-n entries!) 5997758a8cdSBarry Smith o mylocmatvec(void *grid_data, double *in, double *out) 6007758a8cdSBarry Smith 6017758a8cdSBarry Smith ML beliefs/usage: move this to to ML_XYT_factor routine 6027758a8cdSBarry Smith o my_ml holds address of ML struct associated w/E_loc, grid_data, grid_tag 6037758a8cdSBarry Smith o grid_tag, grid_data, my_ml used in 6047758a8cdSBarry Smith ML_Set_CSolve(my_ml, grid_tag, grid_data, ML_Do_CoarseDirect); 6057758a8cdSBarry Smith o grid_data used in 6067758a8cdSBarry Smith A_matvec(grid_data,v,u); 6077758a8cdSBarry Smith 6087758a8cdSBarry Smith Usage: 6097758a8cdSBarry Smith **************************************xyt.h***********************************/ 610*52f87cdaSBarry Smith extern PetscInt XYT_factor(xyt_ADT xyt_handle, /* prev. allocated xyt handle */ 611*52f87cdaSBarry Smith PetscInt *local2global, /* global column mapping */ 612*52f87cdaSBarry Smith PetscInt n, /* local num rows */ 613*52f87cdaSBarry Smith PetscInt m, /* local num cols */ 6147758a8cdSBarry Smith void *mylocmatvec, /* b_loc=A_local.x_loc */ 6157758a8cdSBarry Smith void *grid_data /* grid data for matvec */ 6167758a8cdSBarry Smith ); 6177758a8cdSBarry Smith 6187758a8cdSBarry Smith 6197758a8cdSBarry Smith /*************************************xyt.h************************************ 6207758a8cdSBarry Smith Function: XYT_solve 6217758a8cdSBarry Smith 6227758a8cdSBarry Smith Input : ADT ptr, b (rhs) 6237758a8cdSBarry Smith Output: x (soln) 6247758a8cdSBarry Smith Return: 6257758a8cdSBarry Smith Description: This function performs x = E^-1.b 6267758a8cdSBarry Smith Usage: XYT_solve(xyt_handle, double *x, double *b) 6277758a8cdSBarry Smith **************************************xyt.h***********************************/ 628*52f87cdaSBarry Smith extern PetscInt XYT_solve(xyt_ADT xyt_handle, double *x, double *b); 6297758a8cdSBarry Smith 6307758a8cdSBarry Smith 6317758a8cdSBarry Smith /*************************************xyt.h************************************ 6327758a8cdSBarry Smith Function: XYT_stats 6337758a8cdSBarry Smith 6347758a8cdSBarry Smith Input : handle 6357758a8cdSBarry Smith Output: 6367758a8cdSBarry Smith Return: 6377758a8cdSBarry Smith Description: 6387758a8cdSBarry Smith factor stats 6397758a8cdSBarry Smith **************************************xyt.h***********************************/ 640*52f87cdaSBarry Smith extern PetscInt XYT_stats(xyt_ADT xyt_handle); 6417758a8cdSBarry Smith 6427758a8cdSBarry Smith 6437758a8cdSBarry Smith /*************************************xyt.h************************************ 6447758a8cdSBarry Smith Function: XYT_sp_1() 6457758a8cdSBarry Smith 6467758a8cdSBarry Smith Input : pointer to ADT 6477758a8cdSBarry Smith Output: 6487758a8cdSBarry Smith Return: 6497758a8cdSBarry Smith Description: sets xyt parameter 1 in xyt_handle 6507758a8cdSBarry Smith Usage: implement later 6517758a8cdSBarry Smith 6523fdc5746SBarry Smith PetscErrorCode XYT_sp_1(xyt_handle,parameter 1 value) 6537758a8cdSBarry Smith **************************************xyt.h***********************************/ 6547758a8cdSBarry Smith 6557758a8cdSBarry Smith /********************************bit_mask.h************************************ 6567758a8cdSBarry Smith 6577758a8cdSBarry Smith Author: Henry M. Tufo III 6587758a8cdSBarry Smith 6597758a8cdSBarry Smith e-mail: hmt@cs.brown.edu 6607758a8cdSBarry Smith 6617758a8cdSBarry Smith snail-mail: 6627758a8cdSBarry Smith Division of Applied Mathematics 6637758a8cdSBarry Smith Brown University 6647758a8cdSBarry Smith Providence, RI 02912 6657758a8cdSBarry Smith 6667758a8cdSBarry Smith Last Modification: 6677758a8cdSBarry Smith 11.21.97 6687758a8cdSBarry Smith *********************************bit_mask.h***********************************/ 6697758a8cdSBarry Smith 6707758a8cdSBarry Smith /********************************bit_mask.h************************************ 6717758a8cdSBarry Smith File Description: 6727758a8cdSBarry Smith ----------------- 6737758a8cdSBarry Smith 6747758a8cdSBarry Smith *********************************bit_mask.h***********************************/ 6757758a8cdSBarry Smith 6767758a8cdSBarry Smith 6777758a8cdSBarry Smith /********************************bit_mask.h************************************ 6787758a8cdSBarry Smith Function: 6797758a8cdSBarry Smith 6807758a8cdSBarry Smith Input : 6817758a8cdSBarry Smith Output: 6827758a8cdSBarry Smith Return: 6837758a8cdSBarry Smith Description: 6847758a8cdSBarry Smith Usage: 6857758a8cdSBarry Smith *********************************bit_mask.h***********************************/ 686*52f87cdaSBarry Smith extern PetscInt div_ceil(PetscInt numin, PetscInt denom); 687*52f87cdaSBarry Smith extern PetscErrorCode set_bit_mask(PetscInt *bm, PetscInt len, PetscInt val); 688*52f87cdaSBarry Smith extern PetscInt len_bit_mask(PetscInt num_items); 689*52f87cdaSBarry Smith extern PetscInt ct_bits(char *ptr, PetscInt n); 690*52f87cdaSBarry Smith extern PetscErrorCode bm_to_proc(char *ptr, PetscInt p_mask, PetscInt *msg_list); 691*52f87cdaSBarry Smith extern PetscInt len_buf(PetscInt item_size, PetscInt num_items); 6927758a8cdSBarry Smith 6937758a8cdSBarry Smith #endif 6947758a8cdSBarry Smith 695