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 #include "selnodes.h"
00010 #include "../GEFEL/slesolv.h"
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024 class dpfeti
00025 {
00026 public:
00027 dpfeti (int np,int mr,long nd);
00028 ~dpfeti();
00029 void initiate (selnodes *selnodschur,selnodes *selnodfeti);
00030 void cornodedetection (gtopology *top,long *ltg,long *domproc,FILE *out);
00031 void dpfetiordering (gtopology *top,long *ltg);
00032 void globcnnum_dpfeti (gtopology *top,long *ltg,long *domproc,FILE *out);
00033 void extract_from_local_vector (double *ev,double *lv);
00034 void put_into_local_vector (double *ev,double *lv);
00035 void extract_from_global_vector (double *ev,double *gv,long nsub);
00036 void put_into_global_vector (double *ev,double *gv,long nsub);
00037
00038 void arrmatrix (double *condmat,long *domproc,FILE *out);
00039 void vectors_br_bm (gtopology *top,gmatrix *gm,
00040 double *condmat,double *condvect,double *rhs,long *domproc,FILE *out);
00041
00042 void rhs_dpfeti (gtopology *top,long *domproc,gmatrix *gm,FILE *out);
00043
00044 void matxvect (gtopology *top,long *domproc,
00045 gmatrix *gm,double *input,double *output,FILE *out);
00046
00047 void cg (gtopology *top,long *domproc,gmatrix *gm,long iv,FILE *out);
00048
00049 void corner_displ (gtopology *top,long *domproc,gmatrix *gm,FILE *out);
00050
00051 void compute_displ (gtopology *top,gmatrix *gm,long *domproc,
00052 double *subdispl,double *rhs,FILE *out);
00053
00054 void solve_system (gtopology *top,gmatrix *gm,
00055 long *domproc,double *lhs,double *rhs,FILE *out,long mespr);
00056
00057
00058 int nproc;
00059
00060 int myrank;
00061
00062 long ndom;
00063
00064
00065 long ndof;
00066
00067 long nidof;
00068
00069 long ncdof;
00070
00071 long tncdof;
00072
00073 long tnmdof;
00074
00075
00076
00077 long maxncdof;
00078 long maxnbdof;
00079 long nlbdof;
00080
00081
00082 long *localcn;
00083 long *globalcn;
00084
00085
00086
00087
00088
00089
00090 long nicgdpfeti;
00091
00092 long anicgdpfeti;
00093
00094 double errcgdpfeti;
00095
00096 double aerrcgdpfeti;
00097
00098 double zero;
00099
00100
00101
00102 double *krc;
00103
00104
00105
00106
00107
00108
00109 long **loccn;
00110
00111 long **globcn;
00112
00113 long *ncdofdom;
00114
00115 long *nlbdofdom;
00116
00117
00118
00119 redsystsolver trssol;
00120
00121
00122 storagetype rsmstor;
00123
00124
00125 slesolv *ssle;
00126
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