public override Point3d getRandomBirthPlace(Libutility util)
 {
     /*Point3d retpos = emitter.getRandEmitPos();
      * double u;
      * double v;
      * _srf.ClosestPoint(retpos, out u, out v);
      * return _srf.PointAt(u, v);*/
     return(emitter.getRandEmitPos());
 }
 public override Point3d getRandomBirthPlace(Libutility util)
 {
     /*Point3d retpos = emitter.getRandEmitPos();
     double u;
     double v;
     _srf.ClosestPoint(retpos, out u, out v);
     return _srf.PointAt(u, v);*/
     return emitter.getRandEmitPos();
 }
 public override float getOffsetTrailValue(int x, int y, int z, Vector3d orient, float viewangle, float offsetangle, float offsetsteps, Libutility util)
 {
     /*Point3d pt;
     Vector3d[] vec;
     _srf.Evaluate((x + 0.5) * u_interval, (y + 0.5) * v_interval, 1, out pt, out vec);
     vec[0].Unitize();
     Vector3d oricopy = new Vector3d(orient.X, orient.Y, orient.Z);
     oricopy.Unitize();
     double dotproduc = Vector3d.Multiply(vec[0], oricopy);
     if (dotproduc > 0.3)
         return -1;
     Point3d _origin = uv_positions[x, y, z];
     Plane oriplane = new Plane(_origin, orient);
     float _u = offsetsteps * util.sinlut[(int)viewangle] * util.coslut[(int)offsetangle];
     float _v = offsetsteps * util.sinlut[(int)viewangle] * util.sinlut[(int)offsetangle];
     float _w = offsetsteps * util.coslut[(int)viewangle];
     Point3d target = oriplane.PointAt(_u, _v, _w);
     Point3d targetid = getIndexByPosition(target.X, target.Y, target.Z);
     int tx = (int)targetid.X;
     int ty = (int)targetid.Y;
     int tz = getWIndex(target.Z);
     if (tx < 0 || tx >= u || ty < 0 || ty >= v || tz < 0 || tz >= w)
         return -1;
     return trail[tx, ty, tz];*/
     float viewRad = viewangle * 3.1416F / 180;
     Point3d _ori = new Point3d(x, y, 0);
     Vector3d _tarOri = orient;
     _tarOri.Transform(Transform.Rotation(viewRad, Plane.WorldXY.ZAxis, _ori));
     _tarOri.Unitize();
     _tarOri *= offsetsteps;
     Point3d _tar = _ori + _tarOri;
     //Point3d _tar = Point3d.Add(_ori, _tarOri);
     int tx = (int)_tar.X;
     int ty = (int)_tar.Y;
     if (tx < 0 || tx >= u || ty < 0 || ty >= v)
         return -1;
     return trail[tx, ty, 0];
 }
 public override Point3d getNeighbourhoodFreePosByIndex(int x, int y, int z, int radius, Libutility util)
 {
     Point3d retpt = new Point3d(-1, -1, -1);
     if (radius < 1)
     {
         return retpt;
     }
     int start_x = x - radius > 0 ? x - radius : 0;
     int start_y = y - radius > 0 ? y - radius : 0;
     int start_z = z - radius > 0 ? z - radius : 0;
     int end_x = x + radius < u ? x + radius : u - 1;
     int end_y = y + radius < v ? y + radius : v - 1;
     int end_z = z + radius < w ? z + radius : w - 1;
     int count = (end_x - start_x + 1) * (end_y - start_y + 1) * (end_z - start_z + 1); //
     //int[] ids = new int[count];
     List<Point3d> freePos = new List<Point3d>();
     for (int i = end_x; i >= start_x; i--)
     {
         for (int j = end_y; j >= start_y; j--)
         {
             //for (int k = end_z; k >= start_z; k--)
             //{
                 if (particle_ids[i, j, 0] == -1)
                 {
                     //retpt = new Point3d(i, j, k);
                     freePos.Add(uv_positions[i, j, 0]);
                     //return retpt;
                 }
             //}
         }
     }
     int lens = freePos.Count;
     if (lens == 0)
         return retpt;
     retpt = freePos[util.getRand(lens)];
     return retpt;
 }
Exemple #5
0
 //public abstract void setBirthPlace(List<Point3d> origin);
 public abstract Point3d getRandomBirthPlace(Libutility util);
Exemple #6
0
 public abstract Point3d getNeighbourhoodFreePosByIndex(int u, int v, int w, int radius, Libutility util);//by index
Exemple #7
0
 public abstract float getOffsetTrailValue(int u, int v, int w, Vector3d orient, float viewangle, float offsetangle, float offsetsteps, Libutility util);
        public override Point3d getNeighbourhoodFreePosByIndex(int x, int y, int z, int radius, Libutility util)
        {
            Point3d retpt = new Point3d(-1, -1, -1);

            if (radius < 1)
            {
                return(retpt);
            }
            int start_x = x - radius > 0 ? x - radius : 0;
            int start_y = y - radius > 0 ? y - radius : 0;
            int start_z = z - radius > 0 ? z - radius : 0;
            int end_x   = x + radius < u ? x + radius : u - 1;
            int end_y   = y + radius < v ? y + radius : v - 1;
            int end_z   = z + radius < w ? z + radius : w - 1;
            int count   = (end_x - start_x + 1) * (end_y - start_y + 1) * (end_z - start_z + 1); //
            //int[] ids = new int[count];
            List <Point3d> freePos = new List <Point3d>();

            for (int i = end_x; i >= start_x; i--)
            {
                for (int j = end_y; j >= start_y; j--)
                {
                    //for (int k = end_z; k >= start_z; k--)
                    //{
                    if (particle_ids[i, j, 0] == -1)
                    {
                        //retpt = new Point3d(i, j, k);
                        freePos.Add(uv_positions[i, j, 0]);
                        //return retpt;
                    }
                    //}
                }
            }
            int lens = freePos.Count;

            if (lens == 0)
            {
                return(retpt);
            }
            retpt = freePos[util.getRand(lens)];
            return(retpt);
        }
        public override float getOffsetTrailValue(int x, int y, int z, Vector3d orient, float viewangle, float offsetangle, float offsetsteps, Libutility util)
        {
            /*Point3d pt;
             * Vector3d[] vec;
             * _srf.Evaluate((x + 0.5) * u_interval, (y + 0.5) * v_interval, 1, out pt, out vec);
             * vec[0].Unitize();
             * Vector3d oricopy = new Vector3d(orient.X, orient.Y, orient.Z);
             * oricopy.Unitize();
             * double dotproduc = Vector3d.Multiply(vec[0], oricopy);
             * if (dotproduc > 0.3)
             *  return -1;
             * Point3d _origin = uv_positions[x, y, z];
             * Plane oriplane = new Plane(_origin, orient);
             * float _u = offsetsteps * util.sinlut[(int)viewangle] * util.coslut[(int)offsetangle];
             * float _v = offsetsteps * util.sinlut[(int)viewangle] * util.sinlut[(int)offsetangle];
             * float _w = offsetsteps * util.coslut[(int)viewangle];
             * Point3d target = oriplane.PointAt(_u, _v, _w);
             * Point3d targetid = getIndexByPosition(target.X, target.Y, target.Z);
             * int tx = (int)targetid.X;
             * int ty = (int)targetid.Y;
             * int tz = getWIndex(target.Z);
             * if (tx < 0 || tx >= u || ty < 0 || ty >= v || tz < 0 || tz >= w)
             *  return -1;
             * return trail[tx, ty, tz];*/
            float    viewRad = viewangle * 3.1416F / 180;
            Point3d  _ori    = new Point3d(x, y, 0);
            Vector3d _tarOri = orient;

            _tarOri.Transform(Transform.Rotation(viewRad, Plane.WorldXY.ZAxis, _ori));
            _tarOri.Unitize();
            _tarOri *= offsetsteps;
            Point3d _tar = _ori + _tarOri;
            //Point3d _tar = Point3d.Add(_ori, _tarOri);
            int tx = (int)_tar.X;
            int ty = (int)_tar.Y;

            if (tx < 0 || tx >= u || ty < 0 || ty >= v)
            {
                return(-1);
            }
            return(trail[tx, ty, 0]);
        }
Exemple #10
0
 /*public override void setBirthPlace(List<Point3d> origin)
  * {
  *  _origins = origin;
  *  //foreach (Point3d pt in origin)
  *  //{
  *  //  setGridCellValueByIndex((int) pt.X, (int) pt.Y, (int) pt.Z, 1, 1);
  *  //}
  * }*/
 public override Point3d getRandomBirthPlace(Libutility util)
 {
     //int length = _origins.Count;
     return(emitter.getRandEmitPos());
     //return _origins[util.getRand(length)];
 }
Exemple #11
0
        /*
         * public trailInfo getOffsetTrailValue(int x, int y, int z, Vector3d orient, float viewangle, float offsetangle, float offsetsteps, Libutility util)
         * {
         * Point3f origin = new Point3f(x, y, z);
         * Plane oriplane = new Plane(origin, orient);
         * float u = offsetsteps * util.sinlut[(int) viewangle] * util.coslut[(int) offsetangle];
         * float v = offsetsteps * util.sinlut[(int) viewangle] * util.sinlut[(int) offsetangle];
         * float w = offsetsteps * util.coslut[(int) viewangle];
         * Point3d target = oriplane.PointAt(u, v, w);
         * int tx = (int) Math.Round(target.X);
         * int ty = (int) Math.Round(target.Y);
         * int tz = (int) Math.Round(target.Z);
         * float tr = 0;
         * if (tx < 0 || tx >= _x || ty < 0 || ty >= _y || tz < 0 || tz >= _z)
         *  tr = -1;
         * else
         *  tr = trail[tx, ty, tz];
         * trailInfo newinfo = new trailInfo(new Point3d(tx, ty, tz), tr);
         * return newinfo;
         * }*/

        public override Point3d getNeighbourhoodFreePosByIndex(int x, int y, int z, int radius, Libutility util)
        {
            Point3d retpt = new Point3d(-1, -1, -1);

            if (radius < 1)
            {
                return(retpt);
            }

            /*
             * int times = 0;
             * int tpz;
             * int tpy;
             * int tpx;
             * do
             * {
             * tpx = util.getRand(x - radius, x + radius + 1);
             * tpy = util.getRand(y - radius, y + radius + 1);
             * tpz = util.getRand(z - radius, z + radius + 1);
             * times++;
             * }
             * while (particle_ids[tpx, tpy, tpz] != -1 && times < 20);
             * if (particle_ids[tpx, tpy, tpz] != -1)
             * return retpt;
             * else
             * return new Point3d(tpx, tpy, tpz);
             */

            int start_x = x - radius > 0 ? x - radius : 0;
            int start_y = y - radius > 0 ? y - radius : 0;
            int start_z = z - radius > 0 ? z - radius : 0;
            int end_x   = x + radius < u ? x + radius : u - 1;
            int end_y   = y + radius < v ? y + radius : v - 1;
            int end_z   = z + radius < w ? z + radius : w - 1;
            int count   = (end_x - start_x + 1) * (end_y - start_y + 1) * (end_z - start_z + 1); //
            //int[] ids = new int[count];
            List <Point3d> freePos = new List <Point3d>();

            for (int i = end_x; i >= start_x; i--)
            {
                for (int j = end_y; j >= start_y; j--)
                {
                    for (int k = end_z; k >= start_z; k--)
                    {
                        if (particle_ids[i, j, k] == -1)
                        {
                            //retpt = new Point3d(i, j, k);
                            freePos.Add(positions[i, j, k]);
                            //return retpt;
                        }
                    }
                }
            }
            int lens = freePos.Count;

            if (lens == 0)
            {
                return(retpt);
            }
            retpt = freePos[util.getRand(lens)];
            return(retpt);
        }
Exemple #12
0
        public override float getOffsetTrailValue(int x, int y, int z, Vector3d orient, float viewangle, float offsetangle, float offsetsteps, Libutility util)
        {
            //Point3f origin = new Point3f(x, y, z);
            Point3d _origin  = positions[x, y, z];
            Plane   oriplane = new Plane(_origin, orient);
            float   _u       = offsetsteps * util.sinlut[(int)viewangle] * util.coslut[(int)offsetangle];
            float   _v       = offsetsteps * util.sinlut[(int)viewangle] * util.sinlut[(int)offsetangle];
            float   _w       = offsetsteps * util.coslut[(int)viewangle];
            //float u = (double) offsetsteps * Math.Sin((double) viewangle * 3.1416F / 180) * Math.Cos((double) offsetangle * 3.1416F / 180);
            //float v = (double)offsetsteps * Math.Sin((double) viewangle * 3.1416F / 180) * Math.Sin((double) offsetangle * 3.1416F / 180);
            //float w = (double)offsetsteps * Math.Cos((double) viewangle * 3.1416F / 180);
            Point3d target = oriplane.PointAt(_u, _v, _w);
            //int tx = (int)Math.Round(target.X);
            //int ty = (int)Math.Round(target.Y);
            //int tz = (int)Math.Round(target.Z);
            int tx = getUIndex(target.X);
            int ty = getVIndex(target.Y);
            int tz = getWIndex(target.Z);

            if (tx < 0 || tx >= u || ty < 0 || ty >= v || tz < 0 || tz >= w)
            {
                return(-1);
            }
            return(trail[tx, ty, tz]);
        }