From a961a1f48b85615712eaa35de7feab3db8a1bf50 Mon Sep 17 00:00:00 2001
From: Jukka Aho <ahojukka5@gmail.com>
Date: Sat, 20 Jun 2015 21:47:54 +0300
Subject: [PATCH] xdmf test

---
 test/runtests.jl  |  2 ++
 test/test_xdmf.jl | 33 +++++++++++++++++++++++++++++++++
 2 files changed, 35 insertions(+)
 create mode 100644 test/test_xdmf.jl

diff --git a/test/runtests.jl b/test/runtests.jl
index d221849..e9d0bdc 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -125,6 +125,8 @@ test_get_element()
 
 
 
+include("test_xdfm.jl")
+
 # write your own tests here
 # @test 1 == JuliaFEM.test()
 # @test_approx_eq 1.0 1.0
diff --git a/test/test_xdmf.jl b/test/test_xdmf.jl
new file mode 100644
index 0000000..65d56b8
--- /dev/null
+++ b/test/test_xdmf.jl
@@ -0,0 +1,33 @@
+using JuliaFEM
+using Base.Test
+
+
+function test_write_to_xdfm()
+    m = new_model()
+    nodes = Dict(1 => [0.0, 0.0, 0.0],
+                 2 => [1.0, 0.0, 0.0],
+                 3 => [0.0, 1.0, 0.0],
+                 4 => [0.0, 0.0, 1.0])
+    add_nodes(m, nodes)
+    el1 = Dict("element_type" => 0x6, "node_ids" => [1, 2, 3, 4])
+    elements = [el1]
+    add_elements(elements)
+
+    # create new field "temperature" for nodal points
+    f = new_field(m.nodes, "temperature")
+    f[1] = 0.1
+    f[2] = 0.2
+    f[3] = 0.3
+    # do NOT write to fourth node.
+
+    # export model to xdmf
+    fn = tempname()
+    println("temp file: ", fn)
+    time = 123
+    write_xdmf(m, fn, time; fields = ["temperature"])
+
+    d = h5open(fn)
+    @test d[123]["temperature"][:] = [0.1, 0.2, 0.3]
+
+end
+test_write_to_xdfm()
\ No newline at end of file