00001 #ifndef SEQSCHUR_H
00002 #define SEQSCHUR_H
00003
00004 #include "skyline.h"
00005 #include "seqselnodes.h"
00006 #include "scr.h"
00007 #include "densemat.h"
00008
00009 class gtopology;
00010 class gmatrix;
00011 class slesolv;
00012
00013
00014
00015
00016
00017
00018
00019
00020 class seqschur
00021 {
00022 public:
00023 seqschur ();
00024 ~seqschur ();
00025
00026 void read (gtopology *top,long mespr,XFILE *in);
00027 void print (FILE *out);
00028
00029
00030
00031
00032 void initiate (seqselnodes *selnodschur);
00033
00034
00035
00036 void schurordering_new (gtopology *top, long *nodeidentif);
00037 void globcnnum_pdd (gtopology *top,long *ltg,long *domproc,FILE *out);
00038
00039 void assemble_subdom_unknowns (gtopology *top,FILE *out);
00040 void subdomain_matrices (gmatrix *gm,FILE *out);
00041
00042 void solve_red_sys_iter (double **condmat,double **condvect,FILE *out);
00043 void solve_red_sys_fin (double **condmat,double **condvect,FILE *out);
00044 void solve_red_sys (double **condmat,double **condvect,FILE *out);
00045
00046 void solve_system (gtopology *top,gmatrix *gm,double *lhs,double *rhs,FILE *out);
00047
00048
00049
00050
00051
00052
00053
00054
00055 long ns;
00056
00057
00058 redsystsolver trssol;
00059
00060
00061 slesolv *ssle;
00062
00063
00064 storagetype rsmstor;
00065
00066
00067 long nicg;
00068
00069 long anicg;
00070
00071 double errcg;
00072
00073 double aerrcg;
00074
00075 double zero;
00076
00077
00078
00079 long ndof;
00080
00081 long nidof;
00082
00083 long nbdof;
00084
00085
00086 long maxnbdof;
00087
00088
00089
00090 long ndofcp;
00091
00092
00093 long *nbdofmas;
00094
00095
00096 gtopology *gtop;
00097
00098
00099 skyline *smsky;
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114 long *nsid;
00115
00116
00117
00118
00119 long *ndofdom;
00120
00121
00122
00123 long **cndom;
00124
00125
00126
00127
00128
00129
00130
00131 storagetype dmstor;
00132
00133
00134
00135
00136
00137
00138
00139 gmatrix *arr;
00140
00141 };
00142
00143 #endif