ハート型のコースティクス

2009-10-25
#graphics

ハート型のコースティクス

smallppmに円柱の当たり判定を組み込んでみた。 ちゃんとハート型のコースティクスができて満足（カージオイドというらしい）、かなり鋭いエッジが出てる。 右の青い反射面からのハートもうっすらと出てる。描画は裏で動かしながらで、１億フォトン/71分。

あとはこれみたいに水面のコースティクスとかやってみたいけど、そうするとポリゴンモデルとの当たり判定を作らなきゃいけなくて、高速化の手法を知らないのでどうしたらいいか悩み中。

  • 水面より上にコースティクスが出てないので変だなーと思ってたら動画を見てわかった。水面は屈折だけで反射してないぽい。
  • こっちはより一層いいな！
  • これも

円柱：

struct Cylinder : public Object {
  double rad, h; Vec p;
  Cylinder(double r_,double h_,Vec p_,Vec c_,Refl_t re_) : Object(c_,re_), rad(r_),h(h_),p(p_){}
  virtual double intersect(const Ray &r, Vec* nrm) const { // returns distance
    Vec o = r.o - p;
    double a = sq(r.d.x) + sq(r.d.z);
    if (a == 0)  return 1e20;
    double b = 2 * (o.x * r.d.x + o.z * r.d.z);
    double c = sq(o.x) + sq(o.z) - sq(rad);
    double D = sq(b) - 4 * a * c;
    if (D <= 0) {
      return 1e20;
    } else {
      double det = sqrt(D);
      double t;
      if ((t = (-b - det) / (2 * a)) > 1e-4 && (fabs(o.y + r.d.y*t) <= h)) {
        // hit
      } else if ((t = (-b + det) / (2 * a)) > 1e-4 && (fabs(o.y + r.d.y*t) <= h)) {
        // hit
      } else {
        return 1e20;
      }

      Vec x=r.o+r.d*t;
      Vec n = (x-this->p); n.y = 0;
      *nrm = n.norm();
      return t;
    }
  }
};

関連記事

新着記事