int read_spinor(spinor * const s, spinor * const r, char * filename, const int position_) {
int status = 0, getpos = 0, prec = 0, prop_type, position = position_, rstat=0;
uint64_t bytes = 0;
char * header_type = NULL;
READER * reader = NULL;
DML_Checksum checksum_read;
DML_Checksum checksum;
char * checksum_string = NULL;
int DML_read_flag = 0;
construct_reader(&reader, filename);
/* determine the propagator type */
prop_type = parse_propagator_type(reader);
switch (prop_type) {
case 1:
/* strictly speaking the following depends on whether we read a source or a propagator */
position = 2 * position_ +1;
break;
case 2:
case 3:
return(-2);
case 11:
case 12:
case 13:
return(-3);
case -1:
case 4:
prop_type = 0;
}
/* seek back to beginning of file*/
destruct_reader(reader);
construct_reader(&reader, filename);
/* Find the desired propagator (could be more than one in a file) */
while ((status = ReaderNextRecord(reader)) != LIME_EOF) {
if (status != LIME_SUCCESS) {
fprintf(stderr, "ReaderNextRecord returned status %d.\n", status);
break;
}
header_type = ReaderType(reader);
if (strcmp("scidac-binary-data", header_type) == 0) {
if (getpos == position) {
break;
}
++getpos;
}
}
if (status == LIME_EOF) {
fprintf(stderr, "Unable to find requested LIME record scidac-binary-data in file %s.\nEnd of file reached before record was found.\n", filename);
return(-5);
}
bytes = ReaderBytes(reader);
if ((uint64_t)bytes == LX * g_nproc_x * LY * g_nproc_y * LZ * g_nproc_z * T * g_nproc_t * sizeof(spinor)) {
prec = 64;
}
else {
if ((uint64_t)bytes == LX * g_nproc_x * LY * g_nproc_y * LZ * g_nproc_z * T * g_nproc_t * sizeof(spinor) / 2) {
prec = 32;
}
else {
fprintf(stderr, "Length of scidac-binary-data record in %s does not match input parameters.\n", filename);
fprintf(stderr, "Found %lu bytes.\n", bytes);
return(-6);
}
}
if (g_cart_id == 0 && g_debug_level >= 0) {
printf("# %s precision read (%d bits).\n", (prec == 64 ? "Double" : "Single") ,prec);
}
if(r == NULL) {
if( (rstat = read_binary_spinor_data_l(s, reader, &checksum)) != 0) {
fprintf(stderr, "read_binary_spinor_data_l failed with return value %d", rstat);
return(-7);
}
}
else {
if( (rstat = read_binary_spinor_data(s, r, reader, &checksum)) != 0) {
fprintf(stderr, "read_binary_spinor_data failed with return value %d", rstat);
return(-7);
}
}
if( g_disable_src_IO_checks != 1 ){
// we search for a scidac-checksum directly after the binary data
// but only until more binary data is found
while ((status = ReaderNextRecord(reader)) != LIME_EOF) {
if (status != LIME_SUCCESS) {
fprintf(stderr, "ReaderNextRecord returned status %d.\n", status);
break;
}
header_type = ReaderType(reader);
if (strcmp("scidac-checksum", header_type) == 0) {
read_message(reader, &checksum_string);
DML_read_flag = parse_checksum_xml(checksum_string, &checksum_read);
free(checksum_string);
break;
}
if (strcmp("scidac-binary-data", header_type) == 0 || strcmp("ildg-binary-data", header_type) == 0) {
break;
}
}
if (!DML_read_flag) {
fprintf(stderr, "LIME record with name: \"scidac-checksum\", in gauge file %s either missing or malformed.\n", filename);
fprintf(stderr, "Unable to verify integrity of gauge field data.\n");
return(-1);
}
if (g_cart_id == 0 && g_debug_level >= 0) {
printf("# Scidac checksums for DiracFermion field %s position %d:\n", filename, position);
printf("# Calculated : A = %#010x B = %#010x.\n", checksum.suma, checksum.sumb);
printf("# Read from LIME headers: A = %#010x B = %#010x.\n", checksum_read.suma, checksum_read.sumb);
}
}
destruct_reader(reader);
return(0);
}
int read_spinor(spinor * const s, spinor * const r, char * filename, const int position_) {
int status = 0, getpos = 0, bytes = 0, prec = 0, prop_type, position = position_, rstat=0;
char *header_type = NULL;
READER *reader = NULL;
DML_Checksum checksum;
construct_reader(&reader, filename);
/* determine the propagator type */
prop_type = parse_propagator_type(reader);
switch (prop_type) {
case 1:
/* strictly speaking the following depends on whether we read a source or a propagator */
position = 2 * position_ +1;
break;
case 2:
case 3:
return(-2);
case 11:
case 12:
case 13:
return(-3);
case -1:
case 4:
prop_type = 0;
}
/* seek back to beginning of file*/
destruct_reader(reader);
construct_reader(&reader, filename);
/* Find the desired propagator (could be more than one in a file) */
while ((status = ReaderNextRecord(reader)) != LIME_EOF) {
if (status != LIME_SUCCESS) {
fprintf(stderr, "ReaderNextRecord returned status %d.\n", status);
break;
}
header_type = ReaderType(reader);
if (strcmp("scidac-binary-data", header_type) == 0) {
if (getpos == position) {
break;
}
++getpos;
}
}
if (status == LIME_EOF) {
fprintf(stderr, "Unable to find requested LIME record scidac-binary-data in file %s.\nEnd of file reached before record was found.\n", filename);
return(-5);
}
bytes = ReaderBytes(reader);
if ((int)bytes == LX * g_nproc_x * LY * g_nproc_y * LZ * g_nproc_z * T * g_nproc_t * sizeof(spinor)) {
prec = 64;
}
else {
if ((int)bytes == LX * g_nproc_x * LY * g_nproc_y * LZ * g_nproc_z * T * g_nproc_t * sizeof(spinor) / 2) {
prec = 32;
}
else {
fprintf(stderr, "Length of scidac-binary-data record in %s does not match input parameters.\n", filename);
fprintf(stderr, "Found %d bytes.\n", bytes);
return(-6);
}
}
if (g_cart_id == 0 && g_debug_level >= 0) {
printf("# %s precision read (%d bits).\n", (prec == 64 ? "Double" : "Single") ,prec);
}
if(r == NULL) {
if( (rstat = read_binary_spinor_data_l(s, reader, &checksum)) != 0) {
fprintf(stderr, "read_binary_spinor_data_l failed with return value %d", rstat);
return(-7);
}
}
else {
if( (rstat = read_binary_spinor_data(s, r, reader, &checksum)) != 0) {
fprintf(stderr, "read_binary_spinor_data failed with return value %d", rstat);
return(-7);
}
}
if (g_cart_id == 0 && g_debug_level >= 0) {
printf("# Scidac checksums for DiracFermion field %s position %d:\n", filename, position);
printf("# Calculated : A = %#x B = %#x.\n", checksum.suma, checksum.sumb);
printf("# No Scidac checksum was read from headers, unable to check integrity of file.\n");
}
destruct_reader(reader);
return(0);
}