00001 #include "pspsolvert.h"
00002 #include "pglobalt.h"
00003 #include "seqfilest.h"
00004 #include <string.h>
00005 #include <math.h>
00006 #include "mpi.h"
00007
00008 void par_solve_stationary_problem ()
00009 {
00010 long i;
00011 double *lhs,*rhs;
00012
00013
00014 i=0;
00015
00016
00017 conductivity_matrix (0);
00018
00019 rhs = Lsrst->give_rhs (0);
00020 lhs = Lsrst->give_lhs (0);
00021
00022
00023
00024 trfel_right_hand_side (0,rhs,Ndoft);
00025
00026
00027
00028 Psolt->par_linear_solver (Gtt,Kmat,lhs,rhs,Outt,Mesprt);
00029
00030 print_initt(-1, "wt");
00031 for (i=0;i<Lsrst->nlc;i++){
00032
00033 print_stept(i, 0, 0.0, NULL);
00034 }
00035 print_closet();
00036
00037 }
00038