upgrade AST

1 files changed, 109 insertions, 0 deletions

diff --git a/ast/finterval.c b/ast/finterval.c
new file mode 100644
index 0000000..7023868
--- /dev/null
+++ b/ast/finterval.c
@@ -0,0 +1,109 @@
+/*
+*+
+*  Name:
+*     finterval.c
+
+*  Purpose:
+*     Define a FORTRAN 77 interface to the AST Interval 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 Interval class.
+
+*  Routines Defined:
+*     AST_ISAINTERVAL
+*     AST_INTERVAL
+
+*  Copyright:
+*     Copyright (C) 1997-2006 Council for the Central Laboratory of the
+*     Research Councils
+
+*  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:
+*     31-AUG-2004 (DSB):
+*        Original version.
+*/
+
+/* 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 "interval.h"              /* C interface to the Interval class */
+
+
+F77_LOGICAL_FUNCTION(ast_isainterval)( INTEGER(THIS), INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   F77_LOGICAL_TYPE(RESULT);
+
+   astAt( "AST_ISAINTERVAL", NULL, 0 );
+   astWatchSTATUS(
+      RESULT = astIsAInterval( astI2P( *THIS ) ) ? F77_TRUE : F77_FALSE;
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_interval)( INTEGER(FRAME),
+                                    DOUBLE_ARRAY(LBND),
+                                    DOUBLE_ARRAY(UBND),
+                                    INTEGER(UNC),
+                                    CHARACTER(OPTIONS),
+                                    INTEGER(STATUS)
+                                    TRAIL(OPTIONS) ) {
+   GENPTR_INTEGER(FRAME)
+   GENPTR_DOUBLE_ARRAY(LBND)
+   GENPTR_DOUBLE_ARRAY(UBND)
+   GENPTR_INTEGER(UNC)
+   GENPTR_CHARACTER(OPTIONS)
+   F77_INTEGER_TYPE(RESULT);
+   char *options;
+   int i;
+
+   astAt( "AST_INTERVAL", 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( astInterval( astI2P( *FRAME ), LBND, UBND,
+                                    astI2P( *UNC ), "%s", options ) );
+      astFree( options );
+   )
+   return RESULT;
+}