summaryrefslogtreecommitdiffstats
path: root/ast/ffitstable.c
diff options
context:
space:
mode:
Diffstat (limited to 'ast/ffitstable.c')
-rw-r--r--ast/ffitstable.c234
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 );
+ )
+}
+