equalizescaling()

void datrw::su::Coordinates::equalizescaling

(

)

equalize scaling

Definition at line 411 of file sucomanager.cc.

References DATRW_assert, DATRW_debug, gelev, gx, gy, Mcontrol, datrw::su::options::SpatialSampling::scalco, datrw::su::ScalCoo::scale, datrw::su::ScalCoo::scaletopower(), sdepth, datrw::su::ScalCoo::smallestpower(), sx, and sy.

Referenced by main(), and setvaluesin().

    {
      bool debug=false;
      // adjust horizontal coordinates
      DATRW_debug(debug, "Coordinates::equalizescaling",
                  "scales on entry"
                  << " gelev.scale " << this->gelev.scale
                  << " sdepth.scale " << this->sdepth.scale);
      /* is not required when search for the smallest possible power
      this->sx.adjustscale();
      this->sy.adjustscale();
      this->gx.adjustscale();
      this->gy.adjustscale();
      */
      // search for smallest possible power
      int psx=this->sx.smallestpower(Mcontrol.scalco);
      int psy=this->sy.smallestpower(Mcontrol.scalco);
      int pgx=this->gx.smallestpower(Mcontrol.scalco);
      int pgy=this->gy.smallestpower(Mcontrol.scalco);
      int pnew = psx > psy ? psx : psy;
      pnew = pnew > pgx ? pnew : pgx;
      pnew = pnew > pgy ? pnew : pgy;
      // check against largest possible power
      /* is not necessary
      pnew = pnew < sx.power() ? pnew : sx.power();
      pnew = pnew < sy.power() ? pnew : sy.power();
      pnew = pnew < gx.power() ? pnew : gx.power();
      pnew = pnew < gy.power() ? pnew : gy.power();
      */
      this->sx.scaletopower(pnew);
      this->sy.scaletopower(pnew);
      this->gx.scaletopower(pnew);
      this->gy.scaletopower(pnew);
      DATRW_assert(((this->sx.scale==this->sy.scale) &&
                    (this->sy.scale==this->gx.scale) &&
                    (this->gx.scale==this->gy.scale)),
                   "ERROR: inconsistent coordinate scaling");
      // adjust vertical coordinates
      /* is not required when scaling with the smallest possible power
      this->sdepth.adjustscale();
      this->gelev.adjustscale();
      DATRW_debug(debug, "Coordinates::equalizescaling",
                  "scales after adjust"
                  << " gelev.scale " << this->gelev.scale
                  << " sdepth.scale " << this->sdepth.scale);
      */
      int psdepth=this->sdepth.smallestpower(Mcontrol.scalco);
      int pgelev=this->gelev.smallestpower(Mcontrol.scalco);
      pnew = psdepth > pgelev ? psdepth : pgelev;
      /* is not necessary
      pnew = pnew < sdepth.power() ? pnew : sdepth.power();
      pnew = pnew < gelev.power() ? pnew : gelev.power();
      */
      this->sdepth.scaletopower(pnew);
      this->gelev.scaletopower(pnew);
      DATRW_assert((this->sdepth.scale==this->gelev.scale),
                   "ERROR: inconsistent coordinate scaling");
    } // void Coordinates::equalizescaling()

Here is the call graph for this function:

Here is the caller graph for this function: