upgrade ast 8.7.1

1 files changed, 226 insertions, 0 deletions

diff --git a/ast/fpolygon.c b/ast/fpolygon.c
new file mode 100644
index 0000000..12247d7
--- /dev/null
+++ b/ast/fpolygon.c
@@ -0,0 +1,226 @@
+/*
+*+
+*  Name:
+*     fpolygon.c
+
+*  Purpose:
+*     Define a FORTRAN 77 interface to the AST Polygon class.
+
+*  Type of Module:
+*     C source file.
+
+*  Description:
+*     This file defines FORTRAN 77-callable C functions which provide
+*     a public FORTRAN 77 interface to the Polygon class.
+
+*  Routines Defined:
+*     AST_ISAPOLYGON
+*     AST_POLYGON
+*     AST_DOWNSIZE
+
+*  Copyright:
+*     Copyright (C) 1997-2006 Council for the Central Laboratory of the
+*     Research Councils
+*     Copyright (C) 2009 Science & Technology Facilities Council.
+*     All Rights Reserved.
+
+*  Licence:
+*     This program is free software: you can redistribute it and/or
+*     modify it under the terms of the GNU Lesser General Public
+*     License as published by the Free Software Foundation, either
+*     version 3 of the License, or (at your option) any later
+*     version.
+*
+*     This program is distributed in the hope that it will be useful,
+*     but WITHOUT ANY WARRANTY; without even the implied warranty of
+*     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+*     GNU Lesser General Public License for more details.
+*
+*     You should have received a copy of the GNU Lesser General
+*     License along with this program.  If not, see
+*     <http://www.gnu.org/licenses/>.
+
+*  Authors:
+*     DSB: David S. Berry (Starlink)
+
+*  History:
+*     27-OCT-2004 (DSB):
+*        Original version.
+*     28-MAY-2009 (DSB):
+*        Added AST_DOWNSIZE.
+*/
+
+/* Define the astFORTRAN77 macro which prevents error messages from
+   AST C functions from reporting the file and line number where the
+   error occurred (since these would refer to this file, they would
+   not be useful). */
+#define astFORTRAN77
+
+/* Header files. */
+/* ============= */
+#include "f77.h"                 /* FORTRAN <-> C interface macros (SUN/209) */
+#include "c2f77.h"               /* F77 <-> C support functions/macros */
+#include "error.h"               /* Error reporting facilities */
+#include "memory.h"              /* Memory handling facilities */
+#include "polygon.h"             /* C interface to the Polygon class */
+
+F77_LOGICAL_FUNCTION(ast_isapolygon)( INTEGER(THIS), INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   F77_LOGICAL_TYPE(RESULT);
+
+   astAt( "AST_ISAPOLYGON", NULL, 0 );
+   astWatchSTATUS(
+      RESULT = astIsAPolygon( astI2P( *THIS ) ) ? F77_TRUE : F77_FALSE;
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_polygon)( INTEGER(FRAME),
+                                   INTEGER(NPNT),
+                                   INTEGER(INDIM),
+                                   DOUBLE_ARRAY(POINTS),
+                                   INTEGER(UNC),
+                                   CHARACTER(OPTIONS),
+                                   INTEGER(STATUS)
+                                   TRAIL(OPTIONS) ) {
+   GENPTR_INTEGER(FRAME)
+   GENPTR_INTEGER(NPNT)
+   GENPTR_INTEGER(INDIM)
+   GENPTR_DOUBLE_ARRAY(POINTS)
+   GENPTR_INTEGER(UNC)
+   GENPTR_CHARACTER(OPTIONS)
+   F77_INTEGER_TYPE(RESULT);
+   char *options;
+   int i;
+
+   astAt( "AST_POLYGON", NULL, 0 );
+   astWatchSTATUS(
+      options = astString( OPTIONS, OPTIONS_length );
+
+/* Truncate the options string to exlucde any trailing spaces. */
+      astChrTrunc( options );
+
+/* Change ',' to '\n' (see AST_SET in fobject.c for why). */
+      if ( astOK ) {
+         for ( i = 0; options[ i ]; i++ ) {
+            if ( options[ i ] == ',' ) options[ i ] = '\n';
+         }
+      }
+
+      RESULT = astP2I( astPolygon( astI2P( *FRAME ), *NPNT, *INDIM, POINTS,
+                                   astI2P( *UNC ), "%s", options ) );
+      astFree( options );
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_downsize)( INTEGER(THIS),
+                                    DOUBLE(MAXERR),
+                                    INTEGER(MAXVERT),
+                                    INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_DOUBLE(MAXERR)
+   GENPTR_INTEGER(MAXVERT)
+   F77_INTEGER_TYPE(RESULT);
+
+   astAt( "AST_DOWNSIZE", NULL, 0 );
+   astWatchSTATUS(
+      RESULT = astP2I( astDownsize( astI2P( *THIS ), *MAXERR, *MAXVERT ) );
+   )
+   return RESULT;
+}
+
+
+/* AST_OUTLINE<X> requires a function for each possible data type, so
+   define it via a macro. */
+#define MAKE_AST_OUTLINE(f,F,Ftype,X,Xtype) \
+F77_INTEGER_FUNCTION(ast_outline##f)( Ftype(VALUE), \
+                                      INTEGER(OPER), \
+                                      Ftype##_ARRAY(ARRAY), \
+                                      INTEGER_ARRAY(LBND), \
+                                      INTEGER_ARRAY(UBND), \
+                                      DOUBLE(MAXERR), \
+                                      INTEGER(MAXVERT), \
+                                      INTEGER_ARRAY(INSIDE), \
+                                      LOGICAL(STARPIX), \
+                                      INTEGER(STATUS) ) { \
+   GENPTR_##Ftype(VALUE) \
+   GENPTR_INTEGER(OPER) \
+   GENPTR_##Ftype##_ARRAY(ARRAY) \
+   GENPTR_INTEGER_ARRAY(LBND) \
+   GENPTR_INTEGER_ARRAY(UBND) \
+   GENPTR_DOUBLE(MAXERR) \
+   GENPTR_INTEGER(MAXVERT) \
+   GENPTR_INTEGER_ARRAY(INSIDE) \
+   GENPTR_LOGICAL(STARPIX) \
+   GENPTR_INTEGER(STATUS) \
+\
+   F77_INTEGER_TYPE RESULT; \
+\
+   astAt( "AST_OUTLINE"#F, NULL, 0 ); \
+   astWatchSTATUS( \
+      RESULT = astP2I( astOutline##X( *VALUE, *OPER, (Xtype *) ARRAY, LBND, \
+                                       UBND, *MAXERR, *MAXVERT, INSIDE, \
+                                       F77_ISTRUE( *STARPIX ) ? 1 : 0 ) ); \
+   ) \
+   return RESULT; \
+}
+
+/* Invoke the above macro to define a function for each data
+   type. Include synonyms for some functions. */
+MAKE_AST_OUTLINE(d,D,DOUBLE,D,double)
+MAKE_AST_OUTLINE(r,R,REAL,F,float)
+MAKE_AST_OUTLINE(i,I,INTEGER,I,int)
+MAKE_AST_OUTLINE(ui,UI,INTEGER,UI,unsigned int)
+MAKE_AST_OUTLINE(s,S,WORD,S,short int)
+MAKE_AST_OUTLINE(us,US,UWORD,US,unsigned short int)
+MAKE_AST_OUTLINE(w,W,WORD,S,short int)
+MAKE_AST_OUTLINE(uw,UW,UWORD,US,unsigned short int)
+MAKE_AST_OUTLINE(b,B,BYTE,B,signed char)
+MAKE_AST_OUTLINE(ub,UB,UBYTE,UB,unsigned char)
+#undef MAKE_AST_OUTLINE
+
+
+/* AST_CONVEX<X> requires a function for each possible data type, so
+   define it via a macro. */
+#define MAKE_AST_CONVEX(f,F,Ftype,X,Xtype) \
+F77_INTEGER_FUNCTION(ast_convex##f)( Ftype(VALUE), \
+                                      INTEGER(OPER), \
+                                      Ftype##_ARRAY(ARRAY), \
+                                      INTEGER_ARRAY(LBND), \
+                                      INTEGER_ARRAY(UBND), \
+                                      LOGICAL(STARPIX), \
+                                      INTEGER(STATUS) ) { \
+   GENPTR_##Ftype(VALUE) \
+   GENPTR_INTEGER(OPER) \
+   GENPTR_##Ftype##_ARRAY(ARRAY) \
+   GENPTR_INTEGER_ARRAY(LBND) \
+   GENPTR_INTEGER_ARRAY(UBND) \
+   GENPTR_LOGICAL(STARPIX) \
+   GENPTR_INTEGER(STATUS) \
+\
+   F77_INTEGER_TYPE RESULT; \
+\
+   astAt( "AST_CONVEX"#F, NULL, 0 ); \
+   astWatchSTATUS( \
+      RESULT = astP2I( astConvex##X( *VALUE, *OPER, (Xtype *) ARRAY, LBND, \
+                                     UBND, F77_ISTRUE( *STARPIX ) ? 1 : 0 ) ); \
+   ) \
+   return RESULT; \
+}
+
+/* Invoke the above macro to define a function for each data
+   type. Include synonyms for some functions. */
+MAKE_AST_CONVEX(d,D,DOUBLE,D,double)
+MAKE_AST_CONVEX(r,R,REAL,F,float)
+MAKE_AST_CONVEX(i,I,INTEGER,I,int)
+MAKE_AST_CONVEX(ui,UI,INTEGER,UI,unsigned int)
+MAKE_AST_CONVEX(s,S,WORD,S,short int)
+MAKE_AST_CONVEX(us,US,UWORD,US,unsigned short int)
+MAKE_AST_CONVEX(w,W,WORD,S,short int)
+MAKE_AST_CONVEX(uw,UW,UWORD,US,unsigned short int)
+MAKE_AST_CONVEX(b,B,BYTE,B,signed char)
+MAKE_AST_CONVEX(ub,UB,UBYTE,UB,unsigned char)
+#undef MAKE_AST_CONVEX
+
+