int mumps_io_open_files_for_read(){
int i,j;
mumps_file_struct *mumps_io_pfile_pointer_array;
#ifdef IRIX64_
struct dioattr dio;
#endif
for(j=0;j<mumps_io_nb_file_type;j++){
mumps_io_pfile_pointer_array=(mumps_files+j)->mumps_io_pfile_pointer_array;
for(i=0;i<(mumps_files+j)->mumps_io_nb_file;i++){
#if ! defined( MUMPS_WIN32 )
(mumps_io_pfile_pointer_array+i)->file=open((mumps_io_pfile_pointer_array+i)->name,(mumps_files+j)->mumps_flag_open);
if((mumps_io_pfile_pointer_array+i)->file==-1){
return mumps_io_sys_error(-90,"Problem while opening OOC file");
}
#else
(mumps_io_pfile_pointer_array+i)->file=fopen((mumps_io_pfile_pointer_array+i)->name,(mumps_files+j)->mumps_flag_open);
if((mumps_io_pfile_pointer_array+i)->file==NULL){
return mumps_io_error(-90,"Problem while opening OOC file");
}
(mumps_io_pfile_pointer_array+i)->is_opened=1;
#endif
}
}
return 0;
}
int mumps_io_read_os_buff__(void * file,void * loc_addr,size_t size,int local_offset){
size_t ret_code;
/* printf("Read with buff %d %d %d\n",(int) size, local_offset,*((int *)file)); */
# ifdef WITH_PFUNC
ret_code=pread(*(int *)file,loc_addr,size,local_offset);
# else
lseek(*(int *)file,(long) local_offset,SEEK_SET);
ret_code=read(*(int *)file,loc_addr,size);
# endif
if((int) ret_code==-1){
return mumps_io_sys_error(-90,"Problem with low level read");
}
return 0;
}
void MUMPS_CALL
MUMPS_OOC_REMOVE_FILE_C(MUMPS_INT *ierr, char *name, mumps_ftnlen l1)
{
char buf[296]; /* for error message, count 256 chars for name */
*ierr=(MUMPS_INT)remove(name);
if(*ierr<0){
#if ! defined(MUMPS_WIN32)
sprintf(buf,"Unable to remove OOC file %s",name);
#else
sprintf(buf,"Unable to remove OOC file %s with return value %d",name,*ierr);
#endif
*ierr = -90;
mumps_io_sys_error((int)*ierr,buf);
return;
}
return;
}
int mumps_io_write_os_buff__(void *file, void *loc_addr, size_t write_size, int where){
size_t ret_code;
/* printf("write with buff %d %d %d\n",(int) write_size, where,*((int *)file)); */
# ifdef WITH_PFUNC
ret_code=pwrite(*(int *)file,loc_addr,write_size,where);
# else
/*in this case all the I/O's are made by the I/O thread => we don't
need to protect the file pointer.*/
lseek(*(int *)file,(long)where,SEEK_SET);
ret_code=write(*(int *)file,loc_addr,write_size);
# endif
if((int)ret_code==-1){
return mumps_io_sys_error(-90,"Problem with low level write");
} else if(ret_code!=write_size){
return mumps_io_error(-90,"Error not enough space on disk \n");
}
return 0;
}
int mumps_free_file_pointers(int *step){
int i,j,bound,ierr;
/* Free prefix only for facto */
if (*step == 0) free(mumps_ooc_file_prefix);
if(mumps_files == NULL )
return 0;
#if ! defined( MUMPS_WIN32 )
#endif
bound=mumps_io_nb_file_type;
/* if(*step==0){ */
/* /\* factorization *\/ */
/* bound=NB_FILE_TYPE_FACTO; */
/* }else{ */
/* /\* solve *\/ */
/* bound=NB_FILE_TYPE_SOLVE; */
/* } */
for(j=0;j<bound;j++){
if( mumps_files[j].mumps_io_pfile_pointer_array == NULL ) {
continue;
}
for(i=0;i<(mumps_files+j)->mumps_io_nb_file_opened;i++){
#if ! defined( MUMPS_WIN32 )
ierr=close((((mumps_files+j)->mumps_io_pfile_pointer_array)+i)->file);
if(ierr==-1){
return mumps_io_sys_error(-90,"Problem while closing OOC file");
}
#else
ierr=fclose((((mumps_files+j)->mumps_io_pfile_pointer_array)+i)->file);
if(ierr==-1){
return mumps_io_error(-90,"Problem while closing OOC file\n");
}
#endif
/* free(*(mumps_io_pfile_name+i)); */
}
free((mumps_files+j)->mumps_io_pfile_pointer_array);
}
/* free(mumps_io_pfile_name); */
free(mumps_files);
#if ! defined( MUMPS_WIN32 )
#endif
return 0;
}
/* Functions */
int mumps_set_file(int type,int file_number_arg){
/* Defines the pattern for the file name. The last 6 'X' will be replaced
so as to name were unique */
char name[351];
#if ! defined(_WIN32)
int fd;
char buf[64]; /* for error message */
#endif
mumps_file_struct *mumps_io_pfile_pointer_array;
/* if ((mumps_files+type)->mumps_io_current_file_number >= ((mumps_files+type)->mumps_io_nb_file)-1){*/
if (file_number_arg > ((mumps_files+type)->mumps_io_nb_file)-1){
/* Exception : probably thrown because of a bad estimation
of number of files. */
/* We increase the number of file needed and then realloc. */
((mumps_files+type)->mumps_io_nb_file)++;
(mumps_files+type)->mumps_io_pfile_pointer_array=realloc((void *)(mumps_files+type)->mumps_io_pfile_pointer_array,((mumps_files+type)->mumps_io_nb_file)*sizeof(mumps_file_struct));
if((mumps_files+type)->mumps_io_pfile_pointer_array==NULL){
return mumps_io_error(-13,"Allocation problem in low-level OOC layer\n");
}
((mumps_files+type)->mumps_io_pfile_pointer_array+((mumps_files+type)->mumps_io_nb_file)-1)->is_opened = 0;
}
mumps_io_pfile_pointer_array=(mumps_files+type)->mumps_io_pfile_pointer_array;
/* Do change the current file */
((mumps_files+type)->mumps_io_current_file_number)=file_number_arg;
if((mumps_io_pfile_pointer_array+file_number_arg)->is_opened!=0){
((mumps_files+type)->mumps_io_current_file_number)=file_number_arg;
return 0;
}
/* #if ! defined( MUMPS_WIN32 )*/
/* MinGW does not have a mkstemp function and MinGW defines _WIN32,
* so we also go in the else branch below with MinGW */
#if ! defined(_WIN32)
strcpy(name,mumps_ooc_file_prefix);
fd=mkstemp(name);
/* Note that a file name is built by mkstemp and that the file is
opened. fd hold the file descriptor to access it.
We want to close the file that will be opened later
and might be removed before the end of the processus.
*/
if(fd < 0) {
sprintf(buf,"mkstemp(%s) failed\n", mumps_ooc_file_prefix);
return mumps_io_sys_error(-99,buf);
} else {
close(fd);
}
#else
sprintf(name,"%s_%d_%d",mumps_ooc_file_prefix,((mumps_files+type)->mumps_io_current_file_number)+1,type);
#endif
/* *(mumps_io_pfile_pointer_array+mumps_io_current_file_number)=fopen(name,"w+"); */
/* *(mumps_io_pfile_name+mumps_io_current_file_number)=(char *)malloc((strlen(name)+1)*sizeof(char)); */
/* if(*(mumps_io_pfile_name+mumps_io_current_file_number)==NULL){ */
/* sprintf(error_str,"Allocation problem in low-level OOC layer\n"); */
/* return -13; */
/* } */
strcpy((mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->name,name);
/* See mumps_io_basic.h for comments on the I/O flags passed to open */
#if ! defined( MUMPS_WIN32 )
(mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->file=open(name,(mumps_files+type)->mumps_flag_open,0666);
/*
CPA: for LU factor file:
(mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->file= open(name, O_WRONLY | O_CREAT | O_TRUNC, 0666); */
if((mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->file==-1){
return mumps_io_sys_error(-90,"Unable to open OOC file");
}
#else
(mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->file=fopen(name,(mumps_files+type)->mumps_flag_open);
if((mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number)->file==NULL){
return mumps_io_error(-90,"Problem while opening OOC file");
}
#endif
(mumps_files+type)->mumps_io_current_file=(mumps_io_pfile_pointer_array+(mumps_files+type)->mumps_io_current_file_number);
((mumps_files+type)->mumps_io_nb_file_opened)++;
if((mumps_files+type)->mumps_io_current_file_number>(mumps_files+type)->mumps_io_last_file_opened){
(mumps_files+type)->mumps_io_last_file_opened=(mumps_files+type)->mumps_io_current_file_number;
}
/* if(*(mumps_io_pfile_pointer_array+mumps_io_current_file_number)==NULL) */
((mumps_files+type)->mumps_io_current_file)->write_pos=0;
((mumps_files+type)->mumps_io_current_file)->is_opened=1;
/* printf("new file created -> num = %d \n", ((mumps_files+type)->mumps_io_last_file_opened));*/
/* printf("new file created %d\n",mumps_io_current_file_number);*/
return 0;
}