pool/python-sherpa
SHA256
11
0
Fork 1
Code Issues Pull Requests Activity
Files
devel
python-sherpa/support-gcc-15.patch
Steve Kowalik 755cbffa2f - Update to 4.17.1: 
* enhancements:
    + improved handling for OGIP products that do not follow the standards
    + Models can now be combined using arbitrary ufuncs
    + various updates to XSpec interface and support for XSpec 12.15.0
      including the additional of 17 new models
    + added interface to support ArviZ for analysis of Bayesian models
  * infrastructure and testing:
    + updated tests to work with python 3.13
  * bug fixes:
    + fixed a bug in the caching code that limited the cache size to 1
    + update to use fake_pha for missions which don't use 1 as the first
      channel
- Dropped patches:
  * sherpa-pr2188-np2docstrings.patch
  * sherpa-pr2207-mpl.patch
  * support-pytest-8.3.4.patch
- Refresed sherpa-suse-libdir.patch
- Add patch support-gcc-15.patch:
  * Specify function declarations correctly.

OBS-URL: https://build.opensuse.org/package/show/devel:languages:python:numeric/python-sherpa?expand=0&rev=52
2025-07-30 01:15:55 +00:00

536 lines
18 KiB
Diff
Raw Permalink Blame History

Index: sherpa-4.17.1/extern/region-4.14/src/region_eval.c
===================================================================
--- sherpa-4.17.1.orig/extern/region-4.14/src/region_eval.c
+++ sherpa-4.17.1/extern/region-4.14/src/region_eval.c
@@ -80,30 +80,28 @@ int regCompareRegion( regRegion* Region1
{
regShape* Shape1;
regShape* Shape2;
- int true = 1;
- int false = 0;
Shape1 = Region1->shape;
Shape2 = Region2->shape;
while (Shape1 != NULL )
{
if ( !Shape2 )
- return false;
+ return 0;
if ( Shape1->component != Shape2->component )
- return false;
+ return 0;
if( !Shape1->isEqual(Shape1, Shape2) )
- return false;
+ return 0;
Shape1 = Shape1->next;
Shape2 = Shape2->next;
}
if ( Shape2 )
- return false;
+ return 0;
- return true;
+ return 1;
}
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/imhfile.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/imhfile.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/imhfile.c
@@ -148,27 +148,27 @@
#define LEN_FITSHDR 11520
int check_immagic();
-int irafgeti4();
-float irafgetr4();
-char *irafgetc2();
-char *irafgetc();
-char *iraf2str();
-static char *same_path();
-static void irafputr4();
-static void irafputi4();
-static void irafputc2();
-static void irafputc();
-static void str2iraf();
+int irafgeti4(char *irafheader, int offset);
+float irafgetr4(char *irafheader, int offset);
+char *irafgetc2(char *irafheader, int offset, int nc);
+char *irafgetc(char *irafheader, int offset, int nc);
+char *iraf2str(char *irafstring, int nchar);
+static char *same_path(char *pixname, char *hdrname);
+static void irafputr4(char *irafheader, int offset, float rnum);
+static void irafputi4(char *irafheader, int offset, int inum);
+static void irafputc2(char *string, char *irafheader, int offset, int nc);
+static void irafputc(char *string, char *irafheader, int offset, int nc);
+static void str2iraf(char *string, char *irafstring, int nchar);
static int headswap=-1; /* =1 to swap data bytes of foreign IRAF file */
-static void irafswap();
-static void irafswap2();
-static void irafswap4();
-static void irafswap8();
-int head_version ();
-int pix_version ();
-int irafncmp ();
+static void irafswap(int bitpix, char *string, int nbytes);
+static void irafswap2(char *string, int nbytes);
+static void irafswap4(char *string, int nbytes);
+static void irafswap8(char *string, int nbytes);
+int head_version(char *irafheader);
+int pix_version(char *irafheader);
+int irafncmp(char *irafheader, char *teststring, int nc);
static int machswap();
-static int irafsize();
+static int irafsize(FILE *diskfile);
#define SECONDS_1970_TO_1980 315532800L
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/hget.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/hget.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/hget.c
@@ -74,7 +74,7 @@
static int use_saolib=0;
#endif
-char *hgetc ();
+char *hgetc(const char *hstring, const char *keyword0);
static char val[VLENGTH+1];
static int multiline = 0;
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/hput.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/hput.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/hput.c
@@ -58,7 +58,7 @@
static int verbose=0; /* Set to 1 to print error messages and other info */
-static void fixnegzero();
+static void fixnegzero(char *string);
/* HPUTI4 - Set int keyword = ival in FITS header string */
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/iget.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/iget.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/iget.c
@@ -52,7 +52,7 @@
#define MAX_LVAL 2000
-static char *isearch();
+static char *isearch(const char *hstring, const char *keyword);
static char val[30];
/* Extract long value for variable from IRAF multiline keyword value */
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcscon.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/wcscon.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcscon.c
@@ -73,10 +73,30 @@
#include <string.h>
#include "wcs.h"
-void fk524(), fk524e(), fk524m(), fk524pv();
-void fk425(), fk425e(), fk425m(), fk425pv();
-void fk42gal(), fk52gal(), gal2fk4(), gal2fk5();
-void fk42ecl(), fk52ecl(), ecl2fk4(), ecl2fk5();
+void fk524(double *ra, double *dec);
+void fk524e(double *ra, double *dec, double epoch);
+void fk524m(double *ra, double *dec, double *rapm, double *decpm);
+void fk524pv(double *ra, double *dec, double *rapm, double *decpm, double *parallax, double *rv);
+void fk425(double *ra, double *dec);
+void fk425e(double *ra, double *dec, double epoch);
+void fk425m(double *ra, double *dec, double *rapm, double *decpm);
+void fk425pv(double *ra, double *dec, double *rapm, double *decpm, double *parallax, double *rv);
+void fk42gal(double *dtheta, double *dphi);
+void fk52gal(double *dtheta, double *dphi);
+void gal2fk4(double *dtheta, double *dphi);
+void gal2fk5(double *dtheta, double *dphi);
+void fk42ecl(double *dtheta, double *dphi, double epoch);
+void fk52ecl(double *dtheta, double *dphi, double epoch);
+void ecl2fk4(double *dtheta, double *dphi, double epoch);
+void ecl2fk5(double *dtheta, double *dphi, double epoch);
+void v2s3(double pos[3], double *rra, double *rdec, double *r);
+void s2v3(double rra, double rdec, double r, double pos[3]);
+void fk4prec(double ep0, double ep1, double *ra, double *dec);
+void fk5prec(double ep0, double ep1, double *ra, double *dec);
+void rotmat(int axes, double rot1, double rot2, double rot3, double *matrix);
+void mprecfk4(double bep0, double bep1, double rmatp[9]);
+void mprecfk5(double ep0, double ep1, double rmatp[9]);
+char *eqstrn(double dra, double ddec);
/* Convert from coordinate system sys1 to coordinate system sys2, converting
proper motions, too, and adding them if an epoch is specified */
@@ -100,7 +120,6 @@ double *pphi; /* Latitude or declination
Input in sys1, returned in sys2 */
{
- void fk5prec(), fk4prec();
/* Set equinoxes if 0.0 */
if (eq1 == 0.0) {
@@ -306,7 +325,6 @@ double *px; /* Parallax in arcseconds */
double *rv; /* Radial velocity in km/sec */
{
- void fk5prec(), fk4prec();
/* Set equinoxes if 0.0 */
if (eq1 == 0.0) {
@@ -510,7 +528,6 @@ double *dphi; /* Latitude or declination
double epoch; /* Besselian epoch in years */
{
- void fk5prec(), fk4prec();
/* Set equinoxes if 0.0 */
if (eq1 == 0.0) {
@@ -1441,9 +1458,8 @@ double *dphi; /* B1950.0 FK4 declination
Reference: Blaauw et al, MNRAS,121,123 (1960) */
{
double pos[3],pos1[3],r,dl,db,rl,rb,rra,rdec,dra,ddec;
- void v2s3(),s2v3();
int i;
- char *eqcoor, *eqstrn();
+ char *eqcoor;
dra = *dtheta;
ddec = *dphi;
@@ -1499,8 +1515,7 @@ double *dphi; /* Galactic latitude (b2)
{
double pos[3],pos1[3],r,dl,db,rl,rb,rra,rdec,dra,ddec;
- void v2s3(),s2v3();
- char *eqcoor, *eqstrn();
+ char *eqcoor;
int i;
/* spherical to cartesian */
@@ -1574,8 +1589,7 @@ double *dphi; /* J2000 declination in de
Reference: Blaauw et al, MNRAS,121,123 (1960) */
{
double pos[3],pos1[3],r,dl,db,rl,rb,rra,rdec,dra,ddec;
- void v2s3(),s2v3();
- char *eqcoor, *eqstrn();
+ char *eqcoor;
int i;
/* Spherical to cartesian */
@@ -1627,9 +1641,8 @@ double *dphi; /* Galactic latitude (b2)
{
double pos[3],pos1[3],r,dl,db,rl,rb,rra,rdec,dra,ddec;
- void v2s3(),s2v3();
int i;
- char *eqcoor, *eqstrn();
+ char *eqcoor;
/* Spherical to Cartesian */
dl = *dtheta;
@@ -1727,8 +1740,6 @@ double *dphi; /* B1950 declination in de
double epoch; /* Besselian epoch in years */
{
- void fk425e(), fk52ecl();
-
/* Convert from B1950 to J2000 coordinates */
fk425e (dtheta, dphi, epoch);
@@ -1755,8 +1766,6 @@ double epoch; /* Besselian epoch in year
double v1[3], v2[3], r;
double rmat[9], *rmati; /* Rotation matrix */
- void rotmat(), v2s3(), s2v3(), fk5prec();
-
/* Precess coordinates from J2000 to epoch */
if (epoch != 2000.0)
fk5prec (2000.0, epoch, dtheta, dphi);
@@ -1810,8 +1819,6 @@ double *dphi; /* Galactic latitude (b2)
double epoch; /* Besselian epoch in years */
{
- void ecl2fk5(), fk524e();
-
/* Convert from ecliptic to J2000 coordinates */
ecl2fk5 (dtheta, dphi, epoch);
@@ -1839,7 +1846,6 @@ double epoch; /* Besselian epoch in year
double rtheta, rphi, v1[3], v2[3];
double t, eps0, r;
double rmat[9]; /* Rotation matrix */
- void v2s3(),s2v3(), fk5prec(), rotmat();
rtheta = degrad (*dtheta);
rphi = degrad (*dphi);
@@ -1902,7 +1908,6 @@ double *dec; /* Dec in degrees mean equa
{
int i, j;
double pm[9], *pmi, v1[3], v2[3], rra, rdec, r;
- void v2s3(),s2v3(), mprecfk4();
rra = degrad (*ra);
rdec = degrad (*dec);
@@ -1950,7 +1955,6 @@ double *dec; /* Dec in degrees mean equa
{
int i, j;
double pm[9], *pmi, v1[3], v2[3], rra, rdec, r;
- void v2s3(),s2v3(), mprecfk5();
rra = degrad (*ra);
rdec = degrad (*dec);
@@ -2003,7 +2007,6 @@ double rmatp[9]; /* 3x3 Precession matri
*/
{
double bigt, t, tas2r, w, zeta, z, theta;
- void rotmat();
/* Interval between basic epoch B1850.0 and beginning epoch in TC */
bigt = ( bep0 - 1850.0 ) / 100.0;
@@ -2047,7 +2050,6 @@ double rmatp[9]; /* 3x3 Precession matri
*/
{
double t0, t, tas2r, w, zeta, z, theta;
- void rotmat();
/* Interval between basic epoch J2000.0 and beginning epoch (JC) */
t0 = ( ep0 - 2000.0 ) / 100.0;
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/tnxpos.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/tnxpos.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/tnxpos.c
@@ -47,10 +47,17 @@
#define max_niter 500
#define SZ_ATSTRING 2000
-static void wf_gsclose();
-static void wf_gsb1pol();
-static void wf_gsb1leg();
-static void wf_gsb1cheb();
+static void wf_gsclose(struct IRAFsurface *sf);
+static void wf_gsb1pol(double x, int order, double *basis);
+static void wf_gsb1leg(double x, int order, double k1, double k2, double *basis);
+static void wf_gsb1cheb(double x, int order, double k1, double k2, double *basis);
+struct IRAFsurface *wf_gsopen(char *astr);
+struct IRAFsurface *wf_gsrestore(double *fit);
+struct IRAFsurface *wf_gspset(int xorder, int yorder, int xterms, double *coeff);
+double wf_gseval(struct IRAFsurface *sf, double x, double y);
+double wf_gsder(struct IRAFsurface *sf, double x, double y, int nxd, int nyd);
+
+extern void wcsrotset(struct WorldCoor *wcs);
/* tnxinit -- initialize the gnomonic forward or inverse transform.
* initialization for this transformation consists of, determining which
@@ -75,9 +82,7 @@ tnxinit (header, wcs)
const char *header; /* FITS header */
struct WorldCoor *wcs; /* pointer to WCS structure */
{
- struct IRAFsurface *wf_gsopen();
char *str1, *str2, *lngstr, *latstr;
- extern void wcsrotset();
/* allocate space for the attribute strings */
str1 = malloc (SZ_ATSTRING);
@@ -166,7 +171,6 @@ double *xpos, *ypos; /*o world coordinat
double x, y, r, phi, theta, costhe, sinthe, dphi, cosphi, sinphi, dlng, z;
double colatp, coslatp, sinlatp, longp;
double xs, ys, ra, dec, xp, yp;
- double wf_gseval();
/* Convert from pixels to image coordinates */
xpix = xpix - wcs->crpix[0];
@@ -317,7 +321,6 @@ double *xpix, *ypix; /*o physical coordi
double s, r, dphi, z, dpi, dhalfpi, twopi, tx;
double xm, ym, f, fx, fy, g, gx, gy, denom, dx, dy;
double colatp, coslatp, sinlatp, longp, sphtol;
- double wf_gseval(), wf_gsder();
/* get the axis numbers */
if (wcs->coorflip) {
@@ -562,7 +565,6 @@ char *astr; /* the input mwcs attrib
int npar, szcoeff;
double *coeff;
struct IRAFsurface *gs;
- struct IRAFsurface *wf_gsrestore();
if (astr[1] == 0)
return (NULL);
@@ -737,7 +739,6 @@ int nxd, nyd; /* order of the derivative
struct IRAFsurface *sf2 = 0;
double *ptr1, *ptr2;
double zfit, norm;
- double wf_gseval();
if (sf1 == NULL)
return (0.0);
@@ -1121,7 +1122,6 @@ double *coeff; /* Plate fit coefficient
{
double *ycoeff;
- struct IRAFsurface *wf_gspset ();
wcs->prjcode = WCS_TNX;
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/zpxpos.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/zpxpos.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/zpxpos.c
@@ -50,7 +50,13 @@
#define max_niter 500
#define SZ_ATSTRING 2000
-static void wf_gsclose();
+
+static void wf_gsclose(struct IRAFsurface *sf);
+
+extern struct IRAFsurface *wf_gsopen(char *astr);
+extern double wf_gseval(struct IRAFsurface *sf, double x, double y);
+extern double wf_gsder(struct IRAFsurface *sf, double x, double y, int nxd, int nyd);
+extern void wcsrotset(struct WorldCoor *wcs);
/* zpxinit -- initialize the zenithal/azimuthal polynomial forward or
* inverse transform. initialization for this transformation consists of,
@@ -81,10 +87,8 @@ const char *header; /* FITS header */
struct WorldCoor *wcs; /* pointer to WCS structure */
{
int i, j;
- struct IRAFsurface *wf_gsopen();
char key[8], *str1, *str2, *lngstr, *latstr, *header1;
double zd1, d1, zd2,d2, zd, d, r;
- extern void wcsrotset();
/* allocate space for the attribute strings */
str1 = malloc (SZ_ATSTRING);
@@ -252,7 +256,6 @@ double *xpos, *ypos; /*o world coordinat
double colatp, coslatp, sinlatp, longp;
double xs, ys, ra, dec, xp, yp;
double a, b, c, d, zd, zd1, zd2, r1, r2, rt, lambda;
- double wf_gseval();
/* Convert from pixels to image coordinates */
xpix = xpix - wcs->crpix[0];
@@ -520,7 +523,6 @@ double *xpix, *ypix; /*o physical coordi
double s, r, dphi, z, dpi, dhalfpi, twopi, tx;
double xm, ym, f, fx, fy, g, gx, gy, denom, dx, dy;
double colatp, coslatp, sinlatp, longp, sphtol;
- double wf_gseval(), wf_gsder();
/* get the axis numbers */
if (wcs->coorflip) {
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcsinit.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/wcsinit.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcsinit.c
@@ -47,11 +47,13 @@
#include <stdlib.h>
#endif
-static void wcseq();
-static void wcseqm();
-static void wcsioset();
-void wcsrotset();
-char wcschar();
+void invert_wcs(struct WorldCoor *wcs);
+static void wcseq(char *hstring, struct WorldCoor *wcs);
+static void wcseqm(char *hstring, struct WorldCoor *wcs, char *mchar);
+static void wcsioset(struct WorldCoor *wcs);
+char wcschar(const char *hstring, const char *name);
+
+extern void wcsrotset(struct WorldCoor *wcs);
/* set up a WCS structure from a FITS image header lhstring bytes long
* for a specified WCS name */
@@ -242,11 +244,6 @@ char *wchar; /* Suffix character for on
double ut;
int nax;
int twod;
- extern int tnxinit();
- extern int zpxinit();
- extern int platepos();
- extern int dsspos();
- void invert_wcs();
wcs = (struct WorldCoor *) calloc (1, sizeof(struct WorldCoor));
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcs.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/wcs.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/wcs.c
@@ -89,12 +89,16 @@
static char wcserrmsg[80];
static char wcsfile[256]={""};
-static void wcslibrot();
-void wcsrotset();
static int wcsproj0 = 0;
static int izpix = 0;
static double zpix = 0.0;
+static void wcslibrot(struct WorldCoor *wcs);
+int wcspos(double xpix, double ypix, struct WorldCoor *wcs, double *xpos, double *ypos);
+int wcspix(double xpos, double ypos, struct WorldCoor *wcs, double *xpix, double *ypix);
+
+extern void wcsrotset(struct WorldCoor *wcs);
+
void
wcsfree (wcs)
struct WorldCoor *wcs; /* WCS structure */
@@ -2123,7 +2127,6 @@ double *xpos,*ypos; /* RA and Dec in deg
{
double xpi, ypi, xp, yp;
double eqin, eqout;
- int wcspos();
if (nowcs (wcs))
return;
@@ -2246,7 +2249,6 @@ int *offscl; /* 0 if within bounds, else
double xp, yp, xpi, ypi;
double eqin, eqout;
int sysin;
- int wcspix();
if (nowcs (wcs))
return;
@@ -2351,7 +2353,6 @@ double *ypos; /* y (dec) coor
{
int offscl;
int i;
- int wcsrev();
double wcscrd[4], imgcrd[4], pixcrd[4];
double phi, theta;
@@ -2392,7 +2393,6 @@ double *ypix; /* y pixel numbe
{
int offscl;
- int wcsfwd();
double wcscrd[4], imgcrd[4], pixcrd[4];
double phi, theta;
Index: sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/dateutil.c
===================================================================
--- sherpa-4.17.1.orig/extern/wcssubs/wcssubs-3.9.5/dateutil.c
+++ sherpa-4.17.1/extern/wcssubs/wcssubs-3.9.5/dateutil.c
@@ -313,11 +313,11 @@
#include "wcs.h"
#include "fitsfile.h"
-static double suntl();
-static void fixdate();
-static int caldays();
-static double dint();
-static double dmod();
+static double suntl(double dj, double ra, double dec, int sys);
+static void fixdate(int *iyr, int *imon, int *iday, int *ihr, int *imn, double *sec, int ndsec);
+static int caldays(int year, int month);
+static double dint(double dnum);
+static double dmod(double dnum, double dm);
static double longitude = 0.0; /* longitude of observatory in degrees (+=west) */
void
Reference in New Issue View Git Blame Copy Permalink