###################################
## Tests for AnalyzeFMRI package ##
###################################

library(AnalyzeFMRI, quietly = TRUE)

## Test that reading of examples works

a1 <- f.read.analyze.volume(system.file("example.img", package="AnalyzeFMRI"))

a1[30:40, 30:40, 10, 1]


## Test that writing of examples works

a2 <- array(1:1000, dim = c(10, 10, 10, 1))
tmpdir <- tempdir()
f.write.analyze(a2, file = "test.array", size = "float", path.out = tmpdir)

a3 <- f.read.analyze.volume(file = paste0(tmpdir, "/test.array.img"))
unlink(tmpdir) # tidy up

sum(a2 != a3) ## should return 0