No ASEG option present
 No field output option specified
	 will use default value out_field.mgz
 will cache transform in file /home/decharms/data/freesurfer/round4/9033b/cvs/transform.txt
 reading fixed volume
 reading moving volume
 surface name = /home/decharms/data/freesurfer/round4/9032b/surf/lh.white type = 0
 surface name = /home/decharms/data/freesurfer/round4/9032b/surf/rh.white type = 0
 surface name = /home/decharms/data/freesurfer/round4/9032b/surf/lh.pial type = 1
 surface name = /home/decharms/data/freesurfer/round4/9032b/surf/rh.pial type = 1
 processing APARC
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
 label = unknown -> value = 0
 label = corpuscallosum -> value = 4
 ignored labels = 0 4 
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
 label = unknown -> value = 0
 label = corpuscallosum -> value = 4
 ignored labels = 0 4 
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
 label = unknown -> value = 0
 label = corpuscallosum -> value = 4
 ignored labels = 0 4 
reading colortable from annotation file...
colortable with 36 entries read (originally /autofs/space/terrier_001/users/nicks/freesurfer/average/colortable_desikan_killiany.txt)
 label = unknown -> value = 0
 label = corpuscallosum -> value = 4
 ignored labels = 0 4 
 surface name = /home/decharms/data/freesurfer/round4/9033b/cvs/lh.resample.white type = 0
 surface name = /home/decharms/data/freesurfer/round4/9033b/cvs/rh.resample.white type = 0
 surface name = /home/decharms/data/freesurfer/round4/9033b/cvs/lh.resample.pial type = 1
 surface name = /home/decharms/data/freesurfer/round4/9033b/cvs/rh.resample.pial type = 1
 trying to read transform from /home/decharms/data/freesurfer/round4/9033b/cvs/transform.txt
 done with Powell
 applying linear transform to MOVING surf pos.
 linear transform = 
0.924437	-0.0130839	0.0363887	
0.008065	1.00276	-0.0349886	
-0.0147243	0.0784853	0.979084	
10.9923	-6.04308	12.0756	
 linear transform = 
0.924437	-0.0130839	0.0363887	
0.008065	1.00276	-0.0349886	
-0.0147243	0.0784853	0.979084	
10.9923	-6.04308	12.0756	
 linear transform = 
0.924437	-0.0130839	0.0363887	
0.008065	1.00276	-0.0349886	
-0.0147243	0.0784853	0.979084	
10.9923	-6.04308	12.0756	
 linear transform = 
0.924437	-0.0130839	0.0363887	
0.008065	1.00276	-0.0349886	
-0.0147243	0.0784853	0.979084	
10.9923	-6.04308	12.0756	
	 DONE
 direct min = 63.4216 53.5712 46.9332 direct max = 188.943 168.282 213.078
 bounding box elts = 
 min = 53.4216 43.5712 36.9332
 max = 198.943 178.282 223.078
 DIRTY value = 1
 create BC container - ripCounter = 0
 create BC container - ripCounter = 0
 create BC container - ripCounter = 0
 create BC container - ripCounter = 0
 creating container for the mesher
 writing surface to file lh.dbg_surf.white
 writing surface to file rh.dbg_surf.white
 writing surface to file lh.dbg_surf.pial
 writing surface to file rh.dbg_surf.pial
 srcPoints size = 0
 ======================
 step = 17
===============
elt_vol= 10
Constructing Delaunay tetrahedralization.
Delaunay seconds:  0
Creating surface mesh.
Perturbing vertices.
Delaunizing segments.
Constraining facets.
Segment and facet seconds:  0.01
Removing unwanted tetrahedra.
Hole seconds:  0
Repairing mesh.
Repair seconds:  0
Adding Steiner points to enforce quality.
Quality seconds:  14.58
Smoothing mesh.
Smooth seconds:  0.98

Writing nodes.
Writing elements.
Writing faces.

Output seconds:  0.09
Total running seconds:  15.66

Statistics:

  Input points: 8
  Input facets: 6
  Input holes: 0
  Input regions: 0

  Mesh points: 105023
  Mesh tetrahedra: 634922
  Mesh faces: 1284442
  Mesh subfaces: 29196
  Mesh subsegments: 613

Saving nodes to iteration.node
Saving elements to iteration.ele
Saving faces to iteration.face
 mesh nodes = 105023 mesh elts = 634922
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 InitializeOctantDeterminant
	 count inserted = 100000 elapsed = 0.214 seconds  element count = 885530
	 count inserted = 200000 elapsed = 0.374 seconds  element count = 2697247
	 count inserted = 300000 elapsed = 0.336 seconds  element count = 4401602
	 count inserted = 400000 elapsed = 0.285 seconds  element count = 5980932
	 count inserted = 500000 elapsed = 0.275 seconds  element count = 7533863
	 count inserted = 600000 elapsed = 0.266 seconds  element count = 9086211
 done building octree - total elements = 9632111
 penalty_weight = 1000
setting bc 198.943 178.282 223.078 -> 0 0 0
	 instead 198.933 178.272 223.068 -> norm = 0.0173205
setting bc 198.943 178.282 36.9332 -> 0 0 0
	 instead 198.933 178.272 36.9432 -> norm = 0.0173205
setting bc 198.943 43.5712 223.078 -> 0 0 0
	 instead 198.933 43.5812 223.068 -> norm = 0.0173205
setting bc 198.943 43.5712 36.9332 -> 0 0 0
	 instead 198.933 43.5812 36.9432 -> norm = 0.0173205
setting bc 53.4216 178.282 223.078 -> 0 0 0
	 instead 53.4316 178.272 223.068 -> norm = 0.0173205
setting bc 53.4216 178.282 36.9332 -> 0 0 0
	 instead 53.4316 178.272 36.9432 -> norm = 0.0173205
setting bc 53.4216 43.5712 223.078 -> 0 0 0
	 instead 53.4316 43.5812 223.068 -> norm = 0.0173205
setting bc 53.4216 43.5712 36.9332 -> 0 0 0
	 instead 53.4316 43.5812 36.9432 -> norm = 0.0173205
 computing statistics for the displacement application error
 average norm of error in placement = 0.0173205
 iterating
 done with candidates
 Active BCs = 121375
 Total BCs = 482528
 no-eqs = 315069
	 percentage done= 0
	 percentage done= 5
	 percentage done= 10
	 percentage done= 15
	 percentage done= 20
	 percentage done= 25
	 percentage done= 30
	 percentage done= 35
	 percentage done= 40
	 percentage done= 45
	 percentage done= 50
	 percentage done= 55
	 percentage done= 60
	 percentage done= 65
	 percentage done= 70
	 percentage done= 75
	 percentage done= 80
	 percentage done= 85
	 percentage done= 90
	 percentage done= 95
 done_bc_natural size of container = 482528
 after building the map
 LOAD size = 24
 after setup_load_mfc
Linear solve converged due to CONVERGED_RTOL iterations 750
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=315069, cols=315069
    total: nonzeros=14526249, allocated nonzeros=14526582
      using I-node routines: found 105039 nodes, limit used is 5
 convergence reason = 2
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=315069, cols=315069
    total: nonzeros=14526249, allocated nonzeros=14526582
      using I-node routines: found 105039 nodes, limit used is 5
 Iterations 750
Absolute-Norm of error = 3.0708e-05
 countInvalid = 4 general-count = 482524
 Average of the error norm = 0.0818676
 Initial error = 0.405004
 Average of the error norm conditional = 0.0679995
 Active Residual error = 0.0171124
 Total error = 0.081894
 c0324 time elapsed updating the sources (seconds) = 0
 orientation_pbs
 solving topology problems
 using old topology solver
 using Lame constants given by 
 E = 10
 nu = 0.3
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 177659
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 516
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 171522
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 540
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 174285
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 519
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 139790
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 450
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 152913
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 504
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 189544
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 594
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 146061
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 468
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 191928
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 576
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 before 
 after
 !!!! No residual problems.
	 DONE
 compressing morph
dcds mesh direct image = 53.4216 43.5712 36.9332
 dir img for computed transform = 4.24431e-313 2.36145e-306 8.67399e-307
 dir img for compressed = 4.24431e-313 2.36145e-306 8.67399e-307
dcds total image for source point = 4.24431e-313 2.36145e-306 8.67399e-307
 name = compressed_step_17
   Initial transform 
	 name = 
 FEM error
 Initial displacement = 6.88518
 Final displacement = 6.50568
 Max after displacement = 21.4363
 Error points = 0
 ======================
 step = 16
===============
elt_vol= 9.6875
Constructing Delaunay tetrahedralization.
Delaunay seconds:  0
Creating surface mesh.
Perturbing vertices.
Delaunizing segments.
Constraining facets.
Segment and facet seconds:  0
Removing unwanted tetrahedra.
Hole seconds:  0
Repairing mesh.
Repair seconds:  0
Adding Steiner points to enforce quality.
Quality seconds:  15.11
Smoothing mesh.
Smooth seconds:  1.03

Writing nodes.
Writing elements.
Writing faces.

Output seconds:  0.08
Total running seconds:  16.22

Statistics:

  Input points: 8
  Input facets: 6
  Input holes: 0
  Input regions: 0

  Mesh points: 108208
  Mesh tetrahedra: 654680
  Mesh faces: 1324277
  Mesh subfaces: 29834
  Mesh subsegments: 623

Saving nodes to iteration.node
Saving elements to iteration.ele
Saving faces to iteration.face
 mesh nodes = 108208 mesh elts = 654680
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
	 count inserted = 100000 elapsed = 0.214 seconds  element count = 864703
	 count inserted = 200000 elapsed = 0.368 seconds  element count = 2636900
	 count inserted = 300000 elapsed = 0.337 seconds  element count = 4313927
	 count inserted = 400000 elapsed = 0.284 seconds  element count = 5867587
	 count inserted = 500000 elapsed = 0.27 seconds  element count = 7392855
	 count inserted = 600000 elapsed = 0.269 seconds  element count = 8918326
 done building octree - total elements = 9758637
 penalty_weight = 1000
setting bc 198.943 178.282 223.078 -> 0 0 0
	 instead 198.933 178.272 223.068 -> norm = 0.0173205
setting bc 198.943 178.282 36.9332 -> 0 0 0
	 instead 198.933 178.272 36.9432 -> norm = 0.0173205
setting bc 198.943 43.5712 223.078 -> 0 0 0
	 instead 198.933 43.5812 223.068 -> norm = 0.0173205
setting bc 198.943 43.5712 36.9332 -> 0 0 0
	 instead 198.933 43.5812 36.9432 -> norm = 0.0173205
setting bc 53.4216 178.282 223.078 -> 0 0 0
	 instead 53.4316 178.272 223.068 -> norm = 0.0173205
setting bc 53.4216 178.282 36.9332 -> 0 0 0
	 instead 53.4316 178.272 36.9432 -> norm = 0.0173205
setting bc 53.4216 43.5712 223.078 -> 0 0 0
	 instead 53.4316 43.5812 223.068 -> norm = 0.0173205
setting bc 53.4216 43.5712 36.9332 -> 0 0 0
	 instead 53.4316 43.5812 36.9432 -> norm = 0.0173205
 computing statistics for the displacement application error
 average norm of error in placement = 0.0173205
 iterating
 done with candidates
 Active BCs = 120405
 Total BCs = 482528
 no-eqs = 324624
	 percentage done= 0
	 percentage done= 5
	 percentage done= 10
	 percentage done= 15
	 percentage done= 20
	 percentage done= 25
	 percentage done= 30
	 percentage done= 35
	 percentage done= 40
	 percentage done= 45
	 percentage done= 50
	 percentage done= 55
	 percentage done= 60
	 percentage done= 65
	 percentage done= 70
	 percentage done= 75
	 percentage done= 80
	 percentage done= 85
	 percentage done= 90
	 percentage done= 95
 done_bc_natural size of container = 482528
 after building the map
 LOAD size = 24
 after setup_load_mfc
Linear solve converged due to CONVERGED_RTOL iterations 767
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=324624, cols=324624
    total: nonzeros=14973630, allocated nonzeros=14973963
      using I-node routines: found 108224 nodes, limit used is 5
 convergence reason = 2
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=324624, cols=324624
    total: nonzeros=14973630, allocated nonzeros=14973963
      using I-node routines: found 108224 nodes, limit used is 5
 Iterations 767
Absolute-Norm of error = 3.27848e-05
 countInvalid = 4 general-count = 482524
 Average of the error norm = 0.0837556
 Initial error = 0.406599
 Average of the error norm conditional = 0.0684362
 Active Residual error = 0.0171238
 Total error = 0.0838013
 c0324 time elapsed updating the sources (seconds) = 0
 orientation_pbs
 solving topology problems
 using old topology solver
 using Lame constants given by 
 E = 10
 nu = 0.3
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 198645
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 552
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 150543
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 456
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 157134
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 522
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 150416
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 477
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 192323
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 561
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 160838
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 510
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 165718
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 528
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 before 
 after
 !!!! No residual problems.
	 DONE
 compressing morph
dcds mesh direct image = 53.4216 43.5712 36.9332
 dir img for computed transform = 4.24431e-313 2.36145e-306 8.67399e-307
 dir img for compressed = 4.24431e-313 2.36145e-306 8.67399e-307
dcds total image for source point = 4.24431e-313 2.36145e-306 8.67399e-307
 name = compressed_step_16
   Initial transform 
	 name = compressed_step_17
   Initial transform 
	 name = 
 FEM error
 Initial displacement = 6.88518
 Final displacement = 6.12558
 Max after displacement = 20.359
 Error points = 0
 ======================
 step = 15
===============
elt_vol= 9.375
Constructing Delaunay tetrahedralization.
Delaunay seconds:  0.01
Creating surface mesh.
Perturbing vertices.
Delaunizing segments.
Constraining facets.
Segment and facet seconds:  0
Removing unwanted tetrahedra.
Hole seconds:  0
Repairing mesh.
Repair seconds:  0
Adding Steiner points to enforce quality.
Quality seconds:  15.5
Smoothing mesh.
Smooth seconds:  1.05

Writing nodes.
Writing elements.
Writing faces.

Output seconds:  0.09
Total running seconds:  16.65

Statistics:

  Input points: 8
  Input facets: 6
  Input holes: 0
  Input regions: 0

  Mesh points: 111645
  Mesh tetrahedra: 675953
  Mesh faces: 1367147
  Mesh subfaces: 30482
  Mesh subsegments: 630

Saving nodes to iteration.node
Saving elements to iteration.ele
Saving faces to iteration.face
 mesh nodes = 111645 mesh elts = 675953
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
	 count inserted = 100000 elapsed = 0.215 seconds  element count = 845031
	 count inserted = 200000 elapsed = 0.366 seconds  element count = 2564501
	 count inserted = 300000 elapsed = 0.335 seconds  element count = 4222033
	 count inserted = 400000 elapsed = 0.284 seconds  element count = 5749304
	 count inserted = 500000 elapsed = 0.273 seconds  element count = 7248398
	 count inserted = 600000 elapsed = 0.257 seconds  element count = 8744829
 done building octree - total elements = 9886946
 penalty_weight = 1000
setting bc 198.943 178.282 223.078 -> 0 0 0
	 instead 198.933 178.272 223.068 -> norm = 0.0173205
setting bc 198.943 178.282 36.9332 -> 0 0 0
	 instead 198.933 178.272 36.9432 -> norm = 0.0173205
setting bc 198.943 43.5712 223.078 -> 0 0 0
	 instead 198.933 43.5812 223.068 -> norm = 0.0173205
setting bc 198.943 43.5712 36.9332 -> 0 0 0
	 instead 198.933 43.5812 36.9432 -> norm = 0.0173205
setting bc 53.4216 178.282 223.078 -> 0 0 0
	 instead 53.4316 178.272 223.068 -> norm = 0.0173205
setting bc 53.4216 178.282 36.9332 -> 0 0 0
	 instead 53.4316 178.272 36.9432 -> norm = 0.0173205
setting bc 53.4216 43.5712 223.078 -> 0 0 0
	 instead 53.4316 43.5812 223.068 -> norm = 0.0173205
setting bc 53.4216 43.5712 36.9332 -> 0 0 0
	 instead 53.4316 43.5812 36.9432 -> norm = 0.0173205
 computing statistics for the displacement application error
 average norm of error in placement = 0.0173205
 iterating
 done with candidates
 Active BCs = 119202
 Total BCs = 482528
 no-eqs = 334935
	 percentage done= 0
	 percentage done= 5
	 percentage done= 10
	 percentage done= 15
	 percentage done= 20
	 percentage done= 25
	 percentage done= 30
	 percentage done= 35
	 percentage done= 40
	 percentage done= 45
	 percentage done= 50
	 percentage done= 55
	 percentage done= 60
	 percentage done= 65
	 percentage done= 70
	 percentage done= 75
	 percentage done= 80
	 percentage done= 85
	 percentage done= 90
	 percentage done= 95
 done_bc_natural size of container = 482528
 after building the map
 LOAD size = 24
 after setup_load_mfc
Linear solve converged due to CONVERGED_RTOL iterations 781
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=334935, cols=334935
    total: nonzeros=15455139, allocated nonzeros=15455490
      using I-node routines: found 111661 nodes, limit used is 5
 convergence reason = 2
KSP Object:
  type: cg
  maximum iterations=10000, initial guess is zero
  tolerances:  relative=1e-09, absolute=1e-50, divergence=10000
  left preconditioning
PC Object:
  type: jacobi
  linear system matrix = precond matrix:
  Matrix Object:
    type=seqaij, rows=334935, cols=334935
    total: nonzeros=15455139, allocated nonzeros=15455490
      using I-node routines: found 111661 nodes, limit used is 5
 Iterations 781
Absolute-Norm of error = 3.1232e-05
 countInvalid = 11 general-count = 482517
 Average of the error norm = 0.0855052
 Initial error = 0.408362
 Average of the error norm conditional = 0.0693114
 Active Residual error = 0.0171596
 Total error = 0.0856309
 c0324 time elapsed updating the sources (seconds) = 0
 orientation_pbs
 solving topology problems
 using old topology solver
 using Lame constants given by 
 E = 10
 nu = 0.3
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 179037
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 561
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 159960
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 498
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 202820
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 597
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 191474
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 624
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 176364
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 561
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 237718
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 645
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 173200
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 504
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 164222
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 540
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 142766
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 510
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 172129
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 552
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 122229
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 411
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
 done building octree - total elements = 171002
 solving cropped problem
 computing statistics for the displacement application error
 average norm of error in placement = 0
 LOAD size = 570
DirectSolverConvergence = ksp-converged-rtol
 btest = 1
 before 
 after
 !!!! No residual problems.
	 DONE
 compressing morph
dcds mesh direct image = 53.4216 43.5712 36.9332
 dir img for computed transform = 4.24431e-313 2.36145e-306 8.67399e-307
 dir img for compressed = 4.24431e-313 2.36145e-306 8.67399e-307
dcds total image for source point = 4.24431e-313 2.36145e-306 8.67399e-307
 name = compressed_step_15
   Initial transform 
	 name = compressed_step_16
   Initial transform 
	 name = compressed_step_17
   Initial transform 
	 name = 
 FEM error
 Initial displacement = 6.88518
 Final displacement = 5.74529
 Max after displacement = 19.0309
 Error points = 0
 ======================
 step = 14
===============
elt_vol= 9.0625
Constructing Delaunay tetrahedralization.
Delaunay seconds:  0
Creating surface mesh.
Perturbing vertices.
Delaunizing segments.
Constraining facets.
Segment and facet seconds:  0.01
Removing unwanted tetrahedra.
Hole seconds:  0
Repairing mesh.
Repair seconds:  0
Adding Steiner points to enforce quality.
Quality seconds:  16
Smoothing mesh.
Smooth seconds:  1.09

Writing nodes.
Writing elements.
Writing faces.

Output seconds:  0.09
Total running seconds:  17.19

Statistics:

  Input points: 8
  Input facets: 6
  Input holes: 0
  Input regions: 0

  Mesh points: 115284
  Mesh tetrahedra: 698895
  Mesh faces: 1413305
  Mesh subfaces: 31030
  Mesh subsegments: 633

Saving nodes to iteration.node
Saving elements to iteration.ele
Saving faces to iteration.face
 mesh nodes = 115284 mesh elts = 698895
 building index src
 done building the list
	 count inserted = 0 elapsed = 0 seconds  element count = 1
	 count inserted = 100000 elapsed = 0.228 seconds  element count = 819844
	 count inserted = 200000 elapsed = 0.374 seconds  element count = 2499279
	 count inserted = 300000 elapsed = 0.344 seconds  element count = 4134582
	 count inserted = 400000 elapsed = 0.285 seconds  element count = 5635920
	 count inserted = 500000 elapsed = 0.272 seconds  element count = 7107546
	 count inserted = 600000 elapsed = 0.251 seconds  element count = 8569514
 done building octree - total elements = 10026766
 penalty_weight = 1000
setting bc 198.943 178.282 223.078 -> 0 0 0
	 instead 198.933 178.272 223.068 -> norm = 0.0173205
setting bc 198.943 178.282 36.9332 -> 0 0 0
	 instead 198.933 178.272 36.9432 -> norm = 0.0173205
setting bc 198.943 43.5712 223.078 -> 0 0 0
	 instead 198.933 43.5812 223.068 -> norm = 0.0173205
setting bc 198.943 43.5712 36.9332 -> 0 0 0
	 instead 198.933 43.5812 36.9432 -> norm = 0.0173205
setting bc 53.4216 178.282 223.078 -> 0 0 0
	 instead 53.4316 178.272 223.068 -> norm = 0.0173205
setting bc 53.4216 178.282 36.9332 -> 0 0 0
	 instead 53.4316 178.272 36.9432 -> norm = 0.0173205
setting bc 53.4216 43.5712 223.078 -> 0 0 0
	 instead 53.4316 43.5812 223.068 -> norm = 0.0173205
setting bc 53.4216 43.5712 36.9332 -> 0 0 0
	 instead 53.4316 43.5812 36.9432 -> norm = 0.0173205
 computing statistics for the displacement application error
 average norm of error in placement = 0.0173205
 iterating
 done with candidates
 Active BCs = 117880
 Total BCs = 482528
 no-eqs = 345852
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Out of memory. This could be due to allocating
[0]PETSC ERROR: too large an object or bleeding by not properly
[0]PETSC ERROR: destroying unneeded objects.
[0]PETSC ERROR: Memory allocated 0 Memory used by process -2147483648
[0]PETSC ERROR: Try running with -malloc_dump or -malloc_log for info.
[0]PETSC ERROR: Memory requested 416405812!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 13, Thu May 15 17:29:26 CDT 2008 HG revision: 4466c6289a0922df26e20626fd4a0b4dd03c8124
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: surf2vol on a linux-gnu named penguin by root Fri Nov 12 10:59:45 2010
[0]PETSC ERROR: Libraries linked from /autofs/space/lyon_006/pubsw/Linux2-2.3-i386/packages/petsc/2.3.3-p13/src/petsc-2.3.3-p13/lib/linux-gnu-c-opt
[0]PETSC ERROR: Configure run at Mon Jun 16 12:09:46 2008
[0]PETSC ERROR: Configure options --with-debugging=0 --with-cc=gcc --with-fc=g77 --download-f-blas-lapack=0 --download-mpich=1 --with-mpi=1 --with-x=0 --with-gnu-copyright-code=0 --with-shared=0 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 FOPTFLAGS=-O3
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: PetscMallocAlign() line 61 in src/sys/memory/mal.c
[0]PETSC ERROR: MatSeqAIJSetPreallocation_SeqAIJ() line 2818 in src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: MatCreateSeqAIJ() line 2703 in src/mat/impls/aij/seq/aij.c
[0]PETSC ERROR: TSolver::setup_matrix() line 658 in unknowndirectory/solver.h
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal[0]PETSC ERROR: or try http://valgrind.org on linux or man libgmalloc on Apple to find memory corruption errors
[0]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run 
[0]PETSC ERROR: to get more information on the crash.
[0]PETSC ERROR: --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Signal received!
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 13, Thu May 15 17:29:26 CDT 2008 HG revision: 4466c6289a0922df26e20626fd4a0b4dd03c8124
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: surf2vol on a linux-gnu named penguin by root Fri Nov 12 10:59:45 2010
[0]PETSC ERROR: Libraries linked from /autofs/space/lyon_006/pubsw/Linux2-2.3-i386/packages/petsc/2.3.3-p13/src/petsc-2.3.3-p13/lib/linux-gnu-c-opt
[0]PETSC ERROR: Configure run at Mon Jun 16 12:09:46 2008
[0]PETSC ERROR: Configure options --with-debugging=0 --with-cc=gcc --with-fc=g77 --download-f-blas-lapack=0 --download-mpich=1 --with-mpi=1 --with-x=0 --with-gnu-copyright-code=0 --with-shared=0 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 FOPTFLAGS=-O3
[0]PETSC ERROR: ------------------------------------------------------------------------
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown file
[unset]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0