update ast 8.6.2

1 files changed, 234 insertions, 0 deletions

diff --git a/ast/ffitstable.c b/ast/ffitstable.c
new file mode 100644
index 0000000..9f2db37
--- /dev/null
+++ b/ast/ffitstable.c
@@ -0,0 +1,234 @@
+/*
+*+
+*  Name:
+*     ffitstable.c
+
+*  Purpose:
+*     Define a FORTRAN 77 interface to the AST FitsTable 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 FitsTable class.
+
+*  Routines Defined:
+*     AST_ISAFITSTABLE
+*     AST_FITSTABLE
+
+*  Copyright:
+*     Copyright (C) 2010 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:
+*     25-NOV-2010 (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 "ast_err.h"             /* AST error codes */
+#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 "fitstable.h"           /* C interface to the FitsTable class */
+
+F77_LOGICAL_FUNCTION(ast_isafitstable)( INTEGER(THIS),
+                                        INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   F77_LOGICAL_TYPE(RESULT);
+
+   astWatchSTATUS(
+   astAt( "AST_ISAFITSTABLE", NULL, 0 );
+      RESULT = astIsAFitsTable( astI2P( *THIS ) ) ? F77_TRUE : F77_FALSE;
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_fitstable)( INTEGER(HEADER),
+                                     CHARACTER(OPTIONS),
+                                     INTEGER(STATUS)
+                                     TRAIL(OPTIONS) ) {
+   GENPTR_INTEGER(HEADER)
+   GENPTR_CHARACTER(OPTIONS)
+   F77_INTEGER_TYPE(RESULT);
+   int i;
+   char *options;
+
+   astAt( "AST_FITSTABLE", 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( astFitsTable( astI2P( *HEADER ), "%s", options ) );
+      astFree( options );
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_gettableheader)( INTEGER(THIS),
+                                          INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   F77_INTEGER_TYPE(RESULT);
+
+   astAt( "AST_GETTABLEHEADER", NULL, 0 );
+   astWatchSTATUS(
+      RESULT = astP2I( astGetTableHeader( astI2P( *THIS ) ) );
+   )
+   return RESULT;
+}
+
+
+F77_SUBROUTINE(ast_puttableheader)( INTEGER(THIS),
+                                    INTEGER(HEADER),
+                                    INTEGER(STATUS) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_INTEGER(HEADER)
+
+   astAt( "AST_PUTTABLEHEADER", NULL, 0 );
+   astWatchSTATUS(
+      astPutTableHeader( astI2P( *THIS ), astI2P( *HEADER ) );
+   )
+}
+
+F77_INTEGER_FUNCTION(ast_columnnull)( INTEGER(THIS),
+                                      CHARACTER(COLUMN),
+                                      LOGICAL(SET),
+                                      INTEGER(NEWVAL),
+                                      LOGICAL(WASSET),
+                                      LOGICAL(HASNULL),
+                                      INTEGER(STATUS)
+                                      TRAIL(COLUMN) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_CHARACTER(COLUMN)
+   GENPTR_LOGICAL(SET)
+   GENPTR_INTEGER(NEWVAL)
+   GENPTR_LOGICAL(WASSET)
+   GENPTR_LOGICAL(HASNULL)
+   F77_INTEGER_TYPE(RESULT);
+   int wasset, hasnull;
+   char *column;
+
+   astAt( "AST_COLUMNNULL", NULL, 0 );
+   astWatchSTATUS(
+      column = astString( COLUMN, COLUMN_length );
+      RESULT = astColumnNull( astI2P( *THIS ), column,
+                              F77_ISTRUE( *SET ) ? 1 : 0, *NEWVAL,
+                              &wasset, &hasnull );
+      *WASSET = wasset ? F77_TRUE : F77_FALSE;
+      *HASNULL = hasnull ? F77_TRUE : F77_FALSE;
+      astFree( column );
+   )
+   return RESULT;
+}
+
+F77_INTEGER_FUNCTION(ast_columnsize)( INTEGER(THIS),
+                                      CHARACTER(COLUMN),
+                                      INTEGER(STATUS)
+                                      TRAIL(COLUMN) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_CHARACTER(COLUMN)
+   F77_INTEGER_TYPE(RESULT);
+   char *column;
+   size_t result;
+
+   astAt( "AST_COLUMNSIZE", NULL, 0 );
+   astWatchSTATUS(
+      column = astString( COLUMN, COLUMN_length );
+      result = astColumnSize( astI2P( *THIS ), column );
+      astFree( column );
+
+      RESULT = result;
+      if( (size_t) RESULT != result && astOK ) {
+         astError( AST__BIGTAB, "AST_COLUMNSIZE(FitsTable): The number of bytes in the "
+                   "column is too large to fit in a Fortran INTEGER.", status );
+      }
+
+   )
+   return RESULT;
+}
+
+F77_SUBROUTINE(ast_getcolumndata)( INTEGER(THIS),
+                                   CHARACTER(COLUMN),
+                                   REAL(RNULL),
+                                   DOUBLE(DNULL),
+                                   INTEGER(MXSIZE),
+                                   BYTE_ARRAY(COLDATA),
+                                   INTEGER(NELEM),
+                                   INTEGER(STATUS)
+                                   TRAIL(COLUMN) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_CHARACTER(COLUMN)
+   GENPTR_REAL(RNULL)
+   GENPTR_DOUBLE(DNULL)
+   GENPTR_INTEGER(MXSIZE)
+   GENPTR_BYTE_ARRAY(COLDATA)
+   GENPTR_INTEGER(NELEM)
+   char *column;
+
+   astAt( "AST_GETCOLUMNDATA", NULL, 0 );
+   astWatchSTATUS(
+      column = astString( COLUMN, COLUMN_length );
+      astGetColumnData( astI2P( *THIS ), column, *RNULL, *DNULL, *MXSIZE,
+                        (void *) COLDATA, NELEM );
+      astFree( column );
+   )
+}
+
+F77_SUBROUTINE(ast_putcolumndata)( INTEGER(THIS),
+                                   CHARACTER(COLUMN),
+                                   INTEGER(CLEN),
+                                   INTEGER(SIZE),
+                                   BYTE_ARRAY(COLDATA),
+                                   INTEGER(STATUS)
+                                   TRAIL(COLUMN) ) {
+   GENPTR_INTEGER(THIS)
+   GENPTR_CHARACTER(COLUMN)
+   GENPTR_INTEGER(CLEN)
+   GENPTR_INTEGER(SIZE)
+   GENPTR_BYTE_ARRAY(COLDATA)
+   char *column;
+
+   astAt( "AST_PUTCOLUMNDATA", NULL, 0 );
+   astWatchSTATUS(
+      column = astString( COLUMN, COLUMN_length );
+      astPutColumnData( astI2P( *THIS ), column, *CLEN, *SIZE, (void *) COLDATA );
+      astFree( column );
+   )
+}
+