TREX/trexio
1
0
Fork 0
mirror of https://github.com/TREX-CoE/trexio.git synced 2024-12-23 04:43:57 +01:00
Code Issues Projects Releases Wiki Activity

consecutive testing for hdf5

Browse Source
This commit is contained in:
q-posev 2021-03-02 11:41:25 +01:00
parent cbb936564e
commit a0172e26c1
2 changed files with 5 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
src/Makefile
View File

@ -44,7 +44,6 @@ libtrio.so: $(OBJECT_FILES) $(HEADER_FILES)
test: libtrio.so test.c test: libtrio.so test.c
$(CC) $(CFLAGS) $(INCLUDE) -Wl,-rpath,$(PWD) -L. test.c -ltrio $(LIBS) -o test $(CC) $(CFLAGS) $(INCLUDE) -Wl,-rpath,$(PWD) -L. test.c -ltrio $(LIBS) -o test
rm test_write.h5
clean: clean:
rm -f *.o libtrio.so test test_write.h5 rm -f *.o libtrio.so test test_write.h5

20
src/test.c
View File

@ -29,8 +29,6 @@ int test_h5write() {
uint64_t num = 12; uint64_t num = 12;
//double charge[12] = {6., 6., 6., 6., 6., 6., 1., 1., 1., 1., 1., 1.};
double coord[36] = { double coord[36] = {
0.00000000 , 1.39250319 , 0.00000000 , 0.00000000 , 1.39250319 , 0.00000000 ,
-1.20594314 , 0.69625160 , 0.00000000 , -1.20594314 , 0.69625160 , 0.00000000 ,
@ -50,15 +48,14 @@ int test_h5write() {
// works: try writing info in an empty file // works: try writing info in an empty file
rc = trio_write_nucleus_num(file,num); rc = trio_write_nucleus_num(file,num);
// rc = trio_write_nucleus_charge(file,charge);
rc = trio_write_nucleus_coord(file,coord); rc = trio_write_nucleus_coord(file,coord);
// should not work: try to rewrite the nucleus_num // should not work: try to rewrite the nucleus_num
rc = trio_write_nucleus_num(file,25); rc = trio_write_nucleus_num(file,25);
// works: try to rewrite the nucleus_coord // works: try to rewrite the nucleus_coord
// coord[0] = 666.0; coord[0] = 666.666;
// rc = trio_write_nucleus_coord(file,coord); rc = trio_write_nucleus_coord(file,coord);
if (rc == TRIO_SUCCESS) { if (rc == TRIO_SUCCESS) {
@ -73,25 +70,18 @@ int test_h5write() {
} }
int test_h5read() { int test_h5read() {
const char* file_name = "test.h5"; const char* file_name = "test_write.h5";
trio_t* file = NULL; trio_t* file = NULL;
trio_exit_code rc; trio_exit_code rc;
uint64_t num; uint64_t num;
//double* charge;
double* coord; double* coord;
file = trio_create(file_name, TRIO_HDF5); file = trio_create(file_name, TRIO_HDF5);
rc = trio_read_nucleus_num(file,&num); rc = trio_read_nucleus_num(file,&num);
assert (num == 4); assert (num == 12);
/*
charge = (double*) calloc(num, sizeof(double));
rc = trio_read_nucleus_charge(file,charge);
assert(charge[] == 1.);
*/
coord = (double*) calloc(3*num, sizeof(double)); coord = (double*) calloc(3*num, sizeof(double));
rc = trio_read_nucleus_coord(file,coord); rc = trio_read_nucleus_coord(file,coord);
@ -100,7 +90,7 @@ int test_h5read() {
printf("%lf \n", coord[i]); printf("%lf \n", coord[i]);
}*/ }*/
double x = coord[0] - 1.2; double x = coord[30] - 2.14171677;
assert( x*x < 1.e-12); assert( x*x < 1.e-12);
if (rc == TRIO_SUCCESS) { if (rc == TRIO_SUCCESS) {