DumperVGroupe< T, Dim > Class Template Reference

#include <dumper_vgroupe.h>

Inheritance diagram for DumperVGroupe< T, Dim >:

[legend]Collaboration diagram for DumperVGroupe< T, Dim >:

[legend]List of all members.

Public Member Functions
	DumperVGroupe (DomainInterface &d, int ID)
Protected Member Functions
void	Dump ()
	dump function that dump to file screen or whatever
void	Init ()
	initialisation of the dumper internal data
void	SetParam (char *key, double value)
	setting parameters to object that auto configure fields from key/value call
Private Attributes
double *	deplacements_par_tranche
unsigned int *	nb_par_tranche
unsigned int	nb_tranches
FILE *	output [1]

---

Detailed Description

template<typename T, unsigned int Dim>
class DumperVGroupe< T, Dim >

Definition at line 52 of file dumper_vgroupe.h.

---

Constructor & Destructor Documentation

template<typename T, unsigned int Dim>

DumperVGroupe< T, Dim >::DumperVGroupe	(	DomainInterface &	d,
		int	ID
	)			[inline]

Definition at line 56 of file dumper_vgroupe.h.

                                           :
    Dumper<T,Dim>(static_cast<T &>(d),ID)
    {
    };

---

Member Function Documentation

template<typename Domain, unsigned int Dim>

void DumperVGroupe< Domain, Dim >::Dump

(

)

[inline, protected, virtual]

dump function that dump to file screen or whatever

Implements DumperInterface.

Definition at line 71 of file dumper_atom_vgroup.cpp.

References UnitsConverter::AtomsUnits, cpt, DBG_ALL, DumperVGroupe< T, Dim >::deplacements_par_tranche, DUMP, UnitsConverter::etalon_metre, ipow(), DumperVGroupe< T, Dim >::nb_par_tranche, DumperVGroupe< T, Dim >::nb_tranches, DumperVGroupe< T, Dim >::output, and UnitsConverter::SetReadUnits().

                                    {
  UnitsConverter unit;
  unit.SetReadUnits(UnitsConverter::AtomsUnits);


  //  Boule & b = static_cast<Boule &>(*this->geom);
  unsigned int couche_index=0;
  unsigned int cpt = 0;

  typedef typename Domain::RefDof RefAtom;
  typedef typename Domain::ContenerDofs ContenerAtoms;
  typedef typename Domain::IteratorDofs IteratorAtoms;

  ContenerAtoms & c = Dumper<Domain,Dim>::dom.getContener();
  IteratorAtoms & it = c.GetIterator();
  RefAtom at = it.GetFirst();

  memset(nb_par_tranche,0,sizeof(int)*nb_tranches);
  memset(deplacements_par_tranche,0,sizeof(double)*nb_tranches);
  
  at = it.GetFirst();
  while (!it.end())
    {
      double d;
      if (Dim == 1)
        d = fabs(at.position0(0));
      else if (Dim == 2)
        d = ipow(at.position0(0),2)+ipow(at.position0(1),2);
      else if (Dim == 3)
        d = ipow(at.position0(0),2)+ipow(at.position0(1),2)+ipow(at.position0(2),2);

      /* la distance au centre de la boule est utilise pour caracteriser la projection dans l'espace des phases */

      DUMP("le d trouve = " << d,DBG_ALL);
      //      d = (d - b.Rmin())/(b.Rmax()-b.Rmin());

      DUMP("le d rapporte dans 0..1 = " << d,DBG_ALL);
      /* petit calcul simpliste pour trouver la tranche d'oignon que l'ont veux */
      couche_index = static_cast<unsigned int>(nb_tranches*d);
      if (Dim == 1) couche_index = cpt;
      DUMP("couche_index trouve = " << couche_index,DBG_ALL);

      /* prepare le deplacement */

      double ux = (at.position(0) - at.position0(0))*unit.etalon_metre;
      double uy = (at.position(1) - at.position0(1))*unit.etalon_metre;
      double uz = (at.position(2) - at.position0(2))*unit.etalon_metre;

      double dep = sqrt(ux*ux+uy*uy+uz*uz);
      if (Dim == 1)
        dep = ux;
      DUMP("deplacement calcule = " << dep,DBG_ALL);
      /* on additionne le deplacement de chaque atome dans la couche qui le concerne */
      deplacements_par_tranche[couche_index] += dep;
      /* on augmente le nombre d'atomes dans la couche - il faudra faire la moyenne a la fin !*/
      ++nb_par_tranche[couche_index];
      
      DUMP("deplacement_par_tranche[" << couche_index << "]= " << deplacements_par_tranche[couche_index],DBG_ALL);
      DUMP("nb_par_tranche[" << couche_index << "]= " << nb_par_tranche[couche_index],DBG_ALL);
      at = it.GetNext();
      ++cpt;
    }
  
  delete &it;
  
  for (unsigned int k=0;k < nb_tranches;++k)
    {
      /* on dump directement la moyenne */
      if (nb_par_tranche[couche_index] == 0) continue;
      DUMP("deplacement_moy_par_tranche[" << couche_index << "]= " << deplacements_par_tranche[k] << "/" << nb_par_tranche[couche_index],DBG_ALL);
      fprintf(output[0],"%d\t%d\t%.15e\n",Dumper<Domain,Dim>::dump_step,k+1,deplacements_par_tranche[k]/nb_par_tranche[couche_index]);
    }    
}

Here is the call graph for this function:

template<typename Domain, unsigned int Dim>

void DumperVGroupe< Domain, Dim >::Init

(

)

[inline, protected, virtual]

initialisation of the dumper internal data

Implements DumperInterface.

Definition at line 55 of file dumper_atom_vgroup.cpp.

References DumperVGroupe< T, Dim >::deplacements_par_tranche, DumperVGroupe< T, Dim >::nb_par_tranche, NB_TRANCHES, DumperVGroupe< T, Dim >::nb_tranches, and DumperVGroupe< T, Dim >::output.

                                    {
  /* j'initialise les tranches pour la projection */
  nb_tranches = NB_TRANCHES;
  if (Dim == 1)
    nb_tranches = Dumper<Domain,Dim>::dom.getContener().nbElem();
  
  deplacements_par_tranche = new double[nb_tranches];
  nb_par_tranche = new unsigned int[nb_tranches];
  
  /* je creer le fichier pour le dump de la vitesse de groupe */
  char temp[255];
  snprintf(temp,255,"%s_vgroupe.plot",Dumper<Domain,Dim>::base_name);
  output[0] = fopen(temp,"wb+");    
}

template<typename Domain, unsigned int Dim>

void DumperVGroupe< Domain, Dim >::SetParam	(	char *	key,
		double	value
	)			[inline, protected, virtual]

setting parameters to object that auto configure fields from key/value call

Implements DumperInterface.

Definition at line 50 of file dumper_atom_vgroup.cpp.

                                                              {

}

---

Member Data Documentation

template<typename T, unsigned int Dim>

double* DumperVGroupe< T, Dim >::deplacements_par_tranche [private]

Definition at line 70 of file dumper_vgroupe.h.

Referenced by DumperVGroupe< T, Dim >::Dump(), and DumperVGroupe< T, Dim >::Init().

template<typename T, unsigned int Dim>

unsigned int* DumperVGroupe< T, Dim >::nb_par_tranche [private]

Definition at line 71 of file dumper_vgroupe.h.

Referenced by DumperVGroupe< T, Dim >::Dump(), and DumperVGroupe< T, Dim >::Init().

template<typename T, unsigned int Dim>

unsigned int DumperVGroupe< T, Dim >::nb_tranches [private]

Definition at line 69 of file dumper_vgroupe.h.

Referenced by DumperVGroupe< T, Dim >::Dump(), and DumperVGroupe< T, Dim >::Init().

template<typename T, unsigned int Dim>

FILE* DumperVGroupe< T, Dim >::output[1] [private]

Definition at line 73 of file dumper_vgroupe.h.

Referenced by DumperVGroupe< T, Dim >::Dump(), and DumperVGroupe< T, Dim >::Init().

---

The documentation for this class was generated from the following files:

• dumper/dumper_vgroupe.h
• dumper/dumper_atom_vgroup.cpp

---