← Back to branch summary

~ubuntu-branches/ubuntu/precise/grass/precise

« back to all changes in this revision

Viewing changes to vector/v.vol.rst/main.c

  • Committer: Bazaar Package Importer
  • Author(s): Francesco Paolo Lovergine
  • Date: 2011-04-13 17:08:41 UTC
  • mfrom: (8.1.7 sid)
  • Revision ID: james.westby@ubuntu.com-20110413170841-ss1t9bic0d0uq0gz
Tags: 6.4.1-1
* New upstream version.
* Now build-dep on libjpeg-dev and current libreadline6-dev.
* Removed patch swig: obsolete.
* Policy bumped to 3.9.2, without changes.

Show diffs side-by-side

added added

removed removed

Lines of Context:
vector/v.vol.rst/main.c
53
53
/* pargr */
54
54
double ns_res, ew_res, tb_res;
55
55
double dmin, ertre;
56
 
int nsizr, nsizc, nsizl, KMAX2, KMIN, KMAX;
 
56
int nsizr, nsizc, nsizl, KMAX2, KMIN, KMAX, KMAXPOINTS;
57
57
 
58
58
/* datgr */
59
59
double *az, *adx, *ady, *adxx, *adyy, *adxy, *adz, *adxz, *adyz, *adzz;
172
172
    int per, npmin;
173
173
    int ii, i, n_rows, n_cols, n_levs;
174
174
    double x_orig, y_orig, z_orig;
175
 
    char dminchar[10];
 
175
    char dminchar[64];
176
176
    struct octdata *data;
177
177
    struct octfunc *functions;
178
178
    struct octtree *tree;
184
184
    struct
185
185
    {
186
186
        struct Option *input, *colnum, *scol, *wheresql, *rescalex, *fi,
187
 
            *segmax, *dmin1, *npmin, *wmult, *outz, *rsm, *maskmap, *zmult,
 
187
            *segmax, *dmin1, *npmin, *npmax, *wmult, *outz, *rsm, *maskmap, *zmult,
188
188
            *cvdev, *gradient, *aspect1, *aspect2, *ncurv, *gcurv, *mcurv,
189
189
            *cellinp, *cellout, *devi;
190
190
    } parm;
318
318
        _("Minimum number of points for approximation in a segment (>segmax)");
319
319
    parm.npmin->guisection = _("Settings");
320
320
 
 
321
    parm.npmax = G_define_option();
 
322
    parm.npmax->key = "npmax";
 
323
    parm.npmax->type = TYPE_INTEGER;
 
324
    parm.npmax->answer = MAXPOINTS;
 
325
    parm.npmax->required = NO;
 
326
    parm.npmax->description =
 
327
        _("Maximum number of points for approximation in a segment (>npmin)");
 
328
    parm.npmax->guisection = _("Settings");
 
329
 
321
330
    parm.dmin1 = G_define_option();
322
331
    parm.dmin1->key = "dmin";
323
332
    parm.dmin1->type = TYPE_DOUBLE;
450
459
    sscanf(parm.rsm->answer, "%lf", &rsm);
451
460
    sscanf(parm.segmax->answer, "%d", &KMAX);
452
461
    sscanf(parm.npmin->answer, "%d", &npmin);
 
462
    sscanf(parm.npmax->answer, "%d", &KMAXPOINTS);
453
463
    sscanf(parm.wmult->answer, "%lf", &wmult);
454
464
    sscanf(parm.zmult->answer, "%lf", &zmult);
455
465
 
480
490
    tb_res_in = tb_res;
481
491
    z_orig = z_orig * zmult;
482
492
    tb_res = tb_res * zmult;
483
 
    KMAX2 = MAXPOINTS;
484
493
 
485
494
    /*    fprintf (stderr, "DMIN = %f\n", dmin); */
486
495
 
550
559
    ii = INPUT(&In, parm.colnum->answer, parm.scol->answer,
551
560
               parm.wheresql->answer);
552
561
 
 
562
    /* now NPOINT is available */
 
563
    KMAX2 =(NPOINT >= KMAXPOINTS ? KMAXPOINTS : NPOINT);
 
564
 
553
565
    Vect_close(&In);
554
566
 
555
567
    if (devi != NULL || cvdev != NULL) {
801
813
                    sprintf(hist.edhist[1],
802
814
                            "dnorm=%f, dmin=%f, wmult=%f, zmult=%f", dnorm,
803
815
                            atof(parm.dmin1->answer), wmult, zmult);
804
 
                    sprintf(hist.edhist[2], "segmax=%d, npmin=%d, rmsdevi=%f",
805
 
                            KMAX, npmin, sqrt(ertot / KMAX2));
 
816
                    sprintf(hist.edhist[2], "segmax=%d, npmin=%d, npmax=%d, rmsdevi=%f",
 
817
                            KMAX, npmin, KMAXPOINTS, sqrt(ertot / KMAX2));
806
818
                    sprintf(hist.edhist[3], "wmin_data=%f, wmax_data=%f",
807
819
                            wmin, wmax);
808
820
                    /* ? sprintf (hist.edhist[4], "wmin_int=%f, wmax_int=%f", wminac, wmaxac); */