00001 #ifndef SCHURCOMPL_H
00002 #define SCHURCOMPL_H
00003
00004 #include <stdio.h>
00005 #include "mpi.h"
00006 #include "gtopology.h"
00007 #include "gmatrix.h"
00008 #include "slesolv.h"
00009 #include "pgalias.h"
00010 #include "partop.h"
00011 #include "partopjk.h"
00012 #include "selnodes.h"
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055 class schurcompl
00056 {
00057 public:
00058 schurcompl(int np,int mr,long nd);
00059 ~schurcompl ();
00060
00061
00062
00063 void initiate (partop *ptop,selnodes *selnodschur);
00064 void initiate (partopjk *ptop,FILE *out);
00065
00066 void schurordering_new (gtopology *top, long *nodeidentif);
00067 void globcnnum_pdd (gtopology *top,long *ltg,long *domproc,FILE *out);
00068
00069 void solve_red_sys_iter (long *domproc,double *condmat,double *condvect,FILE *out);
00070 void solve_red_sys_fin (long *domproc,double *condmat,double *condvect,long decomp,FILE *out);
00071 void solve_red_sys (long *domproc,double *condmat,double *condvect,long decomp,FILE *out);
00072
00073 void solve_system (gtopology *top,gmatrix *gm,long *domproc,double *lhs,double *rhs,FILE *out);
00074
00075 void gather_bound_vect (double *lv,double *gv,long *domproc);
00076 double pss_gather_bound_vect (double *lv1,double *gv1,double *lv2,double *gv2,long *domproc);
00077 long selected_norm_calculation (long cid,double err,double thresh,partopjk *ptop,gtopology *top,double *lv1,double *gv1,double *rhs,double *grhs,long *domproc,long n,FILE *out);
00078 double pss_gather_bound_vect_old (double *lv,double *gv,long *domproc);
00079 void scatter_bound_vect (double *lv,double *gv,long *domproc);
00080
00081 double unbalanced_values (double *lv,long *domproc,FILE *out);
00082
00083
00084 int nproc;
00085
00086 int myrank;
00087
00088 long ndom;
00089
00090
00091 long ndof;
00092
00093 long nidof;
00094
00095 long nbdof;
00096
00097
00098
00099
00100
00101
00102
00103
00104 long maxnbdof;
00105
00106
00107 long nicgsch;
00108
00109 long anicgsch;
00110
00111 double errcgsch;
00112
00113 double aerrcgsch;
00114
00115 double zero;
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126 long ndofcp;
00127
00128
00129 redsystsolver trssol;
00130
00131
00132 storagetype rsmstor;
00133
00134
00135 slesolv *ssle;
00136
00137
00138 storagetype dmstor;
00139
00140
00141
00142
00143
00144 long *nbdofmas;
00145
00146 long destr_nbdofmas;
00147
00148 gtopology *gtop;
00149 gmatrix *arr;
00150
00151 };
00152
00153
00154 #endif