00001
00002
00003
00004
00005
00006
00007 #include "radiationmat.h"
00008 #include "stochdrivert.h"
00009
00010 radiationmat::radiationmat (void)
00011 {
00012 k=0.0;
00013 c=0.0;
00014 }
00015 radiationmat::~radiationmat (void)
00016 {}
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 void radiationmat::matcond (matrix &d,long ri,long ci,long ipp)
00032 {
00033 long n;
00034 n = d.n;
00035
00036 switch (n){
00037 case 1:{
00038 matcond1d (d,ri,ci,ipp);
00039 break;
00040 }
00041 case 2:{
00042 matcond2d (d,ri,ci,ipp);
00043 break;
00044 }
00045 case 3:{
00046 matcond3d (d,ri,ci,ipp);
00047 break;
00048 }
00049 default:{
00050 print_err("unknown number of components of conductivity tensor is required",__FILE__,__LINE__,__func__);
00051 }
00052 }
00053 }
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064 void radiationmat::matcond1d (matrix &d,long ri,long ci,long ipp)
00065 {
00066 double kk;
00067
00068 kk = get_k();
00069
00070 d[0][0] = kk;
00071 }
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081 void radiationmat::matcond2d (matrix &d,long ri,long ci,long ipp)
00082 {
00083 double kk;
00084
00085 kk = get_k();
00086
00087 fillm(0.0,d);
00088
00089 d[0][0] = kk; d[0][1] = 0.0;
00090 d[1][0] = 0.0; d[1][1] = kk;
00091 }
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102 void radiationmat::matcond3d (matrix &d,long ri,long ci,long ipp)
00103 {
00104 double kk;
00105
00106 kk = get_k();
00107
00108 fillm(0.0,d);
00109
00110 d[0][0]=kk; d[0][1]=0.0; d[0][2]=0.0;
00111 d[1][0]=0.0; d[1][1]=kk; d[1][2]=0.0;
00112 d[2][0]=0.0; d[2][1]=0.0; d[2][2]=kk;
00113 }
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124 void radiationmat::matcap (double &cc,long ri,long ci,long ipp)
00125 {
00126 cc = 0.0;
00127
00128 cc = get_c();
00129 }
00130
00131
00132
00133
00134
00135
00136
00137 void radiationmat::read (XFILE *in)
00138 {
00139 xfscanf (in,"%lf %lf",&c,&k);
00140 }
00141
00142
00143
00144
00145
00146
00147
00148 void radiationmat::print (FILE *out)
00149 {
00150 fprintf (out," %e %e",c,k);
00151 }
00152
00153
00154
00155
00156
00157
00158
00159 double radiationmat::get_k()
00160 {
00161 return(k);
00162 }
00163
00164
00165
00166
00167
00168
00169 double radiationmat::get_c()
00170 {
00171 return(c);
00172 }
00173
00174
00175
00176
00177
00178
00179 void radiationmat::changeparam (atsel &atm,vector &val)
00180 {
00181 long i;
00182
00183 for (i=0;i<atm.num;i++){
00184 switch (atm.atrib[i]){
00185 case 0:{
00186 c=val[0];
00187 break;
00188 }
00189 case 1:{
00190 k=val[1];
00191 break;
00192 }
00193 default:{
00194 print_err("wrong number of atribute is required",__FILE__,__LINE__,__func__);
00195 }
00196 }
00197 }
00198 }