From f19eec50b9c8cd5850dd1713d3b763a7c02f3127 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Sat, 18 Apr 2009 18:34:17 +0000
Subject: [PATCH] make binary writing of iterates work for all TargetSpace
 types

[[Imported from SVN: r4050]]
---
 src/riemanniantrsolver.cc | 11 ++---------
 1 file changed, 2 insertions(+), 9 deletions(-)

diff --git a/src/riemanniantrsolver.cc b/src/riemanniantrsolver.cc
index 7b09fe79..0e386667 100644
--- a/src/riemanniantrsolver.cc
+++ b/src/riemanniantrsolver.cc
@@ -383,15 +383,8 @@ void RiemannianTrustRegionSolver<GridType,TargetSpace>::solve()
             if (!fpRod)
                 DUNE_THROW(SolverError, "Couldn't open file " << iRodFilename << " for writing");
             
-            for (int j=0; j<x_.size(); j++) {
-
-                for (int k=0; k<3; k++)
-                    fwrite(&x_[j].r[k], sizeof(double), 1, fpRod);
-
-                for (int k=0; k<4; k++)  // 3d hardwired here!
-                    fwrite(&x_[j].q[k], sizeof(double), 1, fpRod);
-
-            }
+            for (int j=0; j<x_.size(); j++)
+                fwrite(&x_[j], sizeof(TargetSpace), 1, fpRod);
 
             fclose(fpRod);
 
-- 
GitLab