Heed::ulsvolume Class Reference

Unlimited surfaces volume. More...

#include <surface.h>

Inheritance diagram for Heed::ulsvolume:

Public Member Functions
	ulsvolume ()
	Default constructor.
	ulsvolume (surface *fsurf[pqqsurf], int fqsurf, char *fname, vfloat fprec)
	ulsvolume (ulsvolume &f)
	ulsvolume (const ulsvolume &fv)
virtual	~ulsvolume ()
	Destructor.
int	check_point_inside (const point &fpt, const vec &dir) const
int	range_ext (trajestep &fts, int s_ext) const
void	ulsvolume_init (surface *fsurf[pqqsurf], int fqsurf, const std::string &fname, vfloat fprec)
virtual void	income (gparticle *)
virtual void	chname (char *nm) const
virtual void	print (std::ostream &file, int l) const
virtual ulsvolume *	copy () const
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, PassivePtr< 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 Member Functions inherited from Heed::RegPassivePtr
	RegPassivePtr (void)
	RegPassivePtr (char fs_ban_del, char fs_ban_sub, char fs_ban_cop=0)
	RegPassivePtr (const RegPassivePtr &f)
RegPassivePtr &	operator= (const RegPassivePtr &f)
CountPP_ns::CountPassivePtr *	book (void) const
void	clear_pointers (void) const
virtual RegPassivePtr *	copy () const
virtual	~RegPassivePtr ()
virtual void	print (std::ostream &file, int l=1) const
void	set_s_ban_del (char fs_ban_del)
char	get_s_ban_del (void) const
void	set_s_ban_sub (char fs_ban_sub)
char	get_s_ban_sub (void) const
void	set_s_ban_cop (char fs_ban_cop)
char	get_s_ban_cop (void) const
void	set_s_allow_del_at_zero_count (char fs_allow_del_at_zero_count)
char	get_s_allow_del_at_zero_count (void) const
long	get_total_number_of_references (void) const

Public Attributes
int	qsurf
ActivePtr< surface >	surf [pqqsurf]
std::string	name
Public Attributes inherited from Heed::absvol
vfloat	prec
bool	s_sensitive

Protected Member Functions
virtual void	get_components (ActivePtr< absref_transmit > &aref_tran)

Protected Attributes
surface *	adrsurf [pqqsurf]

Additional Inherited Members
Static Public Member Functions inherited from Heed::RegPassivePtr
static void	set_s_ban_del_ignore (char fs_ban_del_ignore)
static char	get_s_ban_del_ignore (void)
static void	set_s_print_adr_cpp (char fs_print_adr_cpp)
static char	get_s_print_adr_cpp (void)

Detailed Description

Unlimited surfaces volume.

Definition at line 124 of file surface.h.

Constructor & Destructor Documentation

ulsvolume() [1/4]

Heed::ulsvolume::ulsvolume

(

void

)

Default constructor.

Definition at line 335 of file surface.cpp.

                         : qsurf(0) {
  name = std::string("not inited ulsvolume");
}

Referenced by copy().

ulsvolume() [2/4]

Heed::ulsvolume::ulsvolume	(	surface *	fsurf[pqqsurf],
		int	fqsurf,
		char *	fname,
		vfloat	fprec
	)

Definition at line 352 of file surface.cpp.

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

ulsvolume() [3/4]

Heed::ulsvolume::ulsvolume

(

ulsvolume &

f

)

Definition at line 361 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].put(f.surf[n].get());
}

ulsvolume() [4/4]

Heed::ulsvolume::ulsvolume

(

const ulsvolume &

fv

)

Definition at line 369 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].put(f.surf[n].get());
}

~ulsvolume()

virtual Heed::ulsvolume::~ulsvolume ( )

inlinevirtual

Destructor.

Definition at line 159 of file surface.h.

{}

Member Function Documentation

check_point_inside()

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

virtual

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 194 of file surface.cpp.

                                                                        {
  mfunname("ulsvolume::check_point_inside(...)");
  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()

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

inlinevirtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 170 of file surface.h.

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

Referenced by print().

copy()

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

inlinevirtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 175 of file surface.h.

{ return new ulsvolume(*this); }

get_components()

void Heed::ulsvolume::get_components ( ActivePtr< absref_transmit > &  aref_tran )

protectedvirtual

Reimplemented from Heed::absref.

Definition at line 189 of file surface.cpp.

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

income()

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

inlinevirtual

Reimplemented from Heed::absvol.

Definition at line 169 of file surface.h.

{}

print()

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

virtual

Reimplemented from Heed::absvol.

Reimplemented in Heed::manip_ulsvolume.

Definition at line 377 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::box::print(), Heed::manip_ulsvolume::print(), and range_ext().

range_ext()

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

virtual

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 212 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
  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	(	surface *	fsurf[pqqsurf],
		int	fqsurf,
		const std::string &	fname,
		vfloat	fprec
	)

Definition at line 339 of file surface.cpp.

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

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

Member Data Documentation

adrsurf

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

protected

Definition at line 149 of file surface.h.

Referenced by get_components().

name

std::string Heed::ulsvolume::name

Definition at line 146 of file surface.h.

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

qsurf

int Heed::ulsvolume::qsurf

Definition at line 144 of file surface.h.

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

surf

ActivePtr<surface> Heed::ulsvolume::surf[pqqsurf]

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