#ifndef MICROFIBER_H#define MICROFIBER_H#include "iotools.h"#include "matrix.h"#include "vector.h"class microfiber{ public:  microfiber (void);  ~microfiber (void);  void read (XFILE *in);  void matstiff (matrix &d,long ipp,long ido);  void nlstresses (long ipp,long ido);  void updateval(long ipp,long im, long ido);  //material parameters  long numberOfMicroplanes;  double e,nu; protected:  vector  microplaneWeights; //(numberOfMicroplanes);  matrix projN,projM,projL;  vector kronecker; //(6);  double k1,k2,k3,k4,c3,c20,c1,c2,c4,c5,c6,c7,c8,c9,c10,c11,c12,c13,c14,c15,c16,c17,c18,c19,mu,ev,ed,et;  double k10,k11,k12,k13;  void initializeData (long numberOfMicroplanes );  inline double macbra(double x);  inline double FVplus(double epsV);  inline double FVminus (double epsV);  inline double FDminus(double epsD);  inline double FDplus(double epsD);  inline double FN(double epsN,double sigmaV);  inline double FT (double sigmaN,double epsV);  inline double FFiber(double epsnN);  inline double maxim (double a,double b);  inline double minim (double a,double b);};#endif