00001 #include "plssolver.h"
00002 #include "pglobal.h"
00003 #include "seqfilesm.h"
00004 #include "genfile.h"
00005 #include "vector.h"
00006 #include <string.h>
00007 #include "mpi.h"
00008
00009 void par_solve_linear_statics ()
00010 {
00011 long i;
00012 double norf, *lhs,*rhs;
00013
00014
00015 i=0;
00016
00017 lhs = Lsrs->give_lhs (0);
00018 rhs = Lsrs->give_rhs (0);
00019
00020
00021
00022 stiffness_matrix (0);
00023
00024 Psolm->computation_stat (Gtm,Smat);
00025
00026 rhs = Lsrs->give_rhs (0);
00027
00028
00029
00030 Mb->lc[i].assemble (i,rhs+i*Ndofm,NULL,1.0);
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051 if (Psolm->tdd==layered_plate){
00052 double *th;
00053
00054
00055 th = new double [Mt->nn];
00056 for (i=0;i<Mt->nn;i++){
00057 th[i]=Mc->give_onethickness (Mt->nodes[i].crst,Mt->nodes[i].idcs);
00058 }
00059
00060
00061 Psolm->constr_mat (th,Out);
00062
00063 delete [] th;
00064 }
00065
00066
00067 fflush(stdout);
00068 fflush(stderr);
00069
00070
00071
00072
00073
00074
00075
00076
00077 Psolm->par_linear_solver (Gtm,Smat,lhs,rhs,Out,Mespr);
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113 print_init(-1, "wt",Pmp->fni,Pmp->fei);
00114 print_flush ();
00115 for (i=0;i<Lsrs->nlc;i++){
00116 print_step(i, 0, 0.0, NULL);
00117 print_flush ();
00118 }
00119 print_close();
00120
00121
00122 if (Mp->adaptivityflag)
00123 Ada->run (2,0,0);
00124
00125 Psolm->computation_stat_print (Out);
00126
00127 }
00128