Heed::ulsvolume Class Reference

Unlimited surfaces volume. More...

#include <surface.h>

Inheritance diagram for Heed::ulsvolume:

Public Member Functions
	ulsvolume ()
	Default constructor.
	ulsvolume (const std::vector< std::shared_ptr< surface > > &fsurf, char *fname, vfloat fprec)
	Constructor from surfaces.
	ulsvolume (ulsvolume &f)
	ulsvolume (const ulsvolume &fv)
virtual	~ulsvolume ()
	Destructor.
int	check_point_inside (const point &fpt, const vec &dir) const override
int	range_ext (trajestep &fts, int s_ext) const override
void	ulsvolume_init (const std::vector< std::shared_ptr< surface > > &fsurf, const std::string &fname, vfloat fprec)
void	income (gparticle *) override
void	chname (char *nm) const override
void	print (std::ostream &file, int l) const override
ulsvolume *	copy () const override
Public Member Functions inherited from Heed::absvol
virtual	~absvol ()
	Destructor.
virtual int	check_point_inside (const point &fpt, const vec &dir) const =0
virtual int	find_embed_vol (const point &fpt, const vec &dir, manip_absvol_treeid *atid) const
virtual int	range (trajestep &fts, int s_ext, int &sb, manip_absvol *&faeid) const
virtual int	range_ext (trajestep &fts, int s_ext) const =0
virtual void	income (gparticle *)
virtual void	chname (char *nm) const
virtual void	print (std::ostream &file, int l) const
virtual absvol *	copy () const
virtual std::vector< manip_absvol * >	Gamanip_embed () const
Public Member Functions inherited from Heed::absref
virtual	~absref ()
	Destructor.
virtual void	down (const abssyscoor *fasc)
	Convert numbering representation of object to basical system of fasc.
virtual void	up (const abssyscoor *fasc)
	Convert numbering representation of objects to new system.
virtual void	turn (const vec &dir, vfloat angle)
	Turn around axis doing via center of coordinate system along dir.
virtual void	shift (const vec &dir)

Public Attributes
int	qsurf = 0
std::array< std::shared_ptr< surface >, pqqsurf >	surf
std::string	name = "non-initialized ulsvolume"
Public Attributes inherited from Heed::absvol
vfloat	prec
bool	s_sensitive

Static Public Attributes
static constexpr int	pqqsurf = 10

Protected Member Functions
absref_transmit	get_components () override

Protected Attributes
surface *	adrsurf [pqqsurf]

Detailed Description

Unlimited surfaces volume.

Definition at line 123 of file surface.h.

Constructor & Destructor Documentation

ulsvolume() [1/4]

Heed::ulsvolume::ulsvolume ( )

inline

Default constructor.

Definition at line 156 of file surface.h.

{}

Referenced by copy().

ulsvolume() [2/4]

Heed::ulsvolume::ulsvolume	(	const std::vector< std::shared_ptr< surface > > &	fsurf,
		char *	fname,
		vfloat	fprec
	)

Constructor from surfaces.

Definition at line 347 of file surface.cpp.

    : qsurf(fsurf.size()), name(fname) {
  mfunname("ulsvolume::ulsvolume(...)");
  check_econd12(qsurf, >, pqqsurf, mcerr);
  prec = fprec;
  for (int n = 0; n < qsurf; ++n) surf[n] = fsurf[n];
}

ulsvolume() [3/4]

Heed::ulsvolume::ulsvolume

(

ulsvolume &

f

)

Definition at line 356 of file surface.cpp.

    : absref(f), absvol(f), qsurf(f.qsurf), name(f.name) {
  mfunname("ulsvolume::ulsvolume(...)");
  check_econd12(f.qsurf, >, pqqsurf, mcerr);
  prec = f.prec;
  for (int n = 0; n < qsurf; ++n) surf[n] = f.surf[n];
}

ulsvolume() [4/4]

Heed::ulsvolume::ulsvolume

(

const ulsvolume &

fv

)

Definition at line 364 of file surface.cpp.

    : absref(f), absvol(f), qsurf(f.qsurf), name(f.name) {
  mfunname("ulsvolume::ulsvolume(...)");
  check_econd12(f.qsurf, >, pqqsurf, mcerr);
  prec = f.prec;
  for (int n = 0; n < qsurf; ++n) surf[n] = f.surf[n];
}

~ulsvolume()

virtual Heed::ulsvolume::~ulsvolume ( )

inlinevirtual

Destructor.

Definition at line 162 of file surface.h.

{}

Member Function Documentation

check_point_inside()

int Heed::ulsvolume::check_point_inside ( const point &  fpt, const vec &  dir  ) const

overridevirtual

Check if a point is inside the volume. If two volumes are adjacent, it may happen that a point belongs to both. To avoid this confusion the parameter dir is used. If dir == (0, 0, 0), and the point is exactly on the border, the behaviour is in general not specified. If dir != (0, 0, 0), and the point is on the border with precision prec, the exiting volume is ignored.

Implements Heed::absvol.

Definition at line 192 of file surface.cpp.

                                                                        {
  mfunname("ulsvolume::check_point_inside(const point&, const vec&)");
  check_econd11(qsurf, <= 0, mcerr);
  for (int n = 0; n < qsurf; n++) {
    if (!(surf[n].get()->check_point_inside(fpt, dir, prec))) {
      return 0;
    }
  }
#ifdef TRACE_find_embed_vol
  indn.n++;
  Imcout << "ulsvolume::check_point_inside: the point is in volume\n";
  Imcout << "point:" << fpt;
  print(mcout, 0);
  indn.n--;
#endif
  return 1;
}

Referenced by check_point_inside(), and range_ext().

chname()

void Heed::ulsvolume::chname ( char *  nm ) const

inlineoverridevirtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 174 of file surface.h.

                                       {
    strcpy(nm, "ulsvolume: ");
    strcat(nm, name.c_str());
  }

Referenced by print().

copy()

ulsvolume * Heed::ulsvolume::copy ( ) const

inlineoverridevirtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 179 of file surface.h.

{ return new ulsvolume(*this); }

get_components()

absref_transmit Heed::ulsvolume::get_components ( )

overrideprotectedvirtual

Reimplemented from Heed::absref.

Definition at line 187 of file surface.cpp.

                                          {
  for (int n = 0; n < qsurf; n++) adrsurf[n] = surf[n].get();
  return absref_transmit(qsurf, (absref**)adrsurf);
}

income()

void Heed::ulsvolume::income ( gparticle *  )

inlineoverridevirtual

Reimplemented from Heed::absvol.

Definition at line 173 of file surface.h.

{}

print()

void Heed::ulsvolume::print ( std::ostream &  file, int  l  ) const

overridevirtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 372 of file surface.cpp.

                                                   {
  char s[1000];
  chname(s);
  Ifile << "ulsvolume::print(l=" << l << "): " << s << '\n';
  if (l > 0) {
    indn.n += 2;
    Ifile << "qsurf=" << qsurf << " prec=" << prec << '\n';
    for (int n = 0; n < qsurf; ++n) {
      Ifile << " nsurf=" << n << '\n';
      surf[n].get()->print(file, l);
    }
    absvol::print(file, l);
    indn.n -= 2;
  }
}

Referenced by check_point_inside(), Heed::manip_ulsvolume::print(), Heed::box::print(), and range_ext().

range_ext()

int Heed::ulsvolume::range_ext ( trajestep &  fts, int  s_ext  ) const

overridevirtual

Find cross with current volume ignoring embraced ones. s_ext=1 exit, now point is inside, but embraced volumes are ingnored. s_ext=0 enter, now point is outside

Implements Heed::absvol.

Definition at line 210 of file surface.cpp.

                                                        {
  mfunnamep("int ulsvolume::range_ext(trajestep& fts, int s_ext) const");
  check_econd11(qsurf, <= 0, mcerr);
#ifdef DEBUG_ulsvolume_range_ext
  mcout << "ulsvolume::range_ext, START, s_ext=" << s_ext << " qsurf=" << qsurf
        << '\n';
  mcout << fts;
#endif
  constexpr int pqcrossurf = 4;
  vfloat crange[pqcrossurf];
  point cpt[pqcrossurf];
  int fs_ext[pqcrossurf];
  int n, m, nc;
  int s = 0;  // sign of crossing
  if (s_ext == 1) {
    for (n = 0; n < qsurf; n++) {
      int qc = surf[n].get()->range(fts, crange, cpt, fs_ext);
      for (m = 0; m < qc; m++) {
        if (fs_ext[m] == 1) {
          s = 1;
          // The last minute change, it was 0 somewhy instead of m
          fts.mrange = crange[m];  // reduce the range
          fts.mpoint = cpt[m];
          break;  // take only the first exit point, it should be closest
        } else if (fs_ext[m] == 0) {
          if (!(surf[n].get()->check_point_inside(fts.currpos, fts.dir,
                                                  prec))) {
            funnw.ehdr(mcerr);
            mcerr << "\nshould never happen\n"
                  << "It may happen if you  call this function with s_ext==1\n"
                  << "for point outside the volume\n";
            spexit(mcerr);
          }
        } else if (fs_ext[m] == 2)
          break;  // don't know what to do, safe to ignore
      }
    }
 
    if (s == 1) {
      fts.s_prec = 0;
    }
    return s;
  } else {       // for if(s_ext==1)
    int ss = 0;  // sign that there is cross with any of the surfaces
    for (n = 0; n < qsurf; n++) {
#ifdef DEBUG_ulsvolume_range_ext
      Iprintn(mcout, n);
#endif
      int qc = surf[n].get()->range(fts, crange, cpt, fs_ext);
#ifdef DEBUG_ulsvolume_range_ext
      mcout << "ulsvolume::range_ext: qc=" << qc << "\n";
      surf[n]->print(mcout, 1);
#endif
      for (nc = 0; nc < qc; nc++)  // loop by crossing points
      {
#ifdef DEBUG_ulsvolume_range_ext
        mcout << "nc=" << nc << " fs_ext[nc]=" << fs_ext[nc] << '\n';
#endif
        if (fs_ext[nc] == 0)  // thus ignoring exitted surfaces
        {
          s = 1;
          for (m = 0; m < qsurf; m++)  // scan other surfaces and verify that
          {                            // the crossing point is inside
            if (m != n) {
              if (surf[m].get()->check_point_inside1(cpt[nc], fs_ext[nc],
                                                     prec) == 0) {
#ifdef DEBUG_ulsvolume_range_ext
                mcout << "m=" << m << '\n';
                mcout << "Since the point is outside of the other surface, "
                      << "it can not be border of volume\n";
#endif
                s = 0;
                break;
              }
            }
          }
#ifdef DEBUG_ulsvolume_range_ext
          Iprintn(mcout, s);
#endif
          if (s == 1) {
#ifdef DEBUG_ulsvolume_range_ext
            mcout << "The crossing point is inside all other surfaces, \n"
                  << "so it is good crossing point\n";
#endif
            ss = 1;
            fts.mrange = crange[nc];
            fts.mpoint = cpt[nc];
            break;  // since points are ordered, go to next surface,
                    // may be there is nearer crossing point
          }
        }
      }
    }
    if (ss == 1) {
      fts.s_prec = 0;
    }
#ifdef DEBUG_ulsvolume_range_ext
    mcout << "ulsvolume::range_ext: at the end\n";
    print(mcout, 1);
    mcout << "ss=" << ss << '\n';
#endif
    return ss;
  }
}

Referenced by Heed::box::range_ext().

ulsvolume_init()

void Heed::ulsvolume::ulsvolume_init	(	const std::vector< std::shared_ptr< surface > > &	fsurf,
		const std::string &	fname,
		vfloat	fprec
	)

Definition at line 334 of file surface.cpp.

                                                                     {
  prec = fprec;
  name = fname;
  if (qsurf > 0) {
    for (int n = 0; n < qsurf; ++n) surf[n].reset();
  }
  qsurf = fsurf.size();
  for (int n = 0; n < qsurf; ++n) {
    surf[n] = fsurf[n];
  }
}

Referenced by Heed::box::init_planes().

Member Data Documentation

adrsurf

surface* Heed::ulsvolume::adrsurf[pqqsurf]

protected

Definition at line 150 of file surface.h.

Referenced by get_components().

name

std::string Heed::ulsvolume::name = "non-initialized ulsvolume"

Definition at line 146 of file surface.h.

Referenced by Heed::manip_ulsvolume::chname(), chname(), and ulsvolume_init().

pqqsurf

constexpr int Heed::ulsvolume::pqqsurf = 10

staticconstexpr

Definition at line 143 of file surface.h.

Referenced by ulsvolume().

qsurf

int Heed::ulsvolume::qsurf = 0

Definition at line 144 of file surface.h.

Referenced by check_point_inside(), get_components(), print(), range_ext(), ulsvolume(), and ulsvolume_init().

surf

std::array<std::shared_ptr<surface>, pqqsurf> Heed::ulsvolume::surf

Definition at line 145 of file surface.h.

Referenced by check_point_inside(), get_components(), print(), range_ext(), ulsvolume(), and ulsvolume_init().

---

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

• surface.h
• surface.cpp