00001 #include "spsolvert.h"
00002 #include "globalt.h"
00003 #include "globmatt.h"
00004 #include "transprint.h"
00005 #include <string.h>
00006 #include "npsolvert.h"
00007
00008
00009
00010
00011
00012
00013 void solve_stationary_problem ()
00014 {
00015 long i;
00016 double *rhs;
00017
00018
00019 conductivity_matrix (0);
00020
00021 rhs=Lsrst->give_rhs (0);
00022
00023
00024
00025 i = 0;
00026 trfel_right_hand_side (i,rhs,Ndoft);
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
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069 i = 0;
00070 Tp->ssle->solve_system (Gtt,Kmat,Lsrst->give_lhs(i),Lsrst->give_rhs(i),Outt);
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091 compute_req_valt (0);
00092
00093 print_initt(-1, "wt");
00094
00095
00096 print_stept(i, 0, 0.0, rhs);
00097
00098 print_closet();
00099
00100
00101 if (Tp->adaptivityflag)
00102 Adat->run (2, false);
00103 }
00104
00105
00106
00107
00108
00109
00110
00111 void solve_radiation_stationary_problem ()
00112 {
00113 long i,ni,lcid;
00114 double *rhs;
00115
00116
00117 conductivity_matrix (0);
00118
00119 rhs=Lsrst->give_rhs (0);
00120
00121 lcid = 0;
00122
00123 ni=100;
00124
00125 print_initt(-1, "wt");
00126 for (i=0;i<ni;i++){
00127
00128 nullv (rhs,Ndoft);
00129 trfel_right_hand_side (lcid,rhs,Ndoft);
00130
00131 Tt->edge_temperature ();
00132 Tt->heat_fluxes (rhs,Outt);
00133
00134 fprintf (Outt,"\n\n\n i %ld",i);
00135
00136
00137
00138
00139 fprintf (Outt,"\n\n\n matice vodivosti pred eliminaci, krok %ld",i);
00140 Kmat->printmat (Outt);
00141 Tp->ssle->solve_system (Gtt,Kmat,Lsrst->give_lhs(0),Lsrst->give_rhs(0),Outt);
00142 fprintf (Outt,"\n\n\n matice vodivosti po eliminaci, krok %ld",i);
00143 Kmat->printmat (Outt);
00144
00145 print_stept(lcid, 0, 0.0, rhs);
00146
00147 }
00148
00149
00150
00151
00152
00153 print_closet();
00154 }
00155