[cig-commits] r6872 - in short/3D/PyLith/trunk: . pylith/meshio unittests/pytests/feassemble unittests/pytests/topology

[brad at geodynamics.org]

Author: brad
Date: 2007-05-13 13:25:40 -0700 (Sun, 13 May 2007)
New Revision: 6872

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/pylith/meshio/SolutionIO.py
   short/3D/PyLith/trunk/unittests/pytests/feassemble/TestQuadrature.py
   short/3D/PyLith/trunk/unittests/pytests/topology/TestMeshGenSimple.py
Log:
Fixed a few small bugs in SolutionIO.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2007-05-13 04:32:31 UTC (rev 6871)
+++ short/3D/PyLith/trunk/TODO	2007-05-13 20:25:40 UTC (rev 6872)
@@ -15,6 +15,10 @@
 
 1. Implement cell geometry stuff.
    Compute Jacobian in cell.
+
+   Questions: 
+     (1) Need to get coordinates of "nodes" in reference cell.
+
    b. C++ unit tests
 
 2. Implement Dirichlet boundary conditions

Modified: short/3D/PyLith/trunk/pylith/meshio/SolutionIO.py
===================================================================
--- short/3D/PyLith/trunk/pylith/meshio/SolutionIO.py	2007-05-13 04:32:31 UTC (rev 6871)
+++ short/3D/PyLith/trunk/pylith/meshio/SolutionIO.py	2007-05-13 20:25:40 UTC (rev 6872)
@@ -75,7 +75,7 @@
       raise ValueError, "Coordinate system for output is unknown."
     self.coordsys.initialize()
 
-    assert(cppHandle != None)
+    assert(self.cppHandle != None)
     self.cppHandle.open(mesh.cppHandle)
     return
 
@@ -116,7 +116,7 @@
     """
     Force synchronization between Python and C++.
     """
-    self.cppHandle.coordsys = self.coordsys
+    self.cppHandle.coordsys = self.coordsys.cppHandle
     return
 
 

Modified: short/3D/PyLith/trunk/unittests/pytests/feassemble/TestQuadrature.py
===================================================================
--- short/3D/PyLith/trunk/unittests/pytests/feassemble/TestQuadrature.py	2007-05-13 04:32:31 UTC (rev 6871)
+++ short/3D/PyLith/trunk/unittests/pytests/feassemble/TestQuadrature.py	2007-05-13 20:25:40 UTC (rev 6872)
@@ -90,7 +90,6 @@
     quadrature = Quadrature1D()
     quadrature.cell = cell
 
-    print cell.basisVert
     quadrature.initialize()
 
     from pylith.utils.testarray import test_double

Modified: short/3D/PyLith/trunk/unittests/pytests/topology/TestMeshGenSimple.py
===================================================================
--- short/3D/PyLith/trunk/unittests/pytests/topology/TestMeshGenSimple.py	2007-05-13 04:32:31 UTC (rev 6871)
+++ short/3D/PyLith/trunk/unittests/pytests/topology/TestMeshGenSimple.py	2007-05-13 20:25:40 UTC (rev 6872)
@@ -40,16 +40,16 @@
     generator.interpolate = True
     generator.setBoundary(generator.createCubeBoundary())
     mesh = generator.create()
-    
-    #from pylith.meshio.SolutionIOVTK import SolutionIOVTK
-    #io = SolutionIOVTK()
-    #io.filename = 'mesh.vtk'
-    #import spatialdata.geocoords.CSCart
-    #io.coordsys = spatialdata.geocoords.CSCart.CSCart()
-    #io.open(mesh)
-    #io.writeTopology(mesh)
-    #io.close()
 
+    from pylith.meshio.SolutionIOVTK import SolutionIOVTK
+    io = SolutionIOVTK()
+    io.filename = 'mesh.vtk'
+    from spatialdata.geocoords.CSCart import CSCart
+    io.coordsys = CSCart()
+    io.open(mesh)
+    io.writeTopology(mesh)
+    io.close()
+
     #from pylith.topology.Distributor import Distributor
     #distributor = Distributor()
     #newMesh = distributor.distribute(mesh)