SphericMirrorCalculator.cpp   SphericMirrorCalculator.cpp 
skipping to change at line 117 skipping to change at line 117
prod0 += alphaDeltaPhi[3*i+k]*alphaDeltaPhi[3*j+k]; prod0 += alphaDeltaPhi[3*i+k]*alphaDeltaPhi[3*j+k];
prod1 += alphaDeltaPhi[i+3*k]*alphaDeltaPhi[j+3*k]; prod1 += alphaDeltaPhi[i+3*k]*alphaDeltaPhi[j+3*k];
} }
if (i==j) { if (i==j) {
prod0 -= 1.0; prod0 -= 1.0;
prod1 -= 1.0; prod1 -= 1.0;
} }
if (fabs(prod0)>1e-10) { if (fabs(prod0)>1e-10) {
qDebug() << "i: " << i << ", j: " << j qDebug() << "i: " << i << ", j: " << j
<< ", prod0: " << prod0 << ", prod1: " << prod1; << ", prod0: " << prod0 << ", prod1: " << prod1;
assert(0); Q_ASSERT(0);
} }
if (fabs(prod1)>1e-10) { if (fabs(prod1)>1e-10) {
qDebug << "i: " << i << ", j: " << j qDebug << "i: " << i << ", j: " << j
<< ", prod0: " << prod0 << ", prod1: " << prod1; << ", prod0: " << prod0 << ", prod1: " << prod1;
assert(0); Q_ASSERT(0);
} }
} }
} }
*/ */
} }
bool SphericMirrorCalculator::transform(const Vec3d &v, bool SphericMirrorCalculator::transform(const Vec3d &v,
double &xb,double &yb) const { double &xb,double &yb) const {
const Vec3d S = DomeCenter + (v * (DomeRadius/v.length())); const Vec3d S = DomeCenter + (v * (DomeRadius/v.length()));
const Vec3d SmP = S - P; const Vec3d SmP = S - P;
 End of changes. 2 change blocks. 
2 lines changed or deleted 2 lines changed or added

This html diff was produced by rfcdiff 1.41. The latest version is available from http://tools.ietf.org/tools/rfcdiff/