00001 #ifndef DPFETI_H
00002 #define DPFETI_H
00003
00004 #include <stdio.h>
00005 #include "mpi.h"
00006 #include "../GEFEL/gtopology.h"
00007 #include "../GEFEL/gmatrix.h"
00008 #include "pgalias.h"
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 class dpfeti
00023 {
00024 public:
00025 dpfeti (int np,int mr,int nd);
00026 ~dpfeti();
00027 void dpfetiordering (gtopology *top,long *ltg);
00028 void globcnnum_dpfeti (gtopology *top,long *ltg,long *domproc,FILE *out);
00029 void locglobfeti (gtopology *top,double *gv,double *lv);
00030 void globlocfeti (gtopology *top,double *gv,double *lv);
00031
00032 void arrmatrix (double *condmat);
00033 void vectors_br_bm (gtopology *top,gmatrix *gm,
00034 double *condmat,double *condvect,double *rhs);
00035
00036 void rhs_dpfeti (gtopology *top,long *domproc,gmatrix *gm);
00037
00038 void matxvect (gtopology *top,long *domproc,
00039 gmatrix *gm,double *input,double *output);
00040
00041 void cg (gtopology *top,long *domproc,gmatrix *gm,long iv,FILE *out);
00042
00043 void corner_displ (gtopology *top,long *domproc,gmatrix *gm);
00044
00045 void compute_displ (gtopology *top,gmatrix *gm,long *domproc,
00046 double *subdispl,double *rhs);
00047
00048 void solve_system (gtopology *top,gmatrix *gm,
00049 long *domproc,double *lhs,double *rhs,FILE *out,long mespr);
00050
00051
00052 int nproc;
00053
00054 int myrank;
00055
00056 int ndom;
00057
00058
00059 long ndof;
00060
00061 long nidof;
00062
00063 long ncdof;
00064
00065 long tncdof;
00066
00067 long tnmdof;
00068
00069
00070 long ncn;
00071
00072 long nbn;
00073
00074 long maxncdof;
00075
00076 long totmaxndofn;
00077
00078 long maxinc;
00079
00080
00081 long nicg;
00082
00083 long anicg;
00084
00085 double errcg;
00086
00087 double aerrcg;
00088
00089 double zero;
00090
00091 double limit,omega,indegamma;
00092
00093
00094 long **lcngcn;
00095
00096 long *inc;
00097
00098 long *nodnum;
00099
00100
00101
00102 double *krc;
00103
00104
00105
00106
00107
00108
00109 long *nbndom;
00110
00111 long *ncndom;
00112
00113 long *ncdofdom;
00114
00115 long **masgcn;
00116
00117
00118 redsystsolver trssol;
00119
00120
00121 storagetype rsmstor;
00122
00123
00124 linsolvertype tlinsol;
00125
00126
00127 precondtype tprec;
00128
00129 gtopology *ptop;
00130
00131
00132 double *br;
00133
00134 double *bm;
00135
00136
00137 gmatrix *arr;
00138
00139 double *cgrhs;
00140 double *cglhs;
00141 double *displ;
00142 };
00143
00144
00145 #endif