Material.h

Go to the documentation of this file.

#ifndef MATERIAL_H_
#define MATERIAL_H_

#ifdef __cplusplus
#ifdef SWIG
#include "Python.h"
#endif
#include "constants.h"
#include "log.h"
#include "linalg.h"
#include <sstream>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <malloc.h>
#endif

#ifdef INTEL

#define MM_MALLOC(size,alignment) _mm_malloc(size, alignment)

#define MM_FREE(array) _mm_free(array)

#else

#define MM_MALLOC(size,alignment) memalign(alignment, size)

#define MM_FREE(array) free(array)
#endif


int material_id();
void reset_material_id();
void maximize_material_id(int material_id);


class Material {

private:

  int _id;

  char* _name;

  double _volume;

  int _num_instances;

  int _num_groups;

  FP_PRECISION* _sigma_t;

  FP_PRECISION _max_sigma_t;

  FP_PRECISION* _sigma_s;

  FP_PRECISION* _sigma_a;

  FP_PRECISION* _sigma_f;

  FP_PRECISION* _nu_sigma_f;

  FP_PRECISION* _chi;

  FP_PRECISION* _fiss_matrix;

  bool _fissionable;

  bool _data_aligned;

  int _num_vector_groups;

public:
  Material(int id=0, const char* name="");
  virtual ~Material();

  int getId() const;
  char* getName() const;
  double getVolume();
  int getNumInstances();
  int getNumEnergyGroups() const;
  FP_PRECISION* getSigmaT();
  FP_PRECISION getMaxSigmaT();
  FP_PRECISION* getSigmaS();
  FP_PRECISION* getSigmaA();
  FP_PRECISION* getSigmaF();
  FP_PRECISION* getNuSigmaF();
  FP_PRECISION* getChi();
  FP_PRECISION* getFissionMatrix();
  FP_PRECISION getSigmaTByGroup(int group);
  FP_PRECISION getSigmaSByGroup(int origin, int destination);
  FP_PRECISION getSigmaAByGroup(int group);
  FP_PRECISION getSigmaFByGroup(int group);
  FP_PRECISION getNuSigmaFByGroup(int group);
  FP_PRECISION getChiByGroup(int group);
  FP_PRECISION getFissionMatrixByGroup(int origin, int destination);
  bool isFissionable();
  bool isDataAligned();
  int getNumVectorGroups();

  void setName(const char* name);
  void setVolume(double volume);
  void incrementVolume(double volume);
  void setNumInstances(int num_instances);
  void incrementNumInstances();
  void setNumEnergyGroups(const int num_groups);

  void setSigmaT(double* xs, int num_groups);
  void setSigmaS(double* xs, int num_groups);
  void setSigmaF(double* xs, int num_groups);
  void setNuSigmaF(double* xs, int num_groups);
  void setChi(double* xs, int num_groups);

  void setSigmaTByGroup(double xs, int group);
  void setSigmaFByGroup(double xs, int group);
  void setNuSigmaFByGroup(double xs, int group);
  void setSigmaSByGroup(double xs, int origin, int destination);
  void setChiByGroup(double xs, int group);
  void setSigmaAByGroup(double xs, int group);

  void buildFissionMatrix();
  void transposeProductionMatrices();
  void alignData();
  Material* clone();

  std::string toString();
  void printString();
};


#endif /* MATERIAL_H_ */