Write the output file in main.

diff --git a/src/Main.hs b/src/Main.hs
index 7b50de9ee3c02fab3e40a7d662198f4fca531a41..b8edf9eb2ad1b6d3e88f4641176aaf291aab6912 100644 (file)
--- a/src/Main.hs
+++ b/src/Main.hs
@@ -18,4 +18,4 @@ main = do
   mridata <- read_values_3d mri_shape "./data/mridata.txt"
   let g = make_grid 1 mridata
   let output = zoom g 1
-  print output
+  write_values_1d output "output.txt"

diff --git a/src/Values.hs b/src/Values.hs
index bd7c4d66af25ba443a6afe7f72d7bc6885df7f3b..a33bbba56d3eef9da970d845de285af4c39ac39d 100644 (file)
--- a/src/Values.hs
+++ b/src/Values.hs
@@ -14,9 +14,11 @@ import Data.Array.Repa (
   fromList,
   index,
   reshape,
+  size
   )
 
-import Data.Array.Repa.IO.Vector (readVectorFromTextFile)
+import Data.Array.Repa.IO.Vector (readVectorFromTextFile,
+                                  writeVectorToTextFile)
 import System.FilePath ()
 import Test.QuickCheck (Arbitrary(..), Gen)
 
@@ -54,6 +56,12 @@ read_values_3d sh path = do
   one_d <- read_values_1d path
   return $ reshape sh one_d
 
+write_values_1d :: Values3D -> FilePath -> IO ()
+write_values_1d v3d path = do
+  let size3d = size $ extent v3d
+  let shape1d = (Z :. size3d)
+  let v1d = reshape shape1d v3d
+  writeVectorToTextFile v1d path
 
 empty3d :: Values3D
 empty3d = Data.Array.Repa.fromList (Z :. 0 :. 0 :. 0) []