00001 #include "pllssolver.h"
00002 #include "pglobal.h"
00003 #include "seqfilesm.h"
00004 #include "genfile.h"
00005 #include <string.h>
00006 #include "mpi.h"
00007
00008 void par_solve_layered_linear_statics ()
00009 {
00010 long i;
00011 double *lhs,*rhs,*th;
00012
00013
00014 i=0;
00015
00016 lhs = Lsrs->give_lhs (0);
00017 rhs = Lsrs->give_rhs (0);
00018
00019
00020 stiffness_matrix (0);
00021
00022
00023 th = new double [Mt->nn];
00024 for (i=0;i<Mt->nn;i++){
00025 th[i]=Mc->give_onethickness (Mt->nodes[i].crst,Mt->nodes[i].idcs);
00026 }
00027
00028
00029 Psolm->constr_mat (th,Out);
00030
00031 delete [] th;
00032
00033
00034 rhs = Lsrs->give_rhs (0);
00035
00036 Mb->lc[i].assemble (i,rhs+i*Ndofm,NULL,1.0);
00037
00038
00039
00040 Psolm->par_linear_solver (Gtm,Smat,lhs,rhs,Out,Mespr);
00041
00042
00043
00044 for (i=0;i<Lsrs->nlc;i++){
00045
00046 compute_req_val (i);
00047
00048
00049
00050 }
00051
00052 }
00053