src/materials.cpp

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <assert.h>
#include "myrlog.h"

#include "materials.h"

using namespace colarray;

Materials::Materials(const char* filename) :
  _filename(filename), _numOfMaterials(0)
{

  // Check if filename is empty. If so, then all materials are assumed
  // to be vacuum
  if (strcmp(_filename.c_str(), "") == 0){
    rInfo("Assuming no materials (vacuum).");
    _numOfMaterials = 0;
    
  } else {
    rInfo("Reading material file %s...", filename);
    ifstream matFile(_filename.c_str()); 
    assert(matFile != 0);
    
    matFile >> _numOfMaterials;
    _mu.resize(_numOfMaterials + 1);     // COUNT STARTS FROME ONE!!!
    _eps.resize(_numOfMaterials + 1);

    for(int i=1; i<_numOfMaterials+1; i++){
      matFile >> _mu[i];
      matFile >> _eps[i];
    }
    rInfo("%d materials read.", _numOfMaterials);
    
  }
}




Materials::~Materials()
{
}




void Materials::getParameters(int id, Matrix<double>& xi, Matrix<double>& mu, Matrix<double>& eps)
{
  // Initialise to zero ...
  mu.fill(0.0);
  eps.fill(0.0);

  // If numOfMaterials is zero, no materials (all vacuum) are assumed
  if (_numOfMaterials == 0){
    mu(0,0) = mu(1,1) = mu(2,2) = 1.0;
    eps(0,0) = eps(1,1) = eps(2,2) = 1.0;
  } else {
    mu(0,0) = mu(1,1) = mu(2,2) = _mu[id];
    eps(0,0) = eps(1,1) = eps(2,2) = _eps[id];
  }
}

---