Changeset 1126 for branches/sandbox


Ignore:
Timestamp:
Mar 29, 2011 4:07:16 PM (10 years ago)
Author:
toby
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/sandbox/rigid.tcl

    r1125 r1126  
    1616    source c:/gsas/sandboxexpgui/gsascmds.tcl
    1717    source C:/gsas/sandboxexpgui/rb.tcl
    18     puts beforeread
     18#    puts beforeread
    1919    expload c:/crystals/expgui/rigid/rb6norb.exp
    2020    mapexp
    21     puts after
     21#    puts after
    2222} else {
    2323    source [file join $expgui(scriptdir) rb.tcl]
     
    263263     set temp_car_group ""
    264264     set total ""
    265      puts $::::rb_coord_num($bodytyp,1)
     265#     puts $::::rb_coord_num($bodytyp,1)
    266266     for {set matrixnum 1} {$matrixnum <= $::rb_matrix_num($bodytyp)} {incr matrixnum} {
    267267         lappend temp_mat $::rb_mult($bodytyp,$matrixnum)
     
    275275         lappend temp_car $temp
    276276     }
    277      puts "sites: $::rb_coord_num($bodytyp,1)"
    278      puts "matrix multiplier:  $temp_mat"
    279      puts "cartesian coords:   $temp_car"
     277#     puts "sites: $::rb_coord_num($bodytyp,1)"
     278#     puts "matrix multiplier:  $temp_mat"
     279#     puts "cartesian coords:   $temp_car"
    280280     AddRigidBody $temp_mat $temp_car
    281281
     
    314314     }
    315315#    delete the rigid body
    316      puts "delete rigid body number $bodytyp"
     316#     puts "delete rigid body number $bodytyp"
    317317     DeleteRigidBody $bodytyp
    318      puts "destroy location $location"
     318#     puts "destroy location $location"
    319319     destroy $location
    320320#    increment expgui
     
    582582    set useflags {}
    583583    foreach i $coords {lappend useflags 1}
    584     puts "frcoords $frcoords"
    585     puts "coords $coords"
     584#    puts "frcoords $frcoords"
     585#    puts "coords $coords"
    586586    # do the fit
    587587    foreach {neworigin newEuler rmsdev newfrac rmsbyatom} \
     
    593593    # show deviations
    594594    foreach atom $atoms rms $rmsbyatom {
    595         puts "[atominfo $phase $atom label]\t$rms"
     595#       puts "[atominfo $phase $atom label]\t$rms"
    596596    }
    597597    #puts "CalcBody $Euler $cell $coords $origin"
     
    654654   set origin "$::origin1 $::origin2 $::origin3"
    655655   set euler "$::euler1 $::euler2 $::euler3"
    656    puts "phase = $::rb_phase"
    657    puts "bodytyp = $::body_type"
    658    puts "firstatom = $::rb_finput"
    659    puts "position = $origin"
    660    puts "Euler = $euler"
     656#   puts "phase = $::rb_phase"
     657#   puts "bodytyp = $::body_type"
     658#   puts "firstatom = $::rb_finput"
     659#   puts "position = $origin"
     660#   puts "Euler = $euler"
    661661   MapRigidBody $::rb_phase $::body_type $::rb_finput $origin $euler
    662662   incr ::rb_map($::body_type)
     
    775775
    776776proc RB_unmap_delete {panel x args} {
    777      puts $panel
     777#     puts $panel
    778778     foreach p $::expmap(phaselist) {
    779779        foreach z [RigidBodyMappingList $p $x] {
     
    858858         lappend temp_car $temp
    859859     }
    860      puts "Matrix Update Info = $bodynum $temp_mat $temp_car"
     860#     puts "Matrix Update Info = $bodynum $temp_mat $temp_car"
    861861#     ReplaceRigidBody $bodynum $temp_mat $temp_car
    862862#     incr ::expgui(changed)
     
    959959      destroy $location.display
    960960      set filelist [RB_Import_Data_Type]
    961       puts $filelist
     961#      puts $filelist
    962962#     menubutton $location.but -text "File Type" -menu $location.but.menu
    963963#      grid [frame $location.display -bd 2 -relief groove] -row 1 -column 0
     
    10691069     set start_loc [lsearch $::expmap(atomlist_$phase) $::gcon_start]
    10701070     set ::rb_atom_range [lrange $::expmap(atomlist_$phase) $start_loc [expr $start_loc + $::gcon_atoms - 1]]
    1071      puts "location = $start_loc  range = $::rb_atom_range"
     1071#     puts "location = $start_loc  range = $::rb_atom_range"
    10721072     set rownum 1
    10731073     set colnum 1
     
    10881088             if {[expr $colnum % 4] == 0} {incr rownum; set colnum 1}
    10891089             set atomid [atominfo $phase $atom  label]
    1090              puts $atomid
     1090#             puts $atomid
    10911091             set ::rb_atom_origin_set($atom) 1
    10921092             grid [checkbutton $main.$atom -text "$atomid" -variable ::rb_atom_origin_set($atom)] -row $rownum -column $colnum
     
    11131113  }
    11141114
    1115                   puts $atom_info_list
     1115#                  puts $atom_info_list
    11161116       set ::rb_param_x1 [lindex $atom_list 0]
    11171117       set ::rb_param_x2 [lindex $atom_list 1]
     
    11771177                }
    11781178        }
    1179         puts "Origin list = $::rb_origin_list"
     1179#        puts "Origin list = $::rb_origin_list"
    11801180}
    11811181
     
    12001200                }
    12011201        }
    1202         puts "Origin list = $::gcon_origin_list"
    1203         puts "vector 1 list = $vector1list"
    1204         puts "vector 2 list = $vector2list"
    1205         puts "number atoms = $::gcon_atoms"
    1206         puts "start atom  = $::gcon_start"
     1202#        puts "Origin list = $::gcon_origin_list"
     1203#        puts "vector 1 list = $vector1list"
     1204#        puts "vector 2 list = $vector2list"
     1205#        puts "number atoms = $::gcon_atoms"
     1206#        puts "start atom  = $::gcon_start"
    12071207
    12081208set temp1 [ExtractRigidBody $::rb_phase $::gcon_atoms $::gcon_start $::gcon_origin_list $vector1list $vector2list]
    12091209if {[lindex $temp1 0] == {} || [lindex $temp1 1] == {} || [lindex $temp1 2] == {}} {
    1210    puts "Geometry Crashed"
     1210#   puts "Geometry Crashed"
    12111211   }
    12121212#puts "string 1 = [lindex $temp1 0]"
     
    12161216set cartesian ""
    12171217lappend cartesian [lindex $temp1 2]
    1218 puts "Cartesian = $cartesian"
     1218#puts "Cartesian = $cartesian"
    12191219
    12201220set bodytyp [AddRigidBody 1 $cartesian]
     
    13651365                         lappend ::rb_var_list ::rb_var($bodnum,$mapnum,saa) ::rb_var($bodnum,$mapnum,sbb)
    13661366
    1367                          puts $main
     1367#                         puts $main
    13681368                         grid [button $main.cfefx($bodnum,$mapnum) -command "RB_Con_Button $main.cfefx($bodnum,$mapnum)" -textvariable ::rb_var($bodnum,$mapnum,x) -width 5] -row $row -column 4
    13691369                         grid [button $main.cfefy($bodnum,$mapnum) -command "RB_Con_Button $main.cfefy($bodnum,$mapnum)" -textvariable ::rb_var($bodnum,$mapnum,y) -width 5]  -row $row -column 5
     
    14281428             grid [button $con.con.var -width 20 -text "Set Constrained Variables" -command {RB_Con_But_Proc $::rbaddresses $::rb_var_name}] -row 4 -column 1
    14291429             grid [entry $con.con.vare -textvariable ::rb_var_name -width 5] -row 4 -column 2
    1430 
     1430             grid [button $con.con.save -width 20 -text "Assign Variables and Save" -command RB_Var_Assign] -row 5 -column 1
    14311431
    14321432
     
    14341434
    14351435proc RB_Con_But_Proc {addresses change args} {
    1436      puts "$addresses $change"
     1436#     puts "$addresses $change"
    14371437     foreach address $addresses {
    14381438             set a [eval $address cget -relief]
     
    14581458
    14591459
    1460 proc RB_Var_Assignqw {args} {
    1461      #Determine number of rigid bodies and rigid body mappings
    1462      set ::rb_var_used [RigidBodyGetNavNums]
    1463      set var_count 1
    1464      set rb_num [RigidBodyList]
    1465      set varnames ""
    1466      foreach phasenum $::expmap(phaselist) {
    1467           foreach bodnum $rb_num {
    1468                set rb_map_num($phasenum,$bodnum) [RigidBodyMappingList $phasenum $bodnum]
    1469                for {set mapnum 1} {$mapnum <= $rb_map_num($phasenum,$bodnum)} {incr mapnum} {
    1470 
    1471                }
    1472           }
    1473      }
    1474 }
    14751460
    14761461#procedure to determine next available variable number for GSAS
     
    14831468#procedure to assign variable names to relationships
    14841469proc RB_Var_Assign {args} {
     1470     set ::rb_varlist [RigidBodyGetVarNums]
    14851471     set varcount 1
    14861472     set varlist ""
     
    14911477            } elseif {[set $var] == "free"} {
    14921478                     set $var [RB_Var_Gen $varcount]
    1493                      set $varcount $var
     1479                     set varcount [set $var]
    14941480            } else {
    1495                      lappend varlist $var
    1496             }
    1497 
    1498 
    1499 
    1500      puts "$var = [set $var]"
    1501      puts "list = $varlist"
    1502      }
    1503 }
     1481#                     puts [lsearch $varlist [set $var]]
     1482                     if {[lsearch $varlist [set $var]] == -1} {
     1483                        lappend varlist [set $var]
     1484#                        puts $varlist
     1485                        set rb_variable([set $var]) [RB_Var_Gen $varcount]
     1486                        set $var $rb_variable([set $var])
     1487
     1488                     } else {
     1489                        set $var $rb_variable([set $var])
     1490                     }
     1491             }
     1492     }
     1493     RB_Var_Set
     1494}
     1495
     1496
     1497
     1498#procedure to send variable numbers to EXP file.
     1499proc RB_Var_Set {args} {
     1500     #Determine number of rigid bodies and rigid body mappings
     1501     set rb_num [RigidBodyList]
     1502     foreach phasenum $::expmap(phaselist) {
     1503          foreach bodnum $rb_num {
     1504               set rb_map_num($phasenum,$bodnum) [RigidBodyMappingList $phasenum $bodnum]
     1505                   for {set mapnum 1} {$mapnum <= $rb_map_num($phasenum,$bodnum)} {incr mapnum} {
     1506                       set rb_list "$::rb_var($bodnum,$mapnum,x) \
     1507                           $::rb_var($bodnum,$mapnum,y) $::rb_var($bodnum,$mapnum,z) \
     1508                           $::rb_var($bodnum,$mapnum,e1) $::rb_var($bodnum,$mapnum,e2) \
     1509                           $::rb_var($bodnum,$mapnum,e3) 0 0 0"
     1510                       RigidBodyVary $phasenum $bodnum $mapnum $rb_list
     1511                     
     1512                       set rb_tls "$::rb_var($bodnum,$mapnum,t11) $::rb_var($bodnum,$mapnum,t22) \
     1513                           $::rb_var($bodnum,$mapnum,t33) $::rb_var($bodnum,$mapnum,t12) \
     1514                           $::rb_var($bodnum,$mapnum,t13) $::rb_var($bodnum,$mapnum,t23) \
     1515                           $::rb_var($bodnum,$mapnum,l11) $::rb_var($bodnum,$mapnum,l22) \
     1516                           $::rb_var($bodnum,$mapnum,l33) $::rb_var($bodnum,$mapnum,l12) \
     1517                           $::rb_var($bodnum,$mapnum,l13) $::rb_var($bodnum,$mapnum,l23) \
     1518                           $::rb_var($bodnum,$mapnum,s12) $::rb_var($bodnum,$mapnum,s13) \
     1519                           $::rb_var($bodnum,$mapnum,s21) $::rb_var($bodnum,$mapnum,s23) \
     1520                           $::rb_var($bodnum,$mapnum,s31) $::rb_var($bodnum,$mapnum,s32) \
     1521                           $::rb_var($bodnum,$mapnum,saa) $::rb_var($bodnum,$mapnum,sbb)"
     1522                       RigidBodyTLSVary $phasenum $bodnum $mapnum $rb_tls
     1523
     1524                       puts $rb_list
     1525                       puts $rb_tls
     1526
     1527               }
     1528          }
     1529     }
     1530}
     1531
     1532
     1533#RB_Refine_Con - sets matrix to set refinement constrols
     1534#RB_Var_Assign - assigns variable numbers to refinement parameters
     1535#RB_Var_Set    - sends variable numbers back to EXP file.
Note: See TracChangeset for help on using the changeset viewer.