** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: ctsupport.h,v 1.2 2000/06/19 19:04:05 kevin Exp $
+** $Id: ctsupport.h,v 1.3 2000/06/19 19:10:08 kevin Exp $
**
**
** This program is free software; you can redistribute it and/or modify
//////////////////////////////////////////////////////////////
// clip.cpp
-int clip_rect(double& x1, double& y1, double& x2, double& y2, const double rect[4]);
-int clip_segment(double& x1, double& y1, double& x2, double& y2, const double u, const double v);
-int clip_sector(double& x1, double& y1, double& x2, double& y2, const double u, const double v);
-int clip_circle(double& x1, double& y1, double& x2, double& y2, const double cx, const double cy, const double radius, double t1, double t2);
-int clip_triangle(double& x1, double& y1, double& x2, double& y2, const double u, const double v, const int clip_xaxis);
+int clip_rect (double& x1, double& y1, double& x2, double& y2, const double rect[4]);
+int clip_segment (double& x1, double& y1, double& x2, double& y2, const double u, const double v);
+int clip_sector (double& x1, double& y1, double& x2, double& y2, const double u, const double v);
+int clip_circle (double& x1, double& y1, double& x2, double& y2, const double cx, const double cy, const double radius, double t1, double t2);
+int clip_triangle (double& x1, double& y1, double& x2, double& y2, const double u, const double v, const int clip_xaxis);
// norm_ang.cpp
double norm_ang (double theta);
// xform.cpp
-void indent_mtx2(GRFMTX_2D m);
-void xlat_mtx2(GRFMTX_2D m, const double x, const double y);
-void scale_mtx2(GRFMTX_2D m, const double sx, const double sy);
-void rot_mtx2(GRFMTX_2D m, const double theta);
-void mult_mtx2(GRFMTX_2D m1, GRFMTX_2D m2, GRFMTX_2D result);
-void xform_mtx2(GRFMTX_2D m, double& x, double& y);
-void rotate2d(double x[], double y[], int pts, double angle);
-void xlat2d(double x[], double y[], int pts, double xoffset, double yoffset);
-void scale2d(double x[], double y[], int pts, double xfact, double yfact);
+void indent_mtx2 (GRFMTX_2D m);
+void xlat_mtx2 (GRFMTX_2D m, const double x, const double y);
+void scale_mtx2 (GRFMTX_2D m, const double sx, const double sy);
+void rot_mtx2 (GRFMTX_2D m, const double theta);
+void mult_mtx2 (const GRFMTX_2D m1, const GRFMTX_2D m2, GRFMTX_2D result);
+void xform_mtx2 (const GRFMTX_2D m, double& x, double& y);
+void rotate2d (double x[], double y[], int pts, double angle);
+void xlat2d (double x[], double y[], int pts, double xoffset, double yoffset);
+void scale2d (double x[], double y[], int pts, double xfact, double yfact);
// simpson.cpp
double integrateSimpson (const double xmin, const double xmax, const double *y, const int np);
** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: backprojectors.cpp,v 1.1 2000/06/19 02:59:34 kevin Exp $
+** $Id: backprojectors.cpp,v 1.2 2000/06/19 19:07:33 kevin Exp $
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
{
double theta = - view_angle; // add half PI to view angle to get perpendicular theta angle
-#if SIZEOF_LONG == 8
- long int scale = 1 << 32;
-#else
- long int scale = 1 << 16;
-#endif
+ kint32 scale = 1 << 16;
double dScale = scale;
- long int halfScale = scale / 2;
+ kint32 halfScale = scale / 2;
- long int det_dx = nearest<long int> (xInc * sin (theta) / detInc * scale);
- long int det_dy = nearest<long int> (yInc * cos (theta) / detInc * scale);
+ kint32 det_dx = nearest<kint32> (xInc * sin (theta) / detInc * scale);
+ kint32 det_dy = nearest<kint32> (yInc * cos (theta) / detInc * scale);
// calculate L for first point in image (0, 0)
- long int detPosColStart = nearest<long int> (start_r * cos (theta - start_phi) / detInc * scale);
+ kint32 detPosColStart = nearest<kint32> (start_r * cos (theta - start_phi) / detInc * scale);
for (int ix = 0; ix < nx; ix++, detPosColStart += det_dx) {
- long int curDetPos = detPosColStart;
+ kint32 curDetPos = detPosColStart;
ImageFileColumn pImCol = v[ix];
for (int iy = 0; iy < ny; iy++, curDetPos += det_dy) {
else
*pImCol++ += filteredProj[iDetPos];
} else if (interpType == I_LINEAR) {
- long int detPosFloor = curDetPos / scale;
- long int detPosRemainder = curDetPos % scale;
+ kint32 detPosFloor = curDetPos / scale;
+ kint32 detPosRemainder = curDetPos % scale;
if (detPosRemainder < 0) {
detPosFloor--;
detPosRemainder += scale;
** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: xform.cpp,v 1.2 2000/06/19 19:04:05 kevin Exp $
+** $Id: xform.cpp,v 1.3 2000/06/19 19:10:08 kevin Exp $
**
** This program is free software; you can redistribute it and/or modify
** it under the terms of the GNU General Public License (version 2) as
}
void
-mult_mtx2 (GRFMTX_2D m1, GRFMTX_2D m2, GRFMTX_2D result)
+mult_mtx2 (const GRFMTX_2D m1, const GRFMTX_2D m2, GRFMTX_2D result)
{
GRFMTX_2D temp;
}
void
-xform_mtx2 (GRFMTX_2D m, double& x, double& y)
+xform_mtx2 (const GRFMTX_2D m, double& x, double& y)
{
double xt = x * m[0][0] + y * m[1][0] + m[2][0];
double yt = x * m[0][1] + y * m[1][1] + m[2][1];