vec.cpp File Reference

#include <stdlib.h>
#include <iomanip>
#include <cmath>
#include "wcpplib/geometry/vec.h"
#include "wcpplib/util/emul_new_stand.h"
#include <CLHEP/Random/RandomEngine.h>
#include "geometry/vec.ic"

Go to the source code of this file.

Functions
	macro_copy_body (absref_transmit) void absref_transmit
vfloat	cos2vec (const vec &r1, const vec &r2)
vfloat	ang2vec (const vec &r1, const vec &r2)
vfloat	sin2vec (const vec &r1, const vec &r2)
vec	project_to_plane (const vec &r, const vec &normal)
vfloat	ang2projvec (const vec &r1, const vec &r2, const vec &normal)
std::ostream &	operator<< (std::ostream &file, const vec &v)
std::ostream &	operator<< (std::ostream &file, const vecReg &v)
std::ostream &	operator<< (std::ostream &file, const basis &b)
std::ostream &	operator<< (std::ostream &file, const basisReg &b)
std::ostream &	operator<< (std::ostream &file, const point &p)
std::ostream &	operator<< (std::ostream &file, const pointReg &b)
std::ostream &	operator<< (std::ostream &file, const abssyscoor &f)
std::ostream &	operator<< (std::ostream &file, const fixsyscoor &f)

Variables
HepRandomEngine &	random_engine
int	vecerror = 0
vec	dex (1, 0, 0)
vec	dey (0, 1, 0)
vec	dez (0, 0, 1)
vec	dv0 (0, 0, 0)

Function Documentation

ang2projvec()

vfloat ang2projvec	(	const vec &	r1,
		const vec &	r2,
		const vec &	normal
	)

Definition at line 212 of file vec.cpp.

                                                                    {
  pvecerror(
      "vfloat ang2projvec(const vec& r1, const vec& r2, const vec& normal)");
  vec rt1 = project_to_plane(r1, normal);
  vec rt2 = project_to_plane(r2, normal);
  if (rt1 == dv0 || rt2 == dv0) {
    vecerror = 1;
    return 0;
  }
  vfloat tang = ang2vec(rt1, rt2);
  if (tang == 0) return tang;  // projections are parallel
  vec at = rt1 || rt2;
  int i = check_par(at, normal, 0.0001);
  //mcout<<"r1="<<r1<<"r2="<<r2<<"normal="<<normal
  //     <<"rt1="<<rt1<<"rt2="<<rt2<<"\ntang="<<tang
  //     <<"\nat="<<at<<" i="<<i<<'\n';
  if (i == -1) return 2.0 * M_PI - tang;
  return tang;  // it works if angle <= PI
}

Referenced by Heed::polygon::check_point_in(), and Heed::splane::range().

ang2vec()

vfloat ang2vec	(	const vec &	r1,
		const vec &	r2
	)

Definition at line 165 of file vec.cpp.

                                             {
  // angle between vectors
  // instead of return acos(cos2vec(r1,r2)); which produces NaN on linux at
  // parallel vectors
  vfloat cs = cos2vec(r1, r2);
  if (vecerror != 0) return 0;
  if (cs > 0.707106781187 || cs < -0.707106781187) {  // 1.0/sqrt(2)
    // pass to sin, it will be more exactly
    vfloat sn = sin2vec(r1, r2);
    if (vecerror != 0) return 0;
    if (cs > 0.0)
      return asin(sn);
    else
      return M_PI - asin(sn);
  }
  return acos(cs);
}

cos2vec()

vfloat cos2vec	(	const vec &	r1,
		const vec &	r2
	)

Definition at line 142 of file vec.cpp.

                                             {
  // cosinus of angle between vectors
  // If one of vectors has zero length, it returns 2.
  pvecerror("vfloat cos2vec(const vec& r1, const vec& r2)");
  vfloat lr1 = length2(r1);
  vfloat lr2 = length2(r2);
  //mcout<<"cos2vec:\n";
  //Iprintn(mcout, lr1);
  //Iprintn(mcout, lr2);
  if (lr1 == 0 || lr2 == 0) {
    vecerror = 1;
    return 0;
  }
  vfloat cs = r1 * r2;
  int sign = 1;
  if (cs < 0) sign = -1;
  cs = cs * cs;
  cs = sign * sqrt(cs / (lr1 * lr2));
  //mcout<<"r1="<<r1<<"r2="<<r2<<"cos="<<cs<<'\n';
  return cs;
  //return r1*r2/(lr1*lr2);
}

Referenced by Heed::spquadr::apos(), basis::basis(), Heed::splane::check_point_inside(), Heed::splane::check_point_inside1(), Heed::straight::distance(), Heed::mparticle::new_speed(), and Heed::splane::range().

macro_copy_body()

macro_copy_body

(

absref_transmit

)

Definition at line 33 of file vec.cpp.

                                                         {
  if (l > 0) {
    Ifile << "absref_transmit::print(l=" << l << ") qaref=" << qaref
          << " qaref_pointer=" << qaref_pointer
          << " qaref_other=" << qaref_other << "\n";
    file.flush();
  }
}
 
absref* absref_transmit::get_other(int /*n*/) { return NULL; }

operator<<() [1/8]

std::ostream & operator<<	(	std::ostream &	file,
		const abssyscoor &	f
	)

Definition at line 595 of file vec.cpp.

                                                              {
  f.print(file, 2);
  return file;
}

operator<<() [2/8]

std::ostream & operator<<	(	std::ostream &	file,
		const basis &	b
	)

Definition at line 509 of file vec.cpp.

                                                         {
  Ifile << "basis: name=" << b.name << '\n';
  indn.n += 2;
  int indnsave = indn.n;
  Ifile << "ex: ";
  indn.n = 0;
  file << b.ex;
  indn.n = indnsave;
  Ifile << "ey: ";
  indn.n = 0;
  file << b.ey;
  indn.n = indnsave;
  Ifile << "ez: ";
  indn.n = 0;
  file << b.ez;
  indn.n = indnsave;
  indn.n -= 2;
  //file << '\n';
  return file;
}

operator<<() [3/8]

std::ostream & operator<<	(	std::ostream &	file,
		const basisReg &	b
	)

Definition at line 530 of file vec.cpp.

                                                            {
  Ifile << "basisReg=" << ((basis&)b);
  return file;
}

operator<<() [4/8]

std::ostream & operator<<	(	std::ostream &	file,
		const fixsyscoor &	f
	)

Definition at line 633 of file vec.cpp.

                                                              {
  Ifile << "fixsyscoor:\n";
  f.RegPassivePtr::print(file, 2);
  f.abssyscoor::print(file, 2);
  return file;
}

operator<<() [5/8]

std::ostream & operator<<	(	std::ostream &	file,
		const point &	p
	)

Definition at line 557 of file vec.cpp.

                                                         {
  Ifile << "point:\n";
  indn.n += 2;
  file << p.v;
  indn.n -= 2;
  return file;
}

operator<<() [6/8]

std::ostream & operator<<	(	std::ostream &	file,
		const pointReg &	b
	)

Definition at line 565 of file vec.cpp.

                                                            {
  Ifile << "pointReg=" << ((point&)b);
  return file;
}

operator<<() [7/8]

std::ostream & operator<<	(	std::ostream &	file,
		const vec &	v
	)

Definition at line 346 of file vec.cpp.

                                                       {
  Ifile << "vector=" << std::setw(13) << v.x << std::setw(13) << v.y
        << std::setw(13) << v.z;
  file << '\n';
  file.flush();
  return file;
}

operator<<() [8/8]

std::ostream & operator<<	(	std::ostream &	file,
		const vecReg &	v
	)

Definition at line 359 of file vec.cpp.

                                                          {
  Ifile << "vecReg=" << ((vec&)v);
  return file;
}

project_to_plane()

vec project_to_plane	(	const vec &	r,
		const vec &	normal
	)

Definition at line 200 of file vec.cpp.

                                                      {
  pvecerror("vec project_to_plane(const vec& r, const vec& normal)");
  vec per(normal || r);
  if (per == dv0) {
    // either one of vectors is 0 or they are parallel
    return dv0;
  }
  vec ax = unit_vec(per || normal);
  vfloat v = ax * r;
  return v * ax;
}

Referenced by Heed::mparticle::curvature(), and Heed::mparticle::new_speed().

sin2vec()

vfloat sin2vec	(	const vec &	r1,
		const vec &	r2
	)

Definition at line 183 of file vec.cpp.

                                             {
  // sinus of angle between vectors
  pvecerror("vfloat sin2vec(const vec& r1, const vec& r2)");
  vfloat lr1 = length2(r1);
  vfloat lr2 = length2(r2);
  if (lr1 == 0 || lr2 == 0) {
    vecerror = 1;
    return 0;
  }
  vfloat sn = length(r1 || r2);
  sn = sn * sn;
  sn = sqrt(sn / (lr1 * lr2));
  //mcout<<"r1="<<r1<<"r2="<<r2<<"sin="<<sn<<'\n';
  return sn;
  //return sin(ang2vec(r1,r2));
}

Referenced by Heed::straight::distance().

Variable Documentation

dex

vec dex(1, 0, 0)	(	1	,
		0	,
		0
	)

dey

vec dey(0, 1, 0)	(	0	,
		1	,
		0
	)

dez

vec dez(0, 0, 1)	(	0	,
		0	,
		1
	)

dv0

vec dv0(0, 0, 0)	(	0	,
		0	,
		0
	)

Referenced by gparticle::calc_step_to_bord(), Heed::box::check_point_inside(), Heed::splane::check_point_inside(), Heed::straight::cross(), Heed::mparticle::curvature(), gparticle::gparticle(), Heed::mparticle::mparticle(), Heed::mparticle::new_speed(), gparticle::step(), Heed::mparticle::step(), and trajestep::trajestep().

random_engine

HepRandomEngine& random_engine

extern

Referenced by vec::random_conic_vec(), vec::random_round_vec(), and vec::random_sfer_vec().

vecerror

int vecerror = 0

Definition at line 31 of file vec.cpp.

Referenced by Heed::spquadr::apos(), basis::basis(), Heed::mparticle::check_consistency(), Heed::circumf::cross(), Heed::polyline::cross(), Heed::polygon::cross(), Heed::plane::cross(), Heed::straight::cross(), Heed::mparticle::new_speed(), Heed::spquadr::perpos(), Heed::HeedDeltaElectron::physics_after_new_speed(), Heed::plane::plane(), Heed::polygon::polygon(), Heed::polygon::range(), Heed::splane::range(), and Heed::straight::vecdistance().