/*------------------------------------

 Iwa_DirectionalBlurFx

 ボケ足の伸ばし方を選択でき、参照画像を追加した DirectionalBlur

//------------------------------------*/

#include "iwa_directionalblurfx.h"

#include "tparamuiconcept.h"

enum FILTER_TYPE { Linear = 0, Gaussian, Flat };

/*------------------------------------------------------------

 参照画像の輝度を０〜１に正規化してホストメモリに読み込む

------------------------------------------------------------*/

template <typename pixel="" raster,="" typename=""></typename>

void Iwa_DirectionalBlurFx::setReferenceRaster(const RASTER srcRas,

                                               float *dstMem, TDimensionI dim) {

  float *dst_p = dstMem;

  for (int j = 0; j < dim.ly; j++) {

    PIXEL *pix = srcRas->pixels(j);

    for (int i = 0; i < dim.lx; i++, pix++, dst_p++) {

      (*dst_p) = ((float)pix->r * 0.3f + (float)pix->g * 0.59f +

                  (float)pix->b * 0.11f) /

                 (float)PIXEL::maxChannelValue;

    }

  }

}

/*------------------------------------------------------------

 ソース画像を０〜１に正規化してホストメモリに読み込む

------------------------------------------------------------*/

template <typename pixel="" raster,="" typename=""></typename>

void Iwa_DirectionalBlurFx::setSourceRaster(const RASTER srcRas, float4 *dstMem,

                                            TDimensionI dim) {

  float4 *chann_p = dstMem;

  for (int j = 0; j < dim.ly; j++) {

    PIXEL *pix = srcRas->pixels(j);

    for (int i = 0; i < dim.lx; i++, pix++, chann_p++) {

      (*chann_p).x = (float)pix->r / (float)PIXEL::maxChannelValue;

      (*chann_p).y = (float)pix->g / (float)PIXEL::maxChannelValue;

      (*chann_p).z = (float)pix->b / (float)PIXEL::maxChannelValue;

      (*chann_p).w = (float)pix->m / (float)PIXEL::maxChannelValue;

    }

  }

}

/*------------------------------------------------------------

 出力結果をChannel値に変換してタイルに格納

------------------------------------------------------------*/

template <typename pixel="" raster,="" typename=""></typename>

void Iwa_DirectionalBlurFx::setOutputRaster(float4 *srcMem, const RASTER dstRas,

                                            TDimensionI dim, int2 margin) {

  int out_j = 0;

  for (int j = margin.y; j < dstRas->getLy() + margin.y; j++, out_j++) {

    PIXEL *pix     = dstRas->pixels(out_j);

    float4 *chan_p = srcMem;

    chan_p += j * dim.lx + margin.x;

    for (int i = 0; i < dstRas->getLx(); i++, pix++, chan_p++) {

      float val;

      val    = (*chan_p).x * (float)PIXEL::maxChannelValue + 0.5f;

      pix->r = (typename PIXEL::Channel)((val > (float)PIXEL::maxChannelValue)

                                             ? (float)PIXEL::maxChannelValue

                                             : val);

      val    = (*chan_p).y * (float)PIXEL::maxChannelValue + 0.5f;

      pix->g = (typename PIXEL::Channel)((val > (float)PIXEL::maxChannelValue)

                                             ? (float)PIXEL::maxChannelValue

                                             : val);

      val    = (*chan_p).z * (float)PIXEL::maxChannelValue + 0.5f;

      pix->b = (typename PIXEL::Channel)((val > (float)PIXEL::maxChannelValue)

                                             ? (float)PIXEL::maxChannelValue

                                             : val);

      val    = (*chan_p).w * (float)PIXEL::maxChannelValue + 0.5f;

      pix->m = (typename PIXEL::Channel)((val > (float)PIXEL::maxChannelValue)

                                             ? (float)PIXEL::maxChannelValue

                                             : val);

    }

  }

}

//------------------------------------

Iwa_DirectionalBlurFx::Iwa_DirectionalBlurFx()

    : m_angle(0.0)

    , m_intensity(10.0)

    , m_bidirectional(false)

    , m_filterType(new TIntEnumParam(Linear, "Linear")) {

  m_intensity->setMeasureName("fxLength");

  m_angle->setMeasureName("angle");

  bindParam(this, "angle", m_angle);

  bindParam(this, "intensity", m_intensity);

  bindParam(this, "bidirectional", m_bidirectional);

  bindParam(this, "filterType", m_filterType);

  addInputPort("Source", m_input);

  addInputPort("Reference", m_reference);

  m_intensity->setValueRange(0, (std::numeric_limits<double>::max)());</double>

  m_filterType->addItem(Gaussian, "Gaussian");

  m_filterType->addItem(Flat, "Flat");

}

//------------------------------------

void Iwa_DirectionalBlurFx::doCompute(TTile &tile, double frame,

                                      const TRenderSettings &settings) {

  /*- 接続していない場合は処理しない -*/

  if (!m_input.isConnected()) {

    tile.getRaster()->clear();

    return;

  }

  TPointD blurVector;

  double angle       = m_angle->getValue(frame) * M_PI_180;

  double intensity   = m_intensity->getValue(frame);

  bool bidirectional = m_bidirectional->getValue();

  blurVector.x = intensity * cos(angle);

  blurVector.y = intensity * sin(angle);

  double shrink = (settings.m_shrinkX + settings.m_shrinkY) / 2.0;

  /*- 平行移動成分を消したアフィン変換をかけ、ボケベクトルを変形 -*/

  TAffine aff = settings.m_affine;

  aff.a13 = aff.a23 = 0;

  TPointD blur      = (1.0 / shrink) * (aff * blurVector);

  /*- ボケの実際の長さを求める -*/

  double blurValue = norm(blur);

  /*- ボケの実際の長さがほぼ０なら入力画像をそのまま return -*/

  if (areAlmostEqual(blurValue, 0, 1e-3)) {

    m_input->compute(tile, frame, settings);

    return;

  }

  /*-  表示の範囲を得る -*/

  TRectD bBox =

      TRectD(tile.m_pos /*- Render画像上(Pixel単位)の位置  -*/

             ,

             TDimensionD(/*- Render画像上(Pixel単位)のサイズ  -*/

                         tile.getRaster()->getLx(), tile.getRaster()->getLy()));

  /*- 上下左右のマージンを得る -*/

  double minX, maxX, minY, maxY;

  if (blur.x > 0.0) /*- X成分が正の場合 -*/

  {

    maxX = blur.x;

    minX = (bidirectional) ? -blur.x : 0.0;

  } else /*- X成分が負の場合 -*/

  {

    maxX = (bidirectional) ? -blur.x : 0.0;

    minX = blur.x;

  }

  if (blur.y > 0.0) /*- Y成分が正の場合 -*/

  {

    maxY = blur.y;

    minY = (bidirectional) ? -blur.y : 0.0;

  } else /*-  Y成分が負の場合 -*/

  {

    maxY = (bidirectional) ? -blur.y : 0.0;

    minY = blur.y;

  }

  int marginLeft   = (int)ceil(abs(minX));

  int marginRight  = (int)ceil(abs(maxX));

  int marginTop    = (int)ceil(abs(maxY));

  int marginBottom = (int)ceil(abs(minY));

  /*- マージンは、フィルタの上下左右を反転した寸法になる -*/

  TRectD enlargedBBox(bBox.x0 - (double)marginRight,

                      bBox.y0 - (double)marginTop, bBox.x1 + (double)marginLeft,

                      bBox.y1 + (double)marginBottom);

  std::cout << "Margin Left:" << marginLeft << " Right:" << marginRight

            << " Bottom:" << marginBottom << " Top:" << marginTop << std::endl;

  TDimensionI enlargedDimIn(/*- Pixel単位に四捨五入 -*/

                            (int)(enlargedBBox.getLx() + 0.5),

                            (int)(enlargedBBox.getLy() + 0.5));

  TTile enlarge_tile;

  m_input->allocateAndCompute(enlarge_tile, enlargedBBox.getP00(),

                              enlargedDimIn, tile.getRaster(), frame, settings);

  /*- 参照画像が有ったら、メモリに取り込む -*/

  float *reference_host = 0;

  TRasterGR8P reference_host_ras;

  if (m_reference.isConnected()) {

    TTile reference_tile;

    m_reference->allocateAndCompute(reference_tile, enlargedBBox.getP00(),

                                    enlargedDimIn, tile.getRaster(), frame,

                                    settings);

    /*- ホストのメモリ確保 -*/

    reference_host_ras =

        TRasterGR8P(sizeof(float) * enlargedDimIn.lx, enlargedDimIn.ly);

    reference_host_ras->lock();

    reference_host = (float *)reference_host_ras->getRawData();

    /*- 参照画像の輝度を０〜１に正規化してホストメモリに読み込む -*/

    TRaster32P ras32 = (TRaster32P)reference_tile.getRaster();

    TRaster64P ras64 = (TRaster64P)reference_tile.getRaster();

    if (ras32)

      setReferenceRaster<traster32p, tpixel32="">(ras32, reference_host,</traster32p,>

                                               enlargedDimIn);

    else if (ras64)

      setReferenceRaster<traster64p, tpixel64="">(ras64, reference_host,</traster64p,>

                                               enlargedDimIn);

  }

  //-------------------------------------------------------

  /*- 計算範囲 -*/

  TDimensionI dimOut(tile.getRaster()->getLx(), tile.getRaster()->getLy());

  TDimensionI filterDim(marginLeft + marginRight + 1,

                        marginTop + marginBottom + 1);

  doCompute_CPU(tile, frame, settings, blur, bidirectional, marginLeft,

                marginRight, marginTop, marginBottom, enlargedDimIn,

                enlarge_tile, dimOut, filterDim, reference_host);

  /*-  参照画像が刺さっている場合、メモリを解放する -*/

  if (reference_host) {

    reference_host_ras->unlock();

    reference_host = 0;

  }

}

//------------------------------------

void Iwa_DirectionalBlurFx::doCompute_CPU(

    TTile &tile, double frame, const TRenderSettings &settings, TPointD &blur,

    bool bidirectional, int marginLeft, int marginRight, int marginTop,

    int marginBottom, TDimensionI &enlargedDimIn, TTile &enlarge_tile,

    TDimensionI &dimOut, TDimensionI &filterDim, float *reference_host) {

  /*- メモリ確保 -*/

  TRasterGR8P in_ras(sizeof(float4) * enlargedDimIn.lx, enlargedDimIn.ly);

  in_ras->lock();

  float4 *in = (float4 *)in_ras->getRawData();

  TRasterGR8P out_ras(sizeof(float4) * enlargedDimIn.lx, enlargedDimIn.ly);

  out_ras->lock();

  float4 *out = (float4 *)out_ras->getRawData();

  TRasterGR8P filter_ras(sizeof(float) * filterDim.lx, filterDim.ly);

  filter_ras->lock();

  float *filter = (float *)filter_ras->getRawData();

  /*- ソース画像を０〜１に正規化してホストメモリに読み込む -*/

  TRaster32P ras32 = (TRaster32P)enlarge_tile.getRaster();

  TRaster64P ras64 = (TRaster64P)enlarge_tile.getRaster();

  if (ras32)

    setSourceRaster<traster32p, tpixel32="">(ras32, in, enlargedDimIn);</traster32p,>

  else if (ras64)

    setSourceRaster<traster64p, tpixel64="">(ras64, in, enlargedDimIn);</traster64p,>

  /*- フィルタ作る -*/

  makeDirectionalBlurFilter_CPU(filter, blur, bidirectional, marginLeft,

                                marginRight, marginTop, marginBottom,

                                filterDim);

  if (reference_host) /*- 参照画像がある場合 -*/

  {

    float4 *out_p = out + marginTop * enlargedDimIn.lx;

    float *ref_p  = reference_host + marginTop * enlargedDimIn.lx;

    /*- フィルタリング -*/

    for (int y = marginTop; y < dimOut.ly + marginTop; y++) {

      out_p += marginRight;

      ref_p += marginRight;

      for (int x = marginRight; x < dimOut.lx + marginRight;

           x++, out_p++, ref_p++) {

        /*- 参照画像が黒ならソースをそのまま返す -*/

        if ((*ref_p) == 0.0f) {

          (*out_p) = in[y * enlargedDimIn.lx + x];

          continue;

        }

        /*- 値を積算する入れ物を用意 -*/

        float4 value = {0.0f, 0.0f, 0.0f, 0.0f};

        float *filter_p = filter;

        if ((*ref_p) == 1.0f) {

          /*- フィルタのサイズでループ

             ただし、フィルタはサンプル点の画像を収集するように

                  用いるため、上下左右反転してサンプルする -*/

          for (int fily = -marginBottom; fily < filterDim.ly - marginBottom;

               fily++) {

            /*- サンプル座標 と インデックス の このスキャンラインのうしろ -*/

            int2 samplePos  = {x + marginLeft, y - fily};

            int sampleIndex = samplePos.y * enlargedDimIn.lx + samplePos.x;

            for (int filx = -marginLeft; filx < filterDim.lx - marginLeft;

                 filx++, filter_p++, sampleIndex--) {

              /*- フィルター値が０またはサンプルピクセルが透明ならcontinue -*/

              if ((*filter_p) == 0.0f || in[sampleIndex].w == 0.0f) continue;

              /*- サンプル点の値にフィルタ値を掛けて積算する -*/

              value.x += in[sampleIndex].x * (*filter_p);

              value.y += in[sampleIndex].y * (*filter_p);

              value.z += in[sampleIndex].z * (*filter_p);

              value.w += in[sampleIndex].w * (*filter_p);

            }

          }

        } else {

          for (int fily = -marginBottom; fily < filterDim.ly - marginBottom;

               fily++) {

            for (int filx = -marginLeft; filx < filterDim.lx - marginLeft;

                 filx++, filter_p++) {

              /*- フィルター値が０ならcontinue -*/

              if ((*filter_p) == 0.0f) continue;

              /*- サンプル座標 -*/

              int2 samplePos = {tround((float)x - (float)filx * (*ref_p)),

                                tround((float)y - (float)fily * (*ref_p))};

              int sampleIndex = samplePos.y * enlargedDimIn.lx + samplePos.x;

              /*- サンプルピクセルが透明ならcontinue -*/

              if (in[sampleIndex].w == 0.0f) continue;

              /*- サンプル点の値にフィルタ値を掛けて積算する -*/

              value.x += in[sampleIndex].x * (*filter_p);

              value.y += in[sampleIndex].y * (*filter_p);

              value.z += in[sampleIndex].z * (*filter_p);

              value.w += in[sampleIndex].w * (*filter_p);

            }

          }

        }

        /*- 値を格納 -*/

        (*out_p) = value;

      }

      out_p += marginLeft;

      ref_p += marginLeft;

    }

  } else /*- 参照画像が無い場合 -*/

  {

    float4 *out_p = out + marginTop * enlargedDimIn.lx;

    /*- フィルタリング -*/

    for (int y = marginTop; y < dimOut.ly + marginTop; y++) {

      out_p += marginRight;

      for (int x = marginRight; x < dimOut.lx + marginRight; x++, out_p++) {

        /*- 値を積算する入れ物を用意 -*/

        float4 value = {0.0f, 0.0f, 0.0f, 0.0f};

        /*- フィルタのサイズでループ

           ただし、フィルタはサンプル点の画像を収集するように

                用いるため、上下左右反転してサンプルする -*/

        int filterIndex = 0;

        for (int fily = -marginBottom; fily < filterDim.ly - marginBottom;

             fily++) {

          /*- サンプル座標 と インデックス の このスキャンラインのうしろ -*/

          int2 samplePos  = {x + marginLeft, y - fily};

          int sampleIndex = samplePos.y * enlargedDimIn.lx + samplePos.x;

          for (int filx = -marginLeft; filx < filterDim.lx - marginLeft;

               filx++, filterIndex++, sampleIndex--) {

            /*- フィルター値が０またはサンプルピクセルが透明ならcontinue -*/

            if (filter[filterIndex] == 0.0f || in[sampleIndex].w == 0.0f)

              continue;

            /*- サンプル点の値にフィルタ値を掛けて積算する -*/

            value.x += in[sampleIndex].x * filter[filterIndex];

            value.y += in[sampleIndex].y * filter[filterIndex];

            value.z += in[sampleIndex].z * filter[filterIndex];

            value.w += in[sampleIndex].w * filter[filterIndex];

          }

        }

        /*- 値を格納 -*/

        (*out_p) = value;

      }

      out_p += marginLeft;

    }

  }

  in_ras->unlock();

  filter_ras->unlock();

  /*- ラスタのクリア -*/

  tile.getRaster()->clear();

  TRaster32P outRas32 = (TRaster32P)tile.getRaster();

  TRaster64P outRas64 = (TRaster64P)tile.getRaster();

  int2 margin         = {marginRight, marginTop};

  if (outRas32)

    setOutputRaster<traster32p, tpixel32="">(out, outRas32, enlargedDimIn, margin);</traster32p,>

  else if (outRas64)

    setOutputRaster<traster64p, tpixel64="">(out, outRas64, enlargedDimIn, margin);</traster64p,>

  out_ras->unlock();

}

//------------------------------------

void Iwa_DirectionalBlurFx::makeDirectionalBlurFilter_CPU(

    float *filter, TPointD &blur, bool bidirectional, int marginLeft,

    int marginRight, int marginTop, int marginBottom, TDimensionI &filterDim) {

  /*- 必要なら、ガウスフィルタを前計算 -*/

  FILTER_TYPE filterType = (FILTER_TYPE)m_filterType->getValue();

  std::vector<float> gaussian;</float>

  if (filterType == Gaussian) {

    gaussian.reserve(101);

    for (int g = 0; g < 101; g++) {

      float x = (float)g / 100.0f;

      // sigma == 0.25

      gaussian.push_back(exp(-x * x / 0.125f));

    }

  }

  /*- フィルタを作る -*/

  TPointD p0 =

      (bidirectional) ? TPointD(-blur.x, -blur.y) : TPointD(0.0f, 0.0f);

  TPointD p1       = blur;

  float2 vec_p0_p1 = {static_cast<float>(p1.x - p0.x),</float>

                      static_cast<float>(p1.y - p0.y)};</float>

  float *fil_p        = filter;

  float intensity_sum = 0.0f;

  for (int fy = -marginBottom; fy <= marginTop; fy++) {

    for (int fx = -marginLeft; fx <= marginRight; fx++, fil_p++) {

      /*- 現在の座標とブラー直線の距離を求める -*/

      /*- P0->サンプル点とP0->P1の内積を求める -*/

      float2 vec_p0_sample = {static_cast<float>(fx - p0.x),</float>

                              static_cast<float>(fy - p0.y)};</float>

      float dot = vec_p0_sample.x * vec_p0_p1.x + vec_p0_sample.y * vec_p0_p1.y;

      /*- 軌跡ベクトルの長さの２乗を計算する -*/

      float length2 = vec_p0_p1.x * vec_p0_p1.x + vec_p0_p1.y * vec_p0_p1.y;

      /*- 距離の2乗を求める -*/

      float dist2;

      float framePosRatio;

      /*- P0より手前にある場合 -*/

      if (dot <= 0.0f) {

        dist2 = vec_p0_sample.x * vec_p0_sample.x +

                vec_p0_sample.y * vec_p0_sample.y;

        framePosRatio = 0.0f;

      } else {

        /*- P0〜P1間にある場合 -*/

        if (dot < length2) {

          float p0_sample_dist2 = vec_p0_sample.x * vec_p0_sample.x +

                                  vec_p0_sample.y * vec_p0_sample.y;

          dist2         = p0_sample_dist2 - dot * dot / length2;

          framePosRatio = dot / length2;

        }

        /*- P1より先にある場合 -*/

        else {

          float2 vec_p1_sample = {static_cast<float>(fx - p1.x),</float>

                                  static_cast<float>(fy - p1.y)};</float>

          dist2 = vec_p1_sample.x * vec_p1_sample.x +

                  vec_p1_sample.y * vec_p1_sample.y;

          framePosRatio = 1.0f;

        }

      }

      /*- 距離が(√2 + 1)/2より遠かったらreturn dist2との比較だから2乗している

       * -*/

      if (dist2 > 1.4571f) {

        (*fil_p) = 0.0f;

        continue;

      }

      /*-

         現在のピクセルのサブピクセル(16*16)が、近傍ベクトルからの距離0.5の範囲にどれだけ

               含まれているかをカウントする -*/

      int count = 0;

      for (int yy = 0; yy < 16; yy++) {

        /*- サブピクセルのY座標 -*/

        float subPosY = (float)fy + ((float)yy - 7.5f) / 16.0f;

        for (int xx = 0; xx < 16; xx++) {

          /*- サブピクセルのX座標 -*/

          float subPosX = (float)fx + ((float)xx - 7.5f) / 16.0f;

          float2 vec_p0_sub = {static_cast<float>(subPosX - p0.x),</float>

                               static_cast<float>(subPosY - p0.y)};</float>

          float sub_dot =

              vec_p0_sub.x * vec_p0_p1.x + vec_p0_sub.y * vec_p0_p1.y;

          /*- 距離の2乗を求める -*/

          float dist2;

          /*- P0より手前にある場合 -*/

          if (sub_dot <= 0.0f)

            dist2 = vec_p0_sub.x * vec_p0_sub.x + vec_p0_sub.y * vec_p0_sub.y;

          else {

            /*- P0〜P1間にある場合 -*/

            if (sub_dot < length2) {

              float p0_sub_dist2 =

                  vec_p0_sub.x * vec_p0_sub.x + vec_p0_sub.y * vec_p0_sub.y;

              dist2 = p0_sub_dist2 - sub_dot * sub_dot / length2;

            }

            /*-  P1より先にある場合 -*/

            else {

              float2 vec_p1_sub = {static_cast<float>(subPosX - p1.x),</float>

                                   static_cast<float>(subPosY - p1.y)};</float>

              dist2 = vec_p1_sub.x * vec_p1_sub.x + vec_p1_sub.y * vec_p1_sub.y;

            }

          }

          /*- 距離の２乗が0.25より近ければカウントをインクリメント -*/

          if (dist2 <= 0.25f) count++;

        }

      }

      /*- 保険 カウントが０の場合はフィールド値０でreturn -*/

      if (count == 0) {

        (*fil_p) = 0.0f;

        continue;

      }

      /*- countは Max256 -*/

      float countRatio = (float)count / 256.0f;

      /*- オフセット値を求める -*/

      float offset =

          (bidirectional) ? abs(framePosRatio * 2.0 - 1.0) : framePosRatio;

      /*- フィルタごとに分ける -*/

      float bokeAsiVal;

      switch (filterType) {

      case Linear:

        bokeAsiVal = 1.0f - offset;

        break;

      case Gaussian: {

        int index   = (int)floor(offset * 100.0f);

        float ratio = offset * 100.0f - (float)index;

        bokeAsiVal =

            gaussian[index] * (1.0f - ratio) + gaussian[index + 1] * ratio;

      } break;

      case Flat:

        bokeAsiVal = 1.0f;

        break;

      default:

        bokeAsiVal = 1.0f - offset;

        break;

      }

      /*- フィールド値の格納 -*/

      (*fil_p) = bokeAsiVal * countRatio;

      intensity_sum += (*fil_p);

    }

  }

  /*- 正規化 -*/

  fil_p = filter;

  for (int f = 0; f < filterDim.lx * filterDim.ly; f++, fil_p++) {

    if ((*fil_p) == 0.0f) continue;

    (*fil_p) /= intensity_sum;

  }

}

//------------------------------------

bool Iwa_DirectionalBlurFx::doGetBBox(double frame, TRectD &bBox,

                                      const TRenderSettings &info) {

  if (false == this->m_input.isConnected()) {

    bBox = TRectD();

    return false;

  }

  bool ret = m_input->doGetBBox(frame, bBox, info);

  if (bBox == TConsts::infiniteRectD) return ret;

  TPointD blur;

  double angle       = m_angle->getValue(frame) * M_PI_180;

  double intensity   = m_intensity->getValue(frame);

  bool bidirectional = m_bidirectional->getValue();

  blur.x = intensity * cos(angle);

  blur.y = intensity * sin(angle);

  /*- 上下左右のマージンを得る -*/

  double minX, maxX, minY, maxY;

  if (blur.x > 0.0) /*- X成分が正の場合 -*/

  {

    maxX = blur.x;

    minX = (bidirectional) ? -blur.x : 0.0;

  } else /*- X成分が負の場合 -*/

  {

    maxX = (bidirectional) ? -blur.x : 0.0;

    minX = blur.x;

  }

  if (blur.y > 0.0) /*- Y成分が正の場合 -*/

  {

    maxY = blur.y;

    minY = (bidirectional) ? -blur.y : 0.0;

  } else /*-  Y成分が負の場合 -*/

  {

    maxY = (bidirectional) ? -blur.y : 0.0;

    minY = blur.y;

  }

  int marginLeft   = (int)ceil(abs(minX));

  int marginRight  = (int)ceil(abs(maxX));

  int marginTop    = (int)ceil(abs(maxY));

  int marginBottom = (int)ceil(abs(minY));

  TRectD enlargedBBox(

      bBox.x0 - (double)marginLeft, bBox.y0 - (double)marginBottom,

      bBox.x1 + (double)marginRight, bBox.y1 + (double)marginTop);

  bBox = enlargedBBox;

  return ret;

}

//------------------------------------

bool Iwa_DirectionalBlurFx::canHandle(const TRenderSettings &info,

                                      double frame) {

  return isAlmostIsotropic(info.m_affine) || m_intensity->getValue(frame) == 0;

}

//------------------------------------

void Iwa_DirectionalBlurFx::getParamUIs(TParamUIConcept *&concepts,

                                        int &length) {

  concepts = new TParamUIConcept[length = 1];

  concepts[0].m_type  = TParamUIConcept::POLAR;

  concepts[0].m_label = "Angle and Intensity";

  concepts[0].m_params.push_back(m_angle);

  concepts[0].m_params.push_back(m_intensity);

}

//------------------------------------

FX_PLUGIN_IDENTIFIER(Iwa_DirectionalBlurFx, "iwa_DirectionalBlurFx")