[method] pde = 1 qoi = 1 ;qoi = 4 [qoi_0] init_time_window = [0,122] noffset = 0 numideal = 0 plotoptsolve = true [perfusion] w_0_probe = 48.35e0 ;w_0_healthy = 6.0;12.0;18.0;24.0;30.0;42.0;54.0;66.0;72.0;84.0 ;w_0_healthy = 6.0;18.0;30.0;42.0 ;w_0_healthy = 6.0 w_0_healthy = 0.0;3.0;6.0;12.0 w_n = 0.0e0 w_i = 0.0e0 w_d = 0.0e0 w_2 = 0.0e0 w_ni = 305.0e0 w_id = 315.0e0 [thermal_conductivity] k_0_probe = .5e0 k_0_healthy = .527e0 k_1 = 0.0e0 k_2 = 0.0e0 k_3 = 0.0e0 [compexec] ;meshdata=s129226-files/s129226_smooth.e ;meshdata = s129226-files/s129226_cooledTip.e ;meshdata = generic-mesh/meshTemplateLoRes.e ;meshdata = /home/dfuentes/DDDAS/data/mdacc/canine_feb07/generic-mesh/meshTemplate.e meshdata = /work/00131/fuentes/data/mdacc/canine_feb07/generic-mesh/meshTemplateLowerRes.e;/work/00131/fuentes/data/mdacc/canine_feb07/generic-mesh/meshTemplateLoRes.e;/work/00131/fuentes/data/mdacc/canine_feb07/generic-mesh/meshTemplateNormRes.e;/work/00131/fuentes/data/mdacc/canine_feb07/generic-mesh/meshTemplateHiRes.e powerdata = [14,16,27,39,41,46,52,55,59,74,82,86,90,95,122]@[0.0,6.0,0.0,10.05,0.0,11.25,0.0,11.25,0.0,11.25,0.0,13.2,0.0,12.0,0.0] ;executable = image executable = dddas ;meshrefine = 1 jobid = dog2fem980 ;numproc = 1 solver = superlu_dist numproc = 16;16;32;48 runtime_options = -snes_ls basic -info -info_exclude null,vec,mat,pc,ksp,snes -ksp_rtol 1.e-9 -pc_type bjacobi -snes_monitor -snes_rtol 1.e-6 -snes_converged_reason -tao_max_funcs 0 b0 = .2335 b1 = .234 b2 = -.0335 ; FIRST R1 rotation Then R3 rotation.... (not commutative) ;-90deg R1 ---> 6deg R3 rotation (NEGATIVE OF PARAVIEW) ; | cos \theta sin \theta 0 | | 1 0 0 | ; | -sin \theta cos \theta 0 | | 0 cos \gamma sin \gamma | ; | 0 0 1 | | 0 -sin \gamma cos \gamma | ; ; | cos \theta sin \theta 0 | | 1 0 0 | ; | -sin \theta cos \theta 0 | | 0 0 -1 | ; | 0 0 1 | | 0 1 0 | ; ; | cos \theta 0 -sin \theta | ; | -sin \theta 0 -cos \theta | ; | 0 1 0 | ; a11 = 0.0 a22 = 0.0 ;>>> cos(pi/180.*6.) ; 0.99452189536827329 ;>>> sin(pi/180.*6.) ; 0.10452846326765347 a00 = 0.99452189536827329 a02 = -0.10452846326765347 a10 = -0.10452846326765347 a12 = -0.99452189536827329 a21 = 1.0 [bc] u_dirichletid = 1 u_flux = 0.0 newton_coeff = 0.0 [optical] mu_a_healthy = 5.0e2 mu_a_tumor = 2.38e2 mu_a_probe = 0.0195e0 mu_s_healthy = 140.0e2 mu_s_tumor = 3.5e2 mu_s_probe = 100.0e0 anfact = 0.88e0 [initial_condition] z_init = 7.431e4 ;probe_init = 21.0 probe_init = 21.0;35.9 u_init = 35.9 [output] ;writeinterval = 1 [probe] ;domain = 2 domain = 2;10 ;center x_0 = .2335 y_0 = .234 z_0 = -.0335 [timestep] isteps_per_ideal = 6 ideal_dt = 6.0e0 [mrti] filein = /work/00131/fuentes/data/mdacc/canine_feb07/980brainStat/dog2_980/dog2_98000/mrivis/temperature.%04d.vtk uncertaintyin = /work/00131/fuentes/data/mdacc/canine_feb07/980brainStat/dog2_980/dog2_98000/mrivis/tmapstd.%04d.vtk ;filein = /home/dfuentes/DDDAS/data/mdacc/canine_feb07/980brainStat/dog2_980/dog2_98000/mrivis/temperature.%04d.vtk ;uncertaintyin = /home/dfuentes/DDDAS/data/mdacc/canine_feb07/980brainStat/dog2_980/dog2_98000/mrivis/tmapstd.%04d.vtk exampath = /FUS4/data2/BioTex/Brain_PhaseII/061113_BrainDog2_Treat/e351850 dirid = 353172 ntime = 122 dictionary = $ITK_SOURCE/Utilities/gdcm/Dicts/GEMS.dic ixlo = 144 ixhi = 164 iylo = 113 iyhi = 137 [kalman] method = 2 magix = 65 magiy = 65 roi = [0,255,0,255]