Template:CMI Example get grid values C: Difference between revisions

From CSDMS
Newer edit →
Automated import of articles
 
Automated import of articles *** existing text overwritten ***
Line 1: Line 1:
<syntaxhighlight lang=c>
    struct sidl_double__array* vals = NULL;


    {
      struct edu_csdms_models_Sedflux3D__data *this =
        edu_csdms_models_Sedflux3D__get_data (self);
      if (this && this->state)
      {
        int n_dims;
        int * dimen;
        double* data = BMI_Get_double (this->state, long_var_name, &n_dims,
            &dimen);
        if (data)
        { /* Create an array that borrows these values. */
          int i;
          int * lower = (int*)malloc (sizeof (int)*n_dims);
          int * upper = (int*)malloc (sizeof (int)*n_dims);
          int * stride = (int*)malloc (sizeof (int)*n_dims);
          /* BMI_get_double returns an array that is compact and has
          * unit stride */
          for (i=0; i<n_dims; i++)
          {
            lower[i] = 0;
            upper[i] = dimen[i]-1;
            //upper[i] = dimen[n_dims-i-1]-1;
          }
          //for (i=1, stride[0]=1; i<n_dims; i++)
          //  stride[i] = stride[i-1]*dimen[i-1];
          for (i=n_dims-2, stride[n_dims-1]=1; i>=0; i--)
            stride[i] = stride[i+1]*dimen[i+1];
          if (n_dims==2)
          {
            int i, i_0;
            fprintf (stderr, "%s: Getting values
", CMI_COMPONENT_NAME);
            fprintf (stderr, "%s: Shape is %dx%d
", CMI_COMPONENT_NAME,
                dimen[0], dimen[1]);
            fprintf (stderr, "%s: Stride is %dx%d
", CMI_COMPONENT_NAME,
                stride[0], stride[1]);
            fflush (stderr);
          }
          //vals = sidl_double__array_borrow (data, n_dims, lower, upper,
          //    stride);
          {
            const int _n_dims = 1;
            const int _lower[1] = {0};
            const int _upper[1] = {dimen[0]*dimen[1]-1};
            const int _stride[1] = {1};
            vals = sidl_double__array_borrow (data, _n_dims, _lower,
                _upper, _stride);
          }
          free (stride);
          free (upper);
          free (lower);
        }
      }
    }
    return vals;
</syntaxhighlight>

Revision as of 17:12, 1 February 2012

Retrieved from "https://csdms.colorado.edu/csdms_wiki/index.php?title=Template:CMI_Example_get_grid_values_C&oldid=33188"