Changeset 921 for trunk/source/graphics_reps/src/HepPolyhedron.cc
- Timestamp:
- Feb 16, 2009, 10:14:30 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/graphics_reps/src/HepPolyhedron.cc
r850 r921 25 25 // 26 26 // 27 // $Id: HepPolyhedron.cc,v 1.3 1 2008/04/28 16:06:06 allisonExp $28 // GEANT4 tag $Name: HEAD$27 // $Id: HepPolyhedron.cc,v 1.32 2008/11/13 09:05:27 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-02-cand-01 $ 29 29 // 30 30 // … … 1630 1630 if (sqrtan1<0.||sqrtan2<0.) k += 4; 1631 1631 1632 if (k != 0) { 1632 if (k != 0) 1633 { 1633 1634 std::cerr << "HepPolyhedronHype: error in input parameters"; 1634 1635 if ((k & 1) != 0) std::cerr << " (radiuses)"; … … 1650 1651 double k2 = r2*r2; 1651 1652 1652 double *zz = new double[n + n], *rr = new double[n + n];1653 double *zz = new double[n+n+1], *rr = new double[n+n+1]; 1653 1654 1654 1655 zz[0] = halfZ; 1655 1656 rr[0] = std::sqrt(sqrtan2*halfZ*halfZ+k2); 1656 1657 1657 for(int i = 1; i < n -1; i++)1658 for(int i = 1; i < n-1; i++) 1658 1659 { 1659 1660 zz[i] = zz[i-1] - dz; 1660 rr[i] =std::sqrt(sqrtan2*zz[i]*zz[i]+k2) ; 1661 1661 rr[i] =std::sqrt(sqrtan2*zz[i]*zz[i]+k2); 1662 1662 } 1663 1663 … … 1665 1665 rr[n-1] = rr[0]; 1666 1666 1667 1668 1667 zz[n] = halfZ; 1669 1668 rr[n] = std::sqrt(sqrtan1*halfZ*halfZ+k1); 1670 for(int i = n+1; i < n +n; i++) 1669 1670 for(int i = n+1; i < n+n; i++) 1671 1671 { 1672 1672 zz[i] = zz[i-1] - dz; 1673 rr[i] =std::sqrt(sqrtan1*zz[i]*zz[i]+k1) ; 1674 1673 rr[i] =std::sqrt(sqrtan1*zz[i]*zz[i]+k1); 1675 1674 } 1676 1675 zz[n+n] = -halfZ;
Note: See TracChangeset
for help on using the changeset viewer.