Ejemplo n.º 1
0
  void hdf_sd_getdata_pro( EnvT* e)
  {
    SizeT nParam=e->NParam(2);
    char fieldname[256];
    int32 rank;
    int32 dims[MAXRANK];
    int32 start[MAXRANK];
    int32 edges[MAXRANK];
    int32 stride[MAXRANK];
    int32 dtype;
    int32 nattrs;
    int32 status;

    DLong sds_id;
    e->AssureScalarPar<DLongGDL>( 0, sds_id);

    DLongGDL* startKW = e->IfDefGetKWAs<DLongGDL>( 0);
    DLongGDL* strideKW = e->IfDefGetKWAs<DLongGDL>( 1);
    DLongGDL* countKW = e->IfDefGetKWAs<DLongGDL>( 2);

    status = SDgetinfo(sds_id, fieldname, &rank, dims, &dtype, &nattrs);

    if (status == FAIL)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Invalid SD dataset ID: " + i2s(sds_id));

    for( SizeT i = 0; i < rank; i++) {
      start[i] = 0;
      stride[i] = 1;
      edges[i] = dims[i];
    }

    if (startKW != NULL) {
      if (startKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of START dimensions:"
			    +i2s(startKW->N_Elements()));
      memcpy(start, &(*startKW)[0], rank*sizeof(int32));

      // Switch to C-ordering
      for( SizeT i = 0; i < rank/2; i++) {
	int32 temp = start[i];
	start[i] = start[rank-1];
	start[rank-1] = temp;
      }
    }

    if (countKW != NULL) {
      if (countKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of COUNT dimensions:"
			    +i2s(countKW->N_Elements()));
      memcpy(edges, &(*countKW)[0], rank*sizeof(int32));

      // Switch to C-ordering
      for( SizeT i = 0; i < rank/2; i++) {
	int32 temp = edges[i];
	edges[i] = edges[rank-1];
	edges[rank-1] = temp;
      }
    }

    if (strideKW != NULL) {
      if (strideKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of STRIDE dimensions:"
			    +i2s(strideKW->N_Elements()));
      memcpy(stride, &(*strideKW)[0], rank*sizeof(int32));

      // Switch to C-ordering
      for( SizeT i = 0; i < rank/2; i++) {
	int32 temp = stride[i];
	stride[i] = stride[rank-1];
	stride[rank-1] = temp;
      }
    }

    // Reverse order of edges for setting up GDL array(
    for( SizeT i = 0; i < rank/2; i++) {
      int32 tempedge = edges[i];
      edges[i] = edges[rank-1];
      edges[rank-1] = tempedge;
    }

    dimension dim((DLong *) edges, rank);

    // Return to original edge values for HDF read
    for( SizeT i = 0; i < rank/2; i++) {
      int32 tempedge = edges[i];
      edges[i] = edges[rank-1];
      edges[rank-1] = tempedge;
    }

    switch ( dtype) {

       case DFNT_FLOAT64: {
	 hdf_sd_getdata_template< DDoubleGDL>( e, dim, sds_id,
					       start, edges, stride, strideKW);
	 break;
       }
       case DFNT_FLOAT32: {
	 hdf_sd_getdata_template< DFloatGDL>( e, dim, sds_id,
					      start, edges, stride, strideKW);
	 break;
       }
       case DFNT_UINT32: {
	 hdf_sd_getdata_template< DULongGDL>( e, dim, sds_id,
					      start, edges, stride, strideKW);
	 break;
       }
       case DFNT_INT32: {
	 hdf_sd_getdata_template< DLongGDL>( e, dim, sds_id,
					     start, edges, stride, strideKW);
	 break;
       }
       case DFNT_UINT16: {
	 hdf_sd_getdata_template< DUIntGDL>( e, dim, sds_id,
					     start, edges, stride, strideKW);
	 break;
       }
       case DFNT_INT16: {
	 hdf_sd_getdata_template< DIntGDL>( e, dim, sds_id,
					    start, edges, stride, strideKW);
	 break;
       }
      case DFNT_UINT8:
      case DFNT_INT8: {
	 hdf_sd_getdata_template< DByteGDL>( e, dim, sds_id,
					     start, edges, stride, strideKW);
	 break;
       }
    }

  }
  void device( EnvT* e)
  {

    GraphicsDevice* actDevice = GraphicsDevice::GetDevice();
    //GET functions are examined BEFORE setting functions.
    //GET_CURRENT_FONT ? //TODO
    
    // GET_DECOMPOSED ?
    {
    static int get_decomposedIx = e->KeywordIx( "GET_DECOMPOSED" );
    if ( e->KeywordPresent( get_decomposedIx ) )
      {
        DLong value = actDevice->GetDecomposed( );
        if ( value == -1 )
          e->Throw( "Keyword GET_DECOMPOSED not allowed for call to: DEVICE" );
        else
          e->SetKW( get_decomposedIx, new DLongGDL( value ) );
      }
    }
    
    //GET_FONTNAMES? //TODO
    //GET_FONTNUM? //TODO
    
    // GET_GRAPHICS_FUNCTION
    {
    static int get_graphics_FunctionIx = e->KeywordIx( "GET_GRAPHICS_FUNCTION");
    if( e->KeywordPresent( get_graphics_FunctionIx)) 
      {
        DLong value = actDevice->GetGraphicsFunction();
        if(value == -1)
          e->Throw( "Keyword GET_GRAPHICS_FUNCTION not allowed for call to: DEVICE");
        else 
          e->SetKW( get_graphics_FunctionIx, new DLongGDL( value));
      }
    }
    
    // GET_PAGE_SIZE ?
    {
     static int get_page_sizeIx = e->KeywordIx("GET_PAGE_SIZE");
      if( e->KeywordPresent( get_page_sizeIx)) 
      {
       DIntGDL* value = actDevice->GetPageSize();
       if (value == NULL) 
          e->Throw( "Keyword GET_PAGE_SIZE not allowed for call to: DEVICE");
       else 
          e->SetKW( get_page_sizeIx, value );
      }
    }

    // GET_PIXEL_DEPTH ?
    {
      static int get_pixel_depthIx = e->KeywordIx( "GET_PIXEL_DEPTH");
      if( e->KeywordPresent( get_pixel_depthIx)) 
      {
        DLong value = actDevice->GetPixelDepth();
        if(value == -1)
          e->Throw( "Keyword GET_PIXEL_DEPTH not allowed for call to: DEVICE");
        else 
          e->SetKW( get_pixel_depthIx, new DLongGDL( value ));
      }
    }

    // GET_SCREEN_SIZE
    {
      static int get_screen_sizeIx = e->KeywordIx("GET_SCREEN_SIZE");
      if( e->KeywordPresent( get_screen_sizeIx)) 
      {
       DIntGDL* fvalue=actDevice->GetScreenSize();
       if (fvalue == NULL) 
          e->Throw( "Keyword GET_SCREEN_SIZE not allowed for call to: DEVICE");
       else {
         (*fvalue)[0]=floor((*fvalue)[0]);
         (*fvalue)[1]=floor((*fvalue)[1]);
         e->SetKW( get_screen_sizeIx, fvalue);
       }
      }
    }
   
    // GET_VISUAL_DEPTH ?
    {
      static int get_visual_depthIx = e->KeywordIx( "GET_VISUAL_DEPTH");
      if( e->KeywordPresent( get_visual_depthIx)) 
      {
        DLong value = actDevice->GetVisualDepth();
        if(value == -1)
          e->Throw( "Keyword GET_VISUAL_DEPTH not allowed for call to: DEVICE");
        else 
          e->SetKW( get_visual_depthIx, new DLongGDL( value));
      }
    }
    
    // GET_VISUAL_NAME ?
    {
      static int get_visual_nameIx = e->KeywordIx( "GET_VISUAL_NAME");
      if( e->KeywordPresent( get_visual_nameIx)) 
      {
        DString value = actDevice->GetVisualName();
        if(value == "")
          e->Throw( "Keyword GET_VISUAL_NAME not allowed for call to: DEVICE");
        else 
          e->SetKW( get_visual_nameIx, new DStringGDL( value));
      }
    } 

    // GET_WINDOW_POSITION ?
    {
     static int get_window_positionIx = e->KeywordIx("GET_WINDOW_POSITION");
      if( e->KeywordPresent( get_window_positionIx)) 
      {
       DIntGDL* value = actDevice->GetWindowPosition();
       if (value == NULL) 
          e->Throw( "Keyword GET_WINDOW_POSITION not allowed for call to: DEVICE");
       else 
          e->SetKW( get_window_positionIx, value);
      }
    }

    // GET_WRITE_MASK ? 
    {
    static int get_write_maskIx = e->KeywordIx( "GET_WRITE_MASK");
      if( e->KeywordPresent( get_write_maskIx)) 
      {
        DLong value = actDevice->GetWriteMask();
        if(value == -1)
          e->Throw( "Keyword GET_WRITE_MASK not allowed for call to: DEVICE");
        else 
          e->SetKW( get_write_maskIx, new DLongGDL( value));
      }
    }    

    // WINDOW_STATE ?
    {
    static int window_stateIx = e->KeywordIx( "WINDOW_STATE");
      if( e->KeywordPresent( window_stateIx)) 
      {
        DByteGDL* value = actDevice->WindowState();
        if(value == NULL)
          e->Throw( "Keyword WINDOW_STATE not allowed for call to: DEVICE");
        else 
          e->SetKW( window_stateIx, value);
      }
    } 

    // CLOSE_FILE {{{
    {
      static int closeFileIx = e->KeywordIx( "CLOSE_FILE"); 
      if( e->KeywordSet( closeFileIx))
	{
	  bool success = actDevice->CloseFile();
	  if( !success)
	    e->Throw( "Current device does not support keyword CLOSE_FILE.");
	}
    }

    // Z_BUFFERING 
    {
      static int z_bufferingIx = e->KeywordIx( "Z_BUFFERING"); 
      BaseGDL* z_buffering = e->GetKW( z_bufferingIx);
      if( z_buffering != NULL)
	{
	  bool success = actDevice->ZBuffering( e->KeywordSet( z_bufferingIx));
	  if( !success)
	    e->Throw( "Current device does not support keyword Z_BUFFERING.");
	}
    }

    // SET_RESOLUTION 
    {
      static int set_resolutionIx = e->KeywordIx( "SET_RESOLUTION"); 
      BaseGDL* set_resolution = e->GetKW( set_resolutionIx);
      if( set_resolution != NULL)
	{
	  DLongGDL* resolution = e->GetKWAs<DLongGDL>( set_resolutionIx);
	  if( resolution->N_Elements() != 2)
	    e->Throw( "Keyword array parameter SET_RESOLUTION must have 2 elements.");
	  DLong x = (*resolution)[0];
	  DLong y = (*resolution)[1];

	  if( x<0 || y<0)
	    e->Throw( "Value of Resolution is out of allowed range.");

	  bool success = actDevice->SetResolution( x, y);
	  if( !success)
	    e->Throw( "Current device does not support keyword SET_RESOLUTION.");
	}
    }

    // DECOMPOSED 
    {
      static int decomposedIx = e->KeywordIx( "DECOMPOSED"); 
      BaseGDL* decomposed = e->GetKW( decomposedIx);
      if( decomposed != NULL)
	{
	  bool success = actDevice->Decomposed( e->KeywordSet( decomposedIx));
	  if( !success)
	    e->Throw( "Current device does not support keyword DECOMPOSED.");
	}
    }  

    // SET_GRAPHICS_FUNCTION
    {
      static int set_graphicsFunctionIx = e->KeywordIx( "SET_GRAPHICS_FUNCTION");
      BaseGDL* set_gfunction = e->GetKW( set_graphicsFunctionIx);
      if( set_gfunction != NULL)
	{
	  DLongGDL* gfunction = e->GetKWAs<DLongGDL>( set_graphicsFunctionIx);
	  bool success = actDevice->SetGraphicsFunction((*gfunction)[0]);
	  if( !success)
	    e->Throw( "Current device does not support keyword SET_GRAPHICS_FUNCTION.");
	}
    }
    // CURSOR_STANDARD
    {
      static int cursorStandardIx = e->KeywordIx( "CURSOR_STANDARD");
      BaseGDL* res = e->GetKW( cursorStandardIx);
      if( res != NULL)
	{
	  DLongGDL* val = e->GetKWAs<DLongGDL>( cursorStandardIx);
	  bool success = actDevice->CursorStandard((*val)[0]);
	  if( !success)
	    e->Throw( "Current device does not support keyword CURSOR_STANDARD.");
	}
    }
    // RETAIN
    {
      static int valIx = e->KeywordIx( "RETAIN");
      BaseGDL* res = e->GetKW( valIx);
      if( res != NULL)
	{
	  DLongGDL* val = e->GetKWAs<DLongGDL>( valIx);
	  bool success = actDevice->SetBackingStore((*val)[0]);
	  if( !success)
	    e->Throw( "Current device does not support keyword RETAIN.");
	}
    }
    // CURSOR_CROSSHAIR
    {
      static int valIx = e->KeywordIx( "CURSOR_CROSSHAIR");
      BaseGDL* res = e->GetKW( valIx);
      if( res != NULL)
	{
	  bool success = actDevice->CursorCrosshair();
	  if( !success)
	    e->Throw( "Current device does not support keyword CURSOR_CROSSHAIR.");
	}
    }
    // CURSOR_ORIGINAL (WARNING: SAME CODE AS  CURSOR_CROSSHAIR!)
    {
      static int valIx = e->KeywordIx( "CURSOR_ORIGINAL");
      BaseGDL* res = e->GetKW( valIx);
      if( res != NULL)
	{
	  bool success = actDevice->CursorCrosshair();
	  if( !success)
	    e->Throw( "Current device does not support keyword CURSOR_ORIGINAL.");
	}
    }

    // FILENAME 
    {
      static int fileNameIx = e->KeywordIx( "FILENAME"); 
      BaseGDL* fileName = e->GetKW( fileNameIx);
      if( fileName != NULL)
	{
	  DString fName;
	  e->AssureStringScalarKW( fileNameIx, fName);
	  if( fName == "")
	    e->Throw( "Null filename not allowed.");
	  WordExp(fName);
	  bool success = actDevice->SetFileName( fName);
	  if( !success)
	    e->Throw( "Current device does not support keyword FILENAME.");
	}
    }

    // LANDSCAPE and PORTRAIT need to be executed before XSIZE, YSIZE, XOFFSET and YOFFSET!
    {
      static int portraitIx = e->KeywordIx( "PORTRAIT");
      static int landscapeIx = e->KeywordIx( "LANDSCAPE"); 
    //IDL consider the value of the first typed of the two options, if both are present.
    //I dunoo how to do that with GDL parsing.
      if (e->KeywordSet(portraitIx) && e->KeywordSet(landscapeIx)) 
        Warning("Warning: both PORTRAIT and LANDSCAPE specified!");

      // LANDSCAPE 
      {
        BaseGDL* landscapeKW=e->GetKW(landscapeIx);
        if (landscapeKW != NULL)
        {
	      DLong isLandscape;
	      e->AssureLongScalarKW( landscapeIx, isLandscape);
          if (isLandscape == 0) {
            bool success = actDevice->SetPortrait();
            if (!success) e->Throw("Current device does not support keyword LANDSCAPE");
          } else {
            bool success = actDevice->SetLandscape();
            if (!success) e->Throw("Current device does not support keyword LANDSCAPE");
          }
        }
      }

      // PORTRAIT 
      {
        BaseGDL* portraitKW=e->GetKW(portraitIx);
        if (portraitKW != NULL)
        {
	      DLong isPortrait;
	      e->AssureLongScalarKW( portraitIx, isPortrait);
          if (isPortrait == 0) {
            bool success = actDevice->SetLandscape();
            if (!success) e->Throw("Current device does not support keyword PORTRAIT");
          } else {
            bool success = actDevice->SetPortrait();
            if (!success) e->Throw("Current device does not support keyword PORTRAIT");
          }
        }
      }
    }

    {
      static int inchesIx = e->KeywordIx( "INCHES");
      // XOFFSET 
      {
        static int xOffsetIx = e->KeywordIx( "XOFFSET");
        BaseGDL* xOffsetKW = e->GetKW( xOffsetIx);
        if( xOffsetKW != NULL)
	  {
	    DFloat xOffsetValue;
	    e->AssureFloatScalarKW( xOffsetIx, xOffsetValue);
	    bool success = actDevice->SetXOffset( xOffsetValue 
						  * (e->KeywordSet(inchesIx) ? 100. * GSL_CONST_MKSA_INCH : 1.)
						  );
	    if( !success)
	      e->Throw( "Current device does not support keyword XOFFSET.");
	  } 
      }

      // YOFFSET 
      {
        static int yOffsetIx = e->KeywordIx( "YOFFSET");
        BaseGDL* yOffsetKW = e->GetKW( yOffsetIx);
        if( yOffsetKW != NULL)
	  {
	    DFloat yOffsetValue;
	    e->AssureFloatScalarKW( yOffsetIx, yOffsetValue);
	    bool success = actDevice->SetYOffset( yOffsetValue 
						  * (e->KeywordSet(inchesIx) ? 100. * GSL_CONST_MKSA_INCH : 1.)
						  );
	    if( !success)
	      e->Throw( "Current device does not support keyword YOFFSET.");
	  } 
      }

      // XSIZE 
      {
        static int xSizeIx = e->KeywordIx( "XSIZE");
        BaseGDL* xSizeKW = e->GetKW( xSizeIx);
        if( xSizeKW != NULL)
	  {
	    DFloat xSizeValue;
	    e->AssureFloatScalarKW( xSizeIx, xSizeValue);
	    bool success = actDevice->SetXPageSize( xSizeValue 
						    * (e->KeywordSet(inchesIx) ? 100. * GSL_CONST_MKSA_INCH : 1.)
						    );
	    if( !success)
	      e->Throw( "Current device does not support keyword XSIZE.");
	  } 
      }

      // YSIZE 
      {
        static int ySizeIx = e->KeywordIx( "YSIZE");
        BaseGDL* ySizeKW = e->GetKW( ySizeIx);
        if( ySizeKW != NULL)
	  {
	    DFloat ySizeValue;
	    e->AssureFloatScalarKW( ySizeIx, ySizeValue);
	    bool success = actDevice->SetYPageSize( ySizeValue
						    * (e->KeywordSet(inchesIx) ? 100. * GSL_CONST_MKSA_INCH : 1.)
						    );
	    if( !success)
	      e->Throw( "Current device does not support keyword YSIZE.");
	  } 
      }
    }

    // SCALE_FACTOR 
    {
      static int scaleIx = e->KeywordIx( "SCALE_FACTOR");
      BaseGDL* scaleKW = e->GetKW( scaleIx);
      if( scaleKW != NULL)
	{
	  DFloat scaleValue;
	  e->AssureFloatScalarKW( scaleIx, scaleValue);
	  bool success = actDevice->SetScale( scaleValue);
	  if( !success)
	    e->Throw( "Current device does not support keyword SCALE.");
	} 
    }

    // COLOR 
    {
      // TODO: turn off with COLOR=0?
      static int colorIx = e->KeywordIx( "COLOR");
      BaseGDL* colorKW = e->GetKW( colorIx);
      if( colorKW != NULL)
	{
	  DLong colorValue;
	  e->AssureLongScalarKW( colorIx, colorValue);
	  bool success = actDevice->SetColor(colorValue);
	  if( !success) e->Throw( "Current device does not support keyword COLOR.");
	} 
    }

    // ENCAPSULATED 
    {
      static int encapsulatedIx = e->KeywordIx( "ENCAPSULATED");
      BaseGDL* encapsulatedKW = e->GetKW( encapsulatedIx);
      if( encapsulatedKW != NULL)
	{
	  bool success;
	  if ((*e->GetKWAs<DIntGDL>(encapsulatedIx))[0] == 0)
	    success = actDevice->SetEncapsulated(false);
	  else
	    success = actDevice->SetEncapsulated(true);
	  if (!success) e->Throw( "Current device does not support keyword ENCAPSULATED.");
	} 
    }
    //BITS_PER_PIXEL
    {
      static int bppIx = e->KeywordIx( "BITS_PER_PIXEL");
      BaseGDL* bppKW = e->GetKW( bppIx);
      if( bppKW != NULL)
	{
	  bool success;
      
	  success = actDevice->SetBPP((*e->GetKWAs<DIntGDL>(bppIx))[0]);
	  if (!success) e->Throw( "Current device does not support keyword BITS_PER_PIXEL.");
	} 
    }    
    // COPY
    {
      static int copyIx = e->KeywordIx("COPY");
      if( e->KeywordPresent( copyIx)) 
      {
       DLongGDL* parameters=e->GetKWAs<DLongGDL>(copyIx);
       if (parameters->N_Elements() > 7 || parameters->N_Elements() < 6) 
         e->Throw("Keyword array parameter COPY must have from 6 to 7 elements.");
// WRONG!! device number may also be a gui number in some cases... HAVING TWO X11-like GRAPHIC DEVICES IS NOT GOOD AT ALL! FIXME!
       if (parameters->N_Elements() == 7) {
         if (!actDevice->WState((*parameters)[6])) e->Throw("Window number "+i2s((*parameters)[6])+" out of range or no more windows or known bug with widgets");
       }
       bool doesTheCopy=actDevice->CopyRegion(parameters);
       if (!doesTheCopy) e->Throw( "Keyword COPY not allowed for call to: DEVICE");
       }
    }
 
  }
Ejemplo n.º 3
0
  void hdf_sd_adddata_pro( EnvT* e)
  {
    SizeT nParam=e->NParam();
    char fieldname[256];
    int32 rank;
    int32 dims[MAXRANK];
    int32 start[MAXRANK];
    int32 edges[MAXRANK];
    int32 stride[MAXRANK];
    int32 dtype;
    int32 nattrs;
    int32 status;

    DLong sds_id;
    e->AssureScalarPar<DLongGDL>( 0, sds_id);

    DLongGDL* startKW = e->IfDefGetKWAs<DLongGDL>( 0);
    DLongGDL* strideKW = e->IfDefGetKWAs<DLongGDL>( 1);
    DLongGDL* countKW = e->IfDefGetKWAs<DLongGDL>( 2);

    status = SDgetinfo(sds_id, fieldname, &rank, dims, &dtype, &nattrs);

    if (status == FAIL)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_ADDDATA: Invalid SD dataset ID: " + i2s(sds_id));

    for( SizeT i = 0; i < rank; i++) {
      start[i] = 0;
      stride[i] = 1;
      edges[i] = dims[i];
    }

    if (startKW != NULL) {
      if (startKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of START dimensions:"
			    +i2s(startKW->N_Elements()));
      memcpy(start, &(*startKW)[0], rank*sizeof(int32));
    }

    if (countKW != NULL) {
      if (countKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of COUNT dimensions:"
			    +i2s(countKW->N_Elements()));
      memcpy(edges, &(*countKW)[0], rank*sizeof(int32));
    }

    if (strideKW != NULL) {
      if (strideKW->N_Elements() != rank)
	throw GDLException( e->CallingNode(), 
			    "HDF_SD_GETDATA: Incorrect number of STRIDE dimensions:"
			    +i2s(strideKW->N_Elements()));
      memcpy(stride, &(*strideKW)[0], rank*sizeof(int32));
    }


    BaseGDL* p1 = e->GetParDefined( 1);
    DByteGDL* p1B = static_cast<DByteGDL*>( p1);

    if (strideKW == NULL)
      SDwritedata(sds_id, start, NULL, edges, (VOIDP) &(*p1B)[0]);
    else
      SDwritedata(sds_id, start, stride, edges, (VOIDP) &(*p1B)[0]);
  }