Mille.h

#ifndef MILLE_H
#define MILLE_H

#include <fstream>

class Mille 
{
 public:
  Mille(const char *outFileName, bool asBinary = true, bool writeZero = false);
  ~Mille();

  void mille(int NLC, const float *derLc, int NGL, const float *derGl,
             const int *label, float rMeas, float sigma);
  void special(int nSpecial, const float *floatings, const int *integers);
  void kill();
  void end();

 private:
  void newSet();
  bool checkBufferSize(int nLocal, int nGlobal);

  std::ofstream myOutFile; // C-binary for output
  bool myAsBinary;         // if false output as text
  bool myWriteZero;        // if true also write out derivatives/lables ==0

  enum {myBufferSize = 5000};
  int   myBufferInt[myBufferSize];   // to collect labels etc.
  float myBufferFloat[myBufferSize]; // to collect derivatives etc.
  int   myBufferPos;
  bool  myHasSpecial; // if true, special(..) already called for this record

  enum {myMaxLabel = (0xFFFFFFFF - (1 << 31))}; // largest label allowed: 2^31 - 1
};
#endif