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