Message: a Bug in G4ToroidalSurface::Intersect ?? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question a Bug in G4ToroidalSurface::Intersect ?? 

Forum: Geometry
Date: 20 May, 2014
From: Yuefeng Qiu <Yuefeng Qiu>

I think there is a bug for G4BREPSolid/G4ToroidalSurface::Intersect. The code original code are :

  G4double rnorm = MaxRadius - MinRadius; // ei tietoa onko oikein...
  rho = MinRadius*MinRadius / (rnorm*rnorm);

But according to the book Graphic GEMS II (by James Arvo) , considering a and b is the two small radius of the Elliptical torus, we have

rho = (a*a)/(b*b)

For circular torus, a = b, therefore rho should be 1.0 in the code.

Am I correct?

Best regards,

Yuefeng

 Add Message Add Message
to: "a Bug in G4ToroidalSurface::Intersect ??"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]