00001 #ifndef DPFETI_NEW_H
00002 #define DPFETI_NEW_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 extract_from_local_vector (double *ev,double *lv);
00030 void put_into_local_vector (double *ev,double *lv);
00031 void extract_from_global_vector (double *ev,double *gv,long nsub);
00032 void put_into_global_vector (double *ev,double *gv,long nsub);
00033
00034 void arrmatrix (double *condmat);
00035 void vectors_br_bm (gtopology *top,gmatrix *gm,
00036 double *condmat,double *condvect,double *rhs);
00037
00038 void rhs_dpfeti (gtopology *top,long *domproc,gmatrix *gm);
00039
00040 void matxvect (gtopology *top,long *domproc,
00041 gmatrix *gm,double *input,double *output);
00042
00043 void cg (gtopology *top,long *domproc,gmatrix *gm,long iv,FILE *out);
00044
00045 void corner_displ (gtopology *top,long *domproc,gmatrix *gm);
00046
00047 void compute_displ (gtopology *top,gmatrix *gm,long *domproc,
00048 double *subdispl,double *rhs);
00049
00050 void solve_system (gtopology *top,gmatrix *gm,
00051 long *domproc,double *lhs,double *rhs,FILE *out,long mespr);
00052
00053
00054 int nproc;
00055
00056 int myrank;
00057
00058 int ndom;
00059
00060
00061 long ndof;
00062
00063 long nidof;
00064
00065 long ncdof;
00066
00067 long tncdof;
00068
00069 long tnmdof;
00070
00071
00072
00073 long maxncdof;
00074 long maxnbdof;
00075
00076
00077
00078
00079
00080
00081
00082
00083 long nicg;
00084
00085 long anicg;
00086
00087 double errcg;
00088
00089 double aerrcg;
00090
00091 double zero;
00092
00093 double limit,omega,indegamma;
00094
00095
00096
00097
00098
00099 double *krc;
00100
00101
00102
00103
00104
00105
00106 long **loccn;
00107
00108 long **globcn;
00109
00110 long *nlbdof;
00111
00112 long *ncdofdom;
00113
00114
00115
00116
00117 redsystsolver trssol;
00118
00119
00120 storagetype rsmstor;
00121
00122
00123 linsolvertype tlinsol;
00124
00125
00126 precondtype tprec;
00127
00128 gtopology *ptop;
00129
00130
00131 double *br;
00132
00133 double *bm;
00134
00135
00136 gmatrix *arr;
00137
00138 double *cgrhs;
00139 double *cglhs;
00140 double *displ;
00141 };
00142
00143
00144 #endif