X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fir.h;h=cadc91b8b465f58472e083d9a55f55e8c08a7fe8;hp=7da06ffd46b6da0e7060aedef91bd4d37cd98f0f;hb=2451ac413848718a1dd666ce6f6464e974680f47;hpb=3873067b47875ba7b09d511cb5cdfeaa35e8c27d diff --git a/include/ir.h b/include/ir.h index 7da06ff..cadc91b 100644 --- a/include/ir.h +++ b/include/ir.h @@ -1,75 +1,15 @@ /***************************************************************************** -** This is part of the CTSim program -** Copyright (C) 1983-2000 Kevin Rosenberg -** -** $Id: ir.h,v 1.22 2000/06/07 07:43:39 kevin Exp $ -** $Log: ir.h,v $ -** Revision 1.22 2000/06/07 07:43:39 kevin -** *** empty log message *** -** -** Revision 1.21 2000/06/07 03:49:54 kevin -** *** empty log message *** -** -** Revision 1.20 2000/06/07 02:30:27 kevin -** Added C++ image files -** -** Revision 1.19 2000/06/07 00:59:38 kevin -** added imagefiles -** -** Revision 1.18 2000/06/05 01:32:45 kevin -** Added C++ compatibility -** -** Revision 1.17 2000/06/03 06:29:08 kevin -** *** empty log message *** -** -** Revision 1.16 2000/05/24 22:48:17 kevin -** First functional version of SDF library for X-window -** -** Revision 1.15 2000/05/16 04:33:17 kevin -** Updated documentation -** -** Revision 1.14 2000/05/11 14:07:00 kevin -** Added support for Windows NT -** -** Revision 1.13 2000/05/08 20:00:48 kevin -** ANSI C changes -** -** Revision 1.12 2000/05/07 12:46:19 kevin -** made c++ compatible -** -** Revision 1.11 2000/05/05 02:37:31 kevin -** renamed phmelm to pelm -** -** Revision 1.10 2000/05/04 18:16:34 kevin -** renamed filter definitions -** -** Revision 1.9 2000/05/04 04:29:18 kevin -** *** empty log message *** -** -** Revision 1.8 2000/05/04 04:25:55 kevin -** Renamed phantom and phantom-element functions/variables -** -** Revision 1.7 2000/05/03 19:51:41 kevin -** function renaming for phantoms and phantom elements +** FILE IDENTIFICATION ** -** Revision 1.6 2000/05/03 08:49:49 kevin -** Code cleanup +** Name: ir.h +** Purpose: Master Image Reconstruction Header +** Programmer: Kevin Rosenberg +** Date Started: July 1, 1984 ** -** Revision 1.5 2000/05/02 20:00:25 kevin -** *** empty log message *** -** -** Revision 1.4 2000/05/02 15:31:39 kevin -** code cleaning -** -** Revision 1.3 2000/04/29 23:24:29 kevin -** *** empty log message *** -** -** Revision 1.2 2000/04/28 14:14:16 kevin -** *** empty log message *** -** -** Revision 1.1.1.1 2000/04/28 13:02:43 kevin -** Initial CVS import for first public release +** This is part of the CTSim program +** Copyright (C) 1983-2000 Kevin Rosenberg ** +** $Id: ir.h,v 1.24 2000/06/09 01:35:33 kevin Exp $ ** ** ** This program is free software; you can redistribute it and/or modify @@ -86,11 +26,6 @@ ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA ******************************************************************************/ -/* ir.h Header File for Image Reconstruction System - * Programmer: Kevin Rosenberg - * Date Started: 7-1-84 - */ - #ifndef IR_H #define IR_H @@ -348,47 +283,22 @@ const static int RAYSUM_TRACE_ROW_ATTEN=18; * FUNCTION DECLARATIONS ************************************************************************/ -#ifdef __cplusplus -extern "C" { -#endif /* __cplusplus */ - -#ifdef MPI_CT -#define MPI_MAX_PROCESS 128 -struct mpi_ct_st -{ - int my_rank; - int nproc; - int base_local_work_units; - int remainder_work_units; - int local_work_units[MPI_MAX_PROCESS]; - int start_work_unit[MPI_MAX_PROCESS]; - MPI_Comm comm; -}; - -extern struct mpi_ct_st mpi_ct; -void mpi_ct_calc_work_units(const unsigned int global_work_units); -#endif - - -/* bspline.c */ -int bspline(int samples, int zoom_factor, int spline_order, double input[], double output[]); - /* convolve.c */ -double convolve(const double f1[], const double f2[], const double dx, const int n, const int np, const FunctionSymmetry func_type); -double convolve_both(const double f1[], const double f2[], const double dx, const int n, const int np); +double convolve (const double f1[], const double f2[], const double dx, const int n, const int np, const FunctionSymmetry func_type); +double convolve_both (const double f1[], const double f2[], const double dx, const int n, const int np); /* dialogs.c */ -int phm_add_pelm_kb(PHANTOM *phm); -PHANTOM *phm_select(void); -int interpolation_select(void); -int filter_select(double *filter_param); +int phm_add_pelm_kb (PHANTOM *phm); +PHANTOM *phm_select (void); +int interpolation_select (void); +int filter_select (double *filter_param); /* filter.c */ -double *filter_generate(const FilterType filt_type, double bw, double xmin, double xmax, int n, double param, const DomainType domain, int numint); -double filter_spatial_response_calc(int filt_type, double x, double bw, double param, int n); -double filter_spatial_response_analytic(int filt_type, double x, double bw, double param); -double filter_frequency_response(int filt_type, double u, double bw, double param); -double sinc(double x, double mult); +double *filter_generate (const FilterType filt_type, double bw, double xmin, double xmax, int n, double param, const DomainType domain, int numint); +double filter_spatial_response_calc (int filt_type, double x, double bw, double param, int n); +double filter_spatial_response_analytic (int filt_type, double x, double bw, double param); +double filter_frequency_response (int filt_type, double u, double bw, double param); +double sinc (double x, double mult); double integral_abscos(double u, double w); /* options.c */ @@ -435,11 +345,11 @@ int raysum_collect(RAYSUM *rs, const DETECTOR *det, const PHANTOM *phm, const in void rayview(const PHANTOM *phm, DETARRAY *darray, const DETECTOR *det, const double xd1, const double yd1, const double xd2, const double yd2, const double xs1, const double ys1, const double xs2, const double ys2, const int unit_pulse); double phm_ray_attenuation (const PHANTOM *phm, const double x1, const double y1, const double x2, const double y2); double pelm_ray_attenuation (PELM *pelm, const double x1, const double y1, const double x2, const double y2); -int pelm_clip_line (const PELM *pelm, double *x1, double *y1, double *x2, double *y2); +int pelm_clip_line (const PELM *pelm, double& x1, double& y1, double& x2, double& y2); void raysum_trace_show_param (const char *label, const char *fmt, int row, int color, ...); /* scanner.c */ -DETECTOR *detector_create(const PHANTOM *phm, int geometry, int ndet, int nview, int nsample, const double rot_anglen); +DETECTOR *detector_create(const PHANTOM *phm, const ScannerGeometry geometry, int ndet, int nview, int nsample, const double rot_anglen); void detector_free(DETECTOR *det); /* rayio.c */ @@ -460,9 +370,6 @@ int detarray_read(RAYSUM *rs, DETARRAY *darray, const int view_num); int detarray_write(RAYSUM *rs, const DETARRAY *darray, const int view_num); int raysum_print(const RAYSUM *rs); -#ifdef __cplusplus -} - /* From phm2image.cpp */ void phm_to_imagefile (const PHANTOM *phm, ImageFile& im, const int col_start, const int col_count, const int nsample, const int trace); int pelm_is_point_inside(PELM *obj, const double x, const double y, const CoordType coord_type); @@ -505,8 +412,4 @@ int backproj_calc_id2 (const RAYSUM *rs, ImageFile& im, const double *t, const double view_angle, const int interp_type); void backproj_term_id2 (const RAYSUM *rs, ImageFile& im); - -#endif /* __cplusplus */ - - #endif