summaryrefslogtreecommitdiffstats
path: root/ast/devtools/simtest.c
diff options
context:
space:
mode:
Diffstat (limited to 'ast/devtools/simtest.c')
-rw-r--r--ast/devtools/simtest.c864
1 files changed, 0 insertions, 864 deletions
diff --git a/ast/devtools/simtest.c b/ast/devtools/simtest.c
deleted file mode 100644
index e36634c..0000000
--- a/ast/devtools/simtest.c
+++ /dev/null
@@ -1,864 +0,0 @@
-#define astCLASS simtest
-
-#include <stdio.h>
-#include <string.h>
-#include "globals.h"
-#include "object.h"
-#include "fitschan.h"
-#include "mapping.h"
-#include "cmpmap.h"
-#include "error.h"
-#include "winmap.h"
-#include "tranmap.h"
-#include "frameset.h"
-#include "pointset.h"
-#include "f77.h"
-#include "memory.h"
-#include "switchmap.h"
-#include "unitmap.h"
-
-void DisplayList( AstMapping *, int, int, int * );
-void sink( const char * );
-const char *source( void );
-void tran( AstMapping *, int * );
-void diff( AstMapping *, int * );
-void tranm( AstMapping **, int *, int, int, int * );
-
-extern F77_SUBROUTINE(kpg1_asreg)( INTEGER(status) );
-
-
-FILE *fd=NULL;
-
-main(){
- char file[ 50 ], text[ 100 ];
- const char *class, *enc0;
- char *attrs, *kvalue;
- AstWinMap *winmap;
- AstFitsChan *fc, *fc2;
- AstChannel *ch, *channel;
- AstMapping *map, *map2, *tmap;
- AstObject *obj, *obj2;
- int opt, nin, nout, i;
- DECLARE_INTEGER(lstatus);
- int *status;
-
- lstatus = 0;
- status = &lstatus;
-
-/* F77_CALL(kpg1_asreg)( INTEGER_ARG(&lstatus) ); */
-
- printf( "Input FITS or AST file [| fitschan attrs] > ");
- scanf( "%s", file );
-
- attrs = strchr( file, '|' );
- if( attrs ) {
- *attrs = 0;
- attrs++;
- } else {
- attrs = "Carlin=0";
- }
-
- fd = fopen( file, "r" );
- if( fd ){
- ch = (AstChannel *) astFitsChan( NULL, sink, attrs, status );
- if( astOK ) {
- fc = (AstFitsChan *) ch;
- while( astOK && fgets( text, 100, fd ) ){
- astPutFits( fc, text, 0 );
- }
- astClearCard( fc );
-
- enc0 = astGetC( fc, "encoding" );
- printf("Using %s encoding\n", enc0 );
- obj = astRead( fc );
- }
- if( !obj || !astOK ){
- printf("Attempting to read file %s as an AST file...\n", file );
- astClearStatus;
- ch = astChannel( source, NULL, "", status );
-
- fclose( fd );
- fd = NULL;
- fd = fopen( file, "r" );
- obj = astRead( ch );
- if( obj != NULL ) {
- fc = NULL;
- } else {
- astClearStatus;
- printf("Reverting to FitsChan and using a dummy UnitMap. \n", file );
- obj = (AstObject *) astUnitMap( 1, "", status );
- astAnnul( ch );
- ch = (AstChannel *) fc;
- }
-
- }
- fclose( fd );
- fd = NULL;
- } else {
- astError( 1, "Cannot open file '%s'.", status, file );
- }
-
- class = obj ? astGetClass( obj ) : NULL;
- if( obj && class ){
- printf( "%s %s has been read from file '%s'.\n",
- strspn( class, "AEIOUaeiou" )?"An":"A", class, file );
- }
-
- if( obj && astIsAFitsChan( obj ) ) {
- fc = (AstFitsChan *) obj;
- astClearCard( fc );
- enc0 = astGetC( fc, "encoding" );
- printf("Using %s encoding\n", enc0 );
- obj = astRead( fc );
- }
-
- if( obj && astIsAMapping( obj ) ){
- map = (AstMapping *) obj;
-
- while( astOK ){
- printf("\n"
- "0 - Quit\n"
- "1 - Transform points\n"
- "2 - Simplify\n"
- "3 - Display\n"
- "4 - Write\n"
- "5 - Shift and scale\n"
- "6 - Set mapping attributes\n"
- "7 - Get mapping attributes\n"
- "8 - Dump the FitsChan\n"
- "9 - Differentiate\n"
- "10 - Find a FITS keyword\n"
- "\n"
- "Option > " );
- opt = 0;
- scanf( "%d", &opt );
-
- if( opt == 0 ){
- break;
-
- } else if( opt == 1 ){
- tran( map, status );
-
- } else if( opt == 2 ){
- map2 = astSimplify( map );
- map = astAnnul( map );
- map = map2;
- obj = (AstObject *) map2;
-
- } else if( opt == 3 ){
- if( astIsAFrameSet( map ) ) {
- tmap = astGetMapping( map, AST__BASE, AST__CURRENT );
- } else {
- tmap = astClone( map );
- }
- DisplayList( tmap, 1, astGetInvert( tmap ), status );
- tmap = astAnnul( tmap );
-
- } else if( opt == 4 ){
- printf( "\n Output FITS file > ");
- scanf( "%s", file );
- fd = fopen( file, "w" );
- if( fd ){
- if( fc ) {
- printf(" Encoding (orig=%s)> ",enc0 );
- scanf( "%s", text );
- if( !strcmp( enc0, text ) && fc != NULL ){
- printf("Using existing FitsChan.\n");
- fc2 = astClone( fc );
- fc = astAnnul( fc );
- ch = NULL;
- astClearCard( fc2 );
- astSet( fc2, "encoding=%s", status, text );
- } else {
- printf("Creating new FitsChan.\n");
- fc2 = astFitsChan( NULL, sink, "encoding=%s", status, text );
- astPutFits( fc2, "NAXIS1 = 1", 1 );
- astPutFits( fc2, "NAXIS2 = 1", 1 );
- astPutFits( fc2, "NAXIS3 = 1", 1 );
- astPutFits( fc2, "NAXIS4 = 1", 1 );
- }
- } else {
- printf(" Encoding for new FitsChan > " );
- scanf( "%s", text );
- fc2 = astFitsChan( NULL, sink, "encoding=%s", status, text );
- astPutFits( fc2, "NAXIS1 = 1", 1 );
- astPutFits( fc2, "NAXIS2 = 1", 1 );
- astPutFits( fc2, "NAXIS3 = 1", 1 );
- astPutFits( fc2, "NAXIS4 = 1", 1 );
- }
- if( !astWrite( fc2, obj ) ){
- printf("No object written to file '%s'.\n", file );
- }
- fc2 = astAnnul( fc2 );
- fclose( fd );
- fd = NULL;
- } else {
- astError( 1, "Cannot open file '%s'.", status, file );
- }
-
- } else if( opt == 5 ){
- nin = astGetNin( map );
- winmap = astWinMap( nin, NULL, NULL, NULL, NULL, "", status );
-
- printf( "\n Enter %d axis shift values > ", nin );
- for( i = 0; i < nin; i++ ) scanf( "%lg", (winmap->a) + i );
-
- printf( "\n Enter %d axis scale values > ", nin );
- for( i = 0; i < nin; i++ ) scanf( "%lg", (winmap->b) + i );
-
- if( !strcmp( class, "FrameSet" ) ){
- astInvert( winmap );
- astRemapFrame( obj, AST__BASE, winmap );
- } else {
- map2 = (AstMapping *) astCmpMap( winmap, map, 1, "", status );
- map = astAnnul( map );
- map = map2;
- obj = (AstObject *) map2;
- }
- winmap = astAnnul( winmap );
-
- } else if( opt == 6 ) {
- printf("Enter attribute specifications: ");
- scanf( "%s", text );
- astSet( map, text, status );
- if( !astOK ) astClearStatus;
-
- } else if( opt == 7 ) {
- printf("Enter attribute name: ");
- scanf( "%s", text );
- printf( "%s = %s\n", text, astGetC( map, text ) );
- if( !astOK ) astClearStatus;
-
- } else if( opt == 8 ) {
- printf( "\n Output text file > ");
- scanf( "%s", file );
- fd = fopen( file, "w" );
- if( fd ){
- channel = astChannel( NULL, sink, "Full=1", status );
- astDump( ch, channel );
- channel = astAnnul( channel );
- fclose( fd );
- fd = NULL;
- } else {
- printf("Failed to open output file '%s'.\n", file );
- perror(" ");
- }
-
- } else if( opt == 9 ){
- diff( map, status );
-
- } else if( opt == 10 ){
- if( fc ) {
- printf("Enter FITS keyword: ");
- scanf( "%s", text );
- if( astGetFitsS( fc, text, &kvalue ) ) {
- printf("%s = '%s'\n", text, kvalue );
- } else {
- printf("%s not found\n", text );
- }
- } else {
- printf("No FITS keywords are available\n", text );
- }
- }
- }
-
- } else {
- astError( 1, "Object should be a Mapping).", status );
- }
-
- if( ch ) ch = astAnnul( ch );
- obj = astAnnul( obj );
-
-}
-
-void DisplayList( AstMapping *map0, int series, int invert, int *status ){
- int nmap, i, inv, res, junk;
- AstMapping *mi, *mf, **map_list, *m, *m2, *map, *smap, *tm1, *tm2;
- AstCmpMap *cm;
- AstTranMap *tm;
- AstChannel *channel;
- int invi, invf, *invert_list, opt, old_inv, nin, in[20], *out, nout;
- int isswitch, tfwd, tinv;
- const char *class;
- char file[ 50 ];
- double dv;
-
- nmap = 0;
- map_list = NULL;
- invert_list = NULL;
- map = astClone( map0 );
-
- isswitch = astIsASwitchMap( map );
- if( isswitch ) {
- astSwitchList( (AstSwitchMap *) map, invert, &nmap, &map_list,
- &invert_list, status );
- } else {
- astMapList( map, series, invert, &nmap, &map_list, &invert_list );
- }
-
- while( 1 ){
- if( isswitch ) {
- printf("\nSwitchMap: (fwd selector, inv selector, routes )\n" );
- printf(" -3 - <transform points>\n");
- } else {
- printf("\nMappings in %s:\n", series?"series":"parallel" );
- printf(" -8 - <invert>\n");
- printf(" -7 - <write list to text file>\n");
- printf(" -6 - <split>\n");
- printf(" -5 - <simplify whole list>\n");
- printf(" -4 - <merge a single mapping>\n");
- printf(" -3 - <transform points>\n");
- printf(" -2 - <show to text file>\n");
- printf(" -1 - <show to screen>\n");
- }
- printf(" 0 - <up>\n");
- for( i = 0; i < nmap; i++ ){
- m = map_list[ i ];
- class = m ? astGetClass( m ) :"(NULL)";
- printf( " %d - %s\n", i + 1, class );
- }
- printf("\n Option > ", 0 );
- opt = 0;
- scanf("%d", &opt );
-
- if ( opt == -8 ) {
-
- for( i = 0; i < nmap; i++ ) map_list[ i ] = astAnnul( map_list[ i ] );
- map_list = astFree( map_list );
- invert_list = astFree( invert_list );
- nmap = 0;
-
- invert = ( invert == 0 );
- if( isswitch ) {
- astSwitchList( (AstSwitchMap *) map, invert, &nmap,
- &map_list, &invert_list, status );
- } else {
- astMapList( map, series, invert, &nmap, &map_list, &invert_list );
- }
-
- } else if ( opt == -7 && !isswitch ) {
- printf( "\n Output text file > ");
- scanf( "%s", file );
- fd = fopen( file, "w" );
- if( fd ){
- channel = astChannel( NULL, sink, "", status );
- for( i = 0; i < nmap; i++ ) {
- inv = astGetInvert( map_list[ i ] );
- astSetInvert( map_list[ i ], invert_list[ i ] );
- astDump( map_list[ i ], channel );
- astSetInvert( map_list[ i ], inv );
- }
- channel = astAnnul( channel );
- fclose( fd );
- fd = NULL;
- } else {
- printf("Failed to open output file.\n" );
- }
-
- } else if( opt == -6 && !isswitch ) {
- old_inv = astGetInvert( map );
- astSetInvert( map, invert );
-
- printf( "\n Pick how many input axes (max=%d)? > ",
- astGetNin( map ) );
- scanf( "%lg", &dv );
- nin = dv;
-
- printf( "\n Enter %d zero-based axis indices > ", nin );
- for( i = 0; i < nin; i++ ) {
- scanf( "%lg", &dv );
- in[ i ] = dv;
- }
-
- out = astMapSplit( map, nin, in, &smap );
- astSetInvert( map, old_inv );
-
- if( out ) {
- printf( "Corresponding outputs are: " );
- nout = astGetNout( smap );
- for( i = 0; i < nout; i++ ) printf( "%d ", out[ i ] );
- printf("\n");
- out = astFree( out );
- map = astAnnul( map );
- map = smap;
- } else {
- printf( "Selected inputs do not correspond to a subset of outputs\n");
- }
-
- for( i = 0; i < nmap; i++ ) map_list[ i ] = astAnnul( map_list[ i ] );
- map_list = astFree( map_list );
- invert_list = astFree( invert_list );
-
- nmap = 0;
- astMapList( map, series, invert, &nmap, &map_list, &invert_list );
-
- } else if( opt == -5 && !isswitch ) {
- old_inv = astGetInvert( map );
- astSetInvert( map, invert );
- smap = astSimplify( map );
- if( smap == map ) printf("\n No simplification performed\n" );
- astSetInvert( map, old_inv );
- map = astAnnul( map );
- map = smap;
-
- for( i = 0; i < nmap; i++ ) map_list[ i ] = astAnnul( map_list[ i ] );
- map_list = astFree( map_list );
- invert_list = astFree( invert_list );
-
- nmap = 0;
- astMapList( map, series, invert, &nmap, &map_list, &invert_list );
-
- } else if( opt == -4 && !isswitch ) {
- printf("\n Mapping index > ", 0 );
- scanf("%d", &i );
- i--;
- res = astMapMerge( map_list[ i ], i, series, &nmap,
- &map_list, &invert_list );
- if( res != -1 ) {
- printf( "Index of first modified Mapping: %d\n", res + 1 );
- } else {
- printf( "No mappings modified\n" );
- }
-
- } else if( opt == -3 && series && !isswitch ) {
-
- if( invert == astGetInvert( map ) ) {
- tfwd = astGetTranForward( map );
- tinv = astGetTranInverse( map );
- } else {
- tinv = astGetTranForward( map );
- tfwd = astGetTranInverse( map );
- }
-
- if( !tfwd ) {
- printf( "\n No forward transformation available\n" );
-
- } else if( !tinv ) {
- printf( "\n No inverse transformation available\n" );
- tranm( map_list, invert_list, nmap, 0, status );
-
- } else {
- tranm( map_list, invert_list, nmap, 1, status );
-
- }
-
- } else if( opt == -3 && isswitch ) {
-
- if( invert == astGetInvert( map ) ) {
- tfwd = astGetTranForward( map );
- tinv = astGetTranInverse( map );
- } else {
- tinv = astGetTranForward( map );
- tfwd = astGetTranInverse( map );
- }
-
- if( !tfwd ) {
- printf( "\n No forward transformation available\n" );
-
- } else if( !tinv ) {
- printf( "\n No inverse transformation available\n" );
- tranm( &map, &invert, 1, 0, status );
-
- } else {
- tranm( &map, &invert, 1, 1, status );
-
- }
-
- } else if( opt == -3 && !series && !isswitch ) {
- old_inv = astGetInvert( map );
- astSetInvert( map, invert );
- tran( map, status );
- astSetInvert( map, old_inv );
-
- } else if ( opt == -2 && !isswitch ) {
-
- for( i = 0; i < nmap; i++ ) {
- inv = astGetInvert( map_list[ i ] );
- astSetInvert( map_list[ i ], invert_list[ i ] );
- if( i == 0 ) {
- tm1 = astClone( map_list[ 0 ] );
- } else {
- tm2 = (AstMapping *) astCmpMap( tm1, map_list[ i ], series, "", status );
- tm1 = astAnnul( tm1 );
- tm1 = tm2;
- }
- astSetInvert( map_list[ i ], inv );
- }
-
- printf( "\n Output text file > ");
- scanf( "%s", file );
- fd = fopen( file, "w" );
- if( fd ){
- channel = astChannel( NULL, sink, "", status );
- old_inv = astGetInvert( tm1 );
- astSetInvert( tm1, invert );
- astDump( tm1, channel );
- astSetInvert( tm1, old_inv );
- channel = astAnnul( channel );
- fclose( fd );
- fd = NULL;
- } else {
- printf("Failed to open output file.\n" );
- }
- tm1= astAnnul( tm1 );
-
- } else if ( opt == -1 && !isswitch ) {
-
-
- for( i = 0; i < nmap; i++ ) {
- inv = astGetInvert( map_list[ i ] );
- astSetInvert( map_list[ i ], invert_list[ i ] );
- if( i == 0 ) {
- tm1 = astClone( map_list[ 0 ] );
- } else {
- tm2 = (AstMapping *) astCmpMap( tm1, map_list[ i ],
- series, "", status );
- tm1 = astAnnul( tm1 );
- tm1 = tm2;
- }
- astSetInvert( map_list[ i ], inv );
- }
-
- channel = astChannel( NULL, NULL, "Full=1", status );
- old_inv = astGetInvert( tm1);
- astSetInvert( tm1, invert );
- astDump( tm1, channel );
- astSetInvert( tm1, old_inv );
- channel = astAnnul( channel );
- tm1= astAnnul( tm1 );
-
- } else if( opt > 0 && opt <= nmap && map_list[ opt - 1 ] ){
- m2 = map_list[ opt - 1 ];
- inv = invert_list[ opt - 1 ];
- class = astGetClass( m2 );
-
- if( m2 && astIsAFrameSet( m2 ) ){
- m = astGetMapping( m2, AST__BASE, AST__CURRENT );
- class = astGetClass( m );
-
- } else if( astIsATranMap( m2 ) ){
-
- astDecompose( m2, &mf, &mi, &invf, &invi, &junk );
-
- if( inv ) {
- m = mi;
- inv = invi;
- } else {
- m = mf;
- inv = invf;
- }
- class = astGetClass( m );
-
- } else {
- m = astClone( m2 );
- }
-
- if( strcmp( class, "CmpMap" ) ){
- DisplayList( m, 1, inv, status );
- } else {
- cm = (AstCmpMap *) m;
- DisplayList( m, cm->series, inv, status );
- }
- } else {
- break;
- }
- }
-
- for( i = 0; i < nmap; i++ ) map_list[ i ] = astAnnul( map_list[ i ] );
- map_list = astFree( map_list );
- invert_list = astFree( invert_list );
- nmap = 0;
-
-}
-
-
-void sink( const char *text ){
- if( fd && text ) fprintf( fd, "%s\n", text );
-}
-
-const char *source( void ){
- static char text[2000];
- static char *ptext;
- if( fd ) {
- ptext = fgets( text, 2000, fd );
- } else {
- ptext = NULL;
- }
- return ptext;
-}
-
-void tran( AstMapping *map, int *status ){
- char text[2000];
- AstFrame *infrm, *outfrm;
- AstPointSet *pset1, *pset2;
- double **ptr1, **ptr2;
- int nin, nout, i, ok;
-
- if( !astGetTranForward( map ) ) {
- printf( "\n No forward transformation available\n" );
- return;
- }
-
- if( !astGetTranInverse( map ) ) {
- printf( "\n No inverse transformation available\n" );
- }
-
- nin = astGetNin( map );
- nout = astGetNout( map );
- pset1 = astPointSet( 1, nin, "", status );
- ptr1 = astGetPoints( pset1 );
- pset2 = astPointSet( 1, nout, "", status );
- ptr2 = astGetPoints( pset2 );
-
- if( astIsAFrameSet( map ) ) {
- infrm = astGetFrame( map, AST__BASE );
- outfrm = astGetFrame( map, AST__CURRENT );
- } else if( astIsAFrame( map ) ) {
- infrm = (AstFrame *) map;
- outfrm = (AstFrame *) map;
- } else {
- infrm = NULL;
- outfrm = NULL;
- }
-
- if( astOK ){
- ok = 0;
- if( infrm != NULL ) {
- ok = 1;
- for( i = 0; i < nin; i++ ) {
- printf( "\n Enter %s value (or <bad>) > ", astGetLabel( infrm, i ) );
- scanf( "%s", text );
- if( !astUnformat( infrm, i, text, ptr1[ i ] ) ) {
- ok = 0;
- printf( "Bad formatted value supplied! Now try unformatted axis values...\n");
- break;
- }
- }
- }
-
- if( !ok ) {
- printf( "\n Enter %d axis values (BAD=999) > ", nin );
- for( i = 0; i < nin; i++ ) {
- scanf( "%lg", ptr1[ i ] );
- if( ptr1[ i ][ 0 ] == 999.0 ) ptr1[ i ][ 0 ] = AST__BAD;
- }
- }
-
- astTransform( map, pset1, 1, pset2 );
- printf( "\n Forward: [ " );
- for( i = 0; i < nin; i++ ) {
- if( ptr1[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr1[ i ][ 0 ] );
- }
- }
-
- printf( "] -> [ " );
- for( i = 0; i < nout; i++ ) {
- if( ptr2[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr2[ i ][ 0 ] );
- }
- }
- printf( "]\n" );
-
- if( infrm != NULL && outfrm != NULL ) {
- printf( " ( " );
- for( i = 0; i < nin; i++ ) {
- printf( "%s ", astFormat( infrm, i, ptr1[ i ][ 0 ] ) );
- }
- printf( ") -> ( " );
- for( i = 0; i < nout; i++ ) {
- printf( "%s ", astFormat( outfrm, i, ptr2[ i ][ 0 ] ) );
- }
- printf( ") formatted\n" );
- }
-
-
- if( astGetTranInverse( map ) ) {
- astTransform( map, pset2, 0, pset1 );
- printf( "\n Inverse: [ " );
- for( i = 0; i < nout; i++ ) {
- if( ptr2[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr2[ i ][ 0 ] );
- }
- }
-
- printf( "] -> [ " );
- for( i = 0; i < nin; i++ ) {
- if( ptr1[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr1[ i ][ 0 ] );
- }
- }
- printf( "]\n" );
-
-
- if( infrm != NULL && outfrm != NULL ) {
- printf( " ( " );
- for( i = 0; i < nout; i++ ) {
- printf( "%s ", astFormat( outfrm, i, ptr2[ i ][ 0 ] ) );
- }
- printf( ") -> ( " );
- for( i = 0; i < nin; i++ ) {
- printf( "%s ", astFormat( infrm, i, ptr1[ i ][ 0 ] ) );
- }
- printf( ") formatted\n" );
- }
- }
- }
- pset1 = astAnnul( pset1 );
- pset2 = astAnnul( pset2 );
-}
-
-void tranm( AstMapping **map, int *inv_list, int nmap, int inv, int *status ){
- AstPointSet *pset1, *pset2;
- double **ptr1, **ptr2;
- int nin, nout, i, j, inv0;
-
- inv0 = astGetInvert( map[ 0 ] );
- astSetInvert( map[ 0 ], inv_list[ 0 ] );
- nin = astGetNin( map[ 0 ] );
- astSetInvert( map[ 0 ], inv0 );
-
- pset1 = astPointSet( 1, nin, "", status );
- ptr1 = astGetPoints( pset1 );
-
- if( astOK ){
- printf( "\n Enter %d axis values (BAD=999) > ", nin );
- for( i = 0; i < nin; i++ ) {
- scanf( "%lg", ptr1[ i ] );
- if( ptr1[ i ][ 0 ] == 999.0 ) ptr1[ i ][ 0 ] = AST__BAD;
- }
-
- printf("Forward...\n");
- for( j=0; j<nmap; j++ ) {
- inv0 = astGetInvert( map[ j ] );
- astSetInvert( map[ j ], inv_list[ j ] );
-
- nout = astGetNout( map[ j ] );
- pset2 = astPointSet( 1, nout, "", status );
- ptr2 = astGetPoints( pset2 );
-
- astTransform( map[ j ], pset1, 1, pset2 );
-
- printf( " Mapping %d gave [ ", j+1 );
- for( i = 0; i < nout; i++ ) {
- if( ptr2[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr2[ i ][ 0 ] );
- }
- }
-
- printf( "]\n" );
-
- pset1 = astAnnul( pset1 );
- astSetInvert( map[ j ], inv0 );
-
- pset1 = pset2;
- ptr1 = ptr2;
- nin = nout;
- }
-
- printf("Inverse...\n");
-
- for( j=nmap - 1; j>=0 && inv; j-- ) {
- inv0 = astGetInvert( map[ j ] );
- astSetInvert( map[ j ], inv_list[ j ] );
-
- nin = astGetNin( map[ j ] );
- pset2 = astPointSet( 1, nin, "", status );
- ptr2 = astGetPoints( pset2 );
-
- astTransform( map[ j ], pset1, 0, pset2 );
-
- printf( " Mapping %d gave [ ", j+1 );
- for( i = 0; i < nin; i++ ) {
- if( ptr2[ i ][ 0 ] == AST__BAD ) {
- printf( "999.0 " );
- } else {
- printf( "%.*g ", DBL_DIG, ptr2[ i ][ 0 ] );
- }
- }
- printf( "]\n" );
-
- pset1 = astAnnul( pset1 );
- astSetInvert( map[ j ], inv0 );
-
- pset1 = pset2;
- ptr1 = ptr2;
- nout = nin;
- }
-
- pset1 = astAnnul( pset1 );
-
- }
-}
-
-
-
-
-void diff( AstMapping *map, int *status ){
- char text[200];
- AstFrame *infrm, *outfrm;
- int nin, nout, i, ok, ax1, ax2;
- double at[2000];
-
-
- nin = astGetNin( map );
- nout = astGetNout( map );
-
- if( astIsAFrameSet( map ) ) {
- infrm = astGetFrame( map, AST__BASE );
- outfrm = astGetFrame( map, AST__CURRENT );
- } else if( astIsAFrame( map ) ) {
- infrm = (AstFrame *) map;
- outfrm = (AstFrame *) map;
- } else {
- infrm = NULL;
- outfrm = NULL;
- }
-
- if( astOK ){
- ok = 0;
- if( infrm != NULL ) {
- ok = 1;
- for( i = 0; i < nin; i++ ) {
- printf( "\n Enter %s value (or <bad>) > ", astGetLabel( infrm, i ) );
- scanf( "%s", text );
- if( !astUnformat( infrm, i, text, at + i ) ) {
- ok = 0;
- printf( "Bad formatted value supplied! Now try unformatted axis values...\n");
- break;
- }
- }
- }
-
- if( !ok ) {
- printf( "\n Enter %d axis values (BAD=999) > ", nin );
- for( i = 0; i < nin; i++ ) {
- scanf( "%lg", at + i );
- if( at[ i ] == 999.0 ) at[ i ] = AST__BAD;
- printf("(got %.*g)\n", DBL_DIG, at[ i ] );
- }
-
- }
-
- for( i = 0; i < nin; i++ ) {
- printf("(%d: got %.*g)\n", i, DBL_DIG, at[ i ] );
- }
-
- printf( "\n Enter output to differentiate (zero-based index): > " );
- scanf( "%d", &ax1 );
- printf( "\n Enter input to vary (zero-based index): > " );
- scanf( "%d", &ax2 );
- printf( "%.*g \n", DBL_DIG, astRate( map, at, ax1, ax2 ) );
- }
-}
-
generated by cgit v0.12 at 2024-09-30 00:14:40 (GMT)