Department of Scientific Computing    Institute for Numerical Simulation    University of Bonn

Documentation

Download

Programming References

Bug Reports / Suggestions

FAQ

Authors

           

---

UDFInformation.hpp

//
// collects some function concerning basic
// informations about UDF files
//

#ifndef _UDF_INFORMATION_
# define _UDF_INFORMATION_

# include <stdio.h>
# include <stdlib.h>
# include <string.h>

#include "UDF.h"
#include "Math.hpp"
#include "Extensions.hpp"


template<int DIM>
struct UDFInformation {
  double *X  [20] ; // X[i][0], ... ,X[size[i]-1]: Gridlines for ith coord.dir. 
  int    size[20],  // 
         nel     ;  // size[0]*..*size[DIM-1], actual dimension

  double mi,ma ; // min/max of scalar data

  bool ascii,do_swap,scan,read ;
  FILE *lastfile ;

    UDFInformation() ;
   ~UDFInformation() ;

    void ReadHeader(FILE *fid, Extensions<DIM> *E) ; // fid must be pointer to an open file
    void Scan      (FILE *fid) ;
    void Print     () ;
} ;


//  
//  Implementations
//

template<int DIM>
UDFInformation<DIM>::UDFInformation() {
 for (size_t i=0; i<sizeof(size)/sizeof(size[0]); i++) {
   X[i]   =NULL ;
   size[i]=0 ;
 }
 ma=mi=0 ;
 
 do_swap= (UDF_CPUisBigEndian()==0) ;
 read=false;
 lastfile=NULL ;
 scan=false ; 
}

template<int DIM>
UDFInformation<DIM>::~UDFInformation() {
 for (size_t i=0; i<sizeof(size)/sizeof(size[0]); i++)
   if (X[i]!=NULL) delete []X[i] ;
}


//
// Print UDFinformation
//
template<int DIM>
void UDFInformation<DIM>::Print() {
 int i,j ;
 if (!read) {
   std::cout<< "no header read\n" ;
   exit(-1) ;
 }
 std::cout<< "UDF file statistics\n" ;
 std::cout<< "number of elements: "<<nel<<" = "<<size[0] ;
 for (i=1; i<DIM; i++) std::cout<<'x'<<size[i] ;
 std::cout<<"\n" ;
 
 if (ascii) std::cout<<"ascii\n"; else std::cout<<"binary\n" ;
 if (scan) {
   std::cout<<"Data scaned\n" ;
   std::cout<<"min="<<mi<<" , max="<<ma<<"\n" ;
 }
}

//
// read UDF header
//
template<int DIM>
void UDFInformation<DIM>::ReadHeader(FILE *fid, Extensions<DIM> *Ext) {
 int i,si,j,N,TypeID, st1,st2 ;
 char s[1000],fd[100],c ;
 float  fx ;
 double x  ;

 read=true  ;
 lastfile=fid   ;
 scan=false ;

 fscanf(fid,"%s %s %s %s %s",s,s,s,s,s) ; // # vtk ....
 fscanf(fid,"%s",s) ;                     // UDF-Data 
 // read extensions:
 if (strcmp(s,"UDF-DataSetEXT00")==0) {
   fscanf(fid,"%d %d %d %d",&N,&TypeID,&st1,&st2) ;
   bool ms=false ;
   for (i=0; i<DIM; i++) {
     si=1<<i ;
     Ext->islifting   [i] =  ((st1&si) != 0);
     Ext->ismultiscale[i] =  ((st2&si) != 0);
     if (Ext->ismultiscale[i]) ms=true ;    
     fscanf(fid,"%d %d", &Ext->BC[i][0], &Ext->BC[i][1]) ;
     if ( (Ext->BC[i][1]!=-1) && (Ext->BC[i][1]!=0) && (Ext->BC[i][1]!=1) ) Ext->BC[i][1]=PERIODIC ;
   }

   // if multiscale representation, then check consistency of wavelet basis 
   if (ms && ( (N!=Ext->WC->N)||(TypeID!=Ext->WC->TypeID(Ext->WC->type)) ) ) {
     std::cout<<"File to read contains wavelet coefficients w.r.t wavelets with ID("<<N<<','<<TypeID<<") but current wavelets have ID (";
     std::cout<<Ext->WC->N<<','<<Ext->WC->TypeID(Ext->WC->type)<<"). This may cause wrong results\n" ;
     assert(0) ;
   }

 } else {
   for (i=0; i<DIM; i++) {
     Ext->BC[i][0]=Ext->BC[i][1] = -MAXINT ;
     Ext->islifting   [i]=false ;
     Ext->ismultiscale[i]=true  ;
   }
 }
 
 fscanf(fid,"%s",s) ;                     // ASCII/BINARY

 ascii=strcmp(s,"ASCII")==0 ;

 fscanf(fid,"%s %s",s,s) ; // DATASET RECTI....
 fscanf(fid,"%s",s) ;      // DIMENSION

 // Dimensions
 nel=1 ;
 for (i=0; i<max(DIM , 3); i++) {
   fscanf(fid,"%d",&si) ;

   if (size[i]>0) {
     if (X[i]==NULL) { std::cout<<"size="<<size[i]<<"\n"; exit(-1) ;}
     if (si>size[i]) {
       delete []X[i] ;
       X[i]=new double[si];
     }
   } else {
     X[i]=new double[si];
   }
     
   size[i]=si ;
   nel *=si   ;

 }

 // coordinates
 for (i=0; i<max(DIM,3); i++) {
   fscanf(fid,"%s %s %s",s,s,fd) ; // X_COORDINATES xx float / double
   fscanf(fid,"%c",&c)     ; // EOL

   if (i<DIM) { Ext->XA[i]=HUGE_VAL, Ext->XE[i]=-HUGE_VAL ;}
   for (j=0; j<size[i]; j++) { // read all values
     if (ascii) {
       fscanf(fid,"%le",&x) ;
     } else {
       if (strcmp(fd,"float")==0) {
         fread(&fx,sizeof(float),1,fid) ;
         if (do_swap) UDF_swap(&fx,sizeof(float),1) ;
         x=(double)fx ;
       } else {
         fread(&x,sizeof(double),1,fid) ;
         if (do_swap) UDF_swap(&x,sizeof(double),1) ;
       }
     } 
     if (i<DIM) {
       X[i][j]=x ;
       Ext->XA[i]=min(Ext->XA[i],x) ;
       Ext->XE[i]=max(Ext->XE[i],x) ;
     }
   } //

 }
}

//
// just read the data and extract some information about the scalars
//
template<int DIM>
void UDFInformation<DIM>::Scan(FILE *fid) {
 int i,n ;
 char s[1000], fd[100],c ;
 float fx ;
 double x ;

 if (lastfile==fid) scan=true ;

  // Read Data
 fscanf(fid,"%s %s",s,s)       ; // POINT_DATA xx
 fscanf(fid,"%s %s %s",s,s,fd) ; // SCALARS bla
 fscanf(fid,"%s %s",s,s)       ; // LOOKUP_TABLE default
 fscanf(fid,"%c",&c)           ; // EOL 
  
 mi= HUGE_VAL ;
 ma=-HUGE_VAL ;
 for (i=0; i<nel; i++) {

   if (ascii) {
     fscanf(fid,"%le",&x) ;
    } else {
     if (strcmp(fd,"float")==0) { 
       fread(&fx,sizeof(float),1,fid) ;
       if (do_swap) UDF_swap(&fx,sizeof(float),1) ; 
       x=(double)fx ;
     } else {
       fread(&x,sizeof(double),1,fid) ;
       if (do_swap) UDF_swap(&x,sizeof(double),1) ;
     }
    } // binary

   mi=min(mi,x) ;
   ma=max(ma,x) ;
 } // i
 fclose(fid) ;
}

#endif

---