diff --git a/rod3d.cc b/rod3d.cc
index 3664cc7ee4b3ca193f0711072cdf959e240432d3..b00351f2e4113af4c0e7d78bd385f14bc658e679 100644
--- a/rod3d.cc
+++ b/rod3d.cc
@@ -124,6 +124,7 @@ int main (int argc, char *argv[]) try
     RodAssembler<GridType> rodAssembler(grid, &localStiffness);
 
     RiemannianTrustRegionSolver<GridType,RigidBodyMotion<3> > rodSolver;
+#if 0
     rodSolver.setup(grid, 
                     &rodAssembler,
                     x,
@@ -137,6 +138,18 @@ int main (int argc, char *argv[]) try
                     baseIterations,
                     baseTolerance,
                     instrumented);
+#else
+    rodSolver.setupTCG(grid, 
+                       &rodAssembler,
+                       x,
+                       dirichletNodes,
+                       tolerance,
+                       maxTrustRegionSteps,
+                       initialTrustRegionRadius,
+                       multigridIterations,
+                       mgTolerance,
+                       instrumented);
+#endif
 
     // /////////////////////////////////////////////////////
     //   Solve!