Beispiel #1
0
        public void generateEFieldLines()
        {
            if (charges.Count == 0)
            {
                eFieldLineVerts = null;
                return;
            }

            int numLinesPerCharge = 25;

            List <VertexPositionColor> tempVerts = new List <VertexPositionColor>();

            Random random = new Random();
            int    lineSegmentsPerLine = 65535 / (charges.Count * numLinesPerCharge);
            float  minLon = 0;
            float  maxLon = (float)Math.PI;

            float minLat = (float)-Math.PI;
            float maxLat = (float)Math.PI;

            int numLoops = (int)Math.Sqrt(numLinesPerCharge);

            for (int i = 0; i < charges.Count; i++)
            {
                // Non-random line starting algorithm
                PointCharge charge = charges[i];
                for (int lonIndex = 0; lonIndex < numLoops; lonIndex++)
                {
                    float lon = MathHelper.Lerp(minLon, maxLon, (float)lonIndex / (float)numLoops);
                    for (int latIndex = 0; latIndex < numLoops; latIndex++)
                    {
                        float   lat         = MathHelper.Lerp(minLat, maxLat, (float)latIndex / (float)numLoops);
                        Vector3 startOffset = 2 * charge.radius * Vector3.Transform(Vector3.UnitX, Matrix.CreateRotationY(lon) * Matrix.CreateRotationX(lat));
                        Vector3 startPos    = charge.center + startOffset;

                        VertexPositionColor[] line = calculateVerticesForEFieldLineFromStartPosition(startPos, charge.charge < 0, lineSegmentsPerLine);
                        tempVerts.AddRange(line);
                    }
                }


                // Random line starting algorithm

                /*PointCharge charge = charges[i];
                 * for (int n = 0; n < numLinesPerCharge; n++)
                 * {
                 *  float lon = (float)(random.NextDouble() * Math.PI);
                 *  float lat = (float)(random.NextDouble() * MathHelper.TwoPi - Math.PI);
                 *
                 *  Vector3 startOffset = 2 * charge.radius * Vector3.Transform(Vector3.UnitX, Matrix.CreateRotationY(lon) * Matrix.CreateRotationX(lat));
                 *  Vector3 startPos = charge.center + startOffset;
                 *
                 *  VertexPositionColor[] line = calculateVerticesForEFieldLineFromStartPosition(startPos, charge.charge < 0, lineSegmentsPerLine);
                 *  tempVerts.AddRange(line);
                 * }*/
            }

            eFieldLineVerts = tempVerts.ToArray();
        }
Beispiel #2
0
 public PointCharge GetChargeWithID(int id)
 {
     for (int i = 0; i < charges.Count; i++)
     {
         PointCharge c = charges[i];
         if (c.id == id)
         {
             return(c);
         }
     }
     return(null);
 }
Beispiel #3
0
        public void AddCharge(PointCharge charge)
        {
            idCounter++;
            charges.Add(charge);

            float minX = float.MaxValue;
            float minY = float.MaxValue;
            float minZ = float.MaxValue;

            float maxX = float.MinValue;
            float maxY = float.MinValue;
            float maxZ = float.MinValue;

            for (int i = 0; i < charges.Count; i++)
            {
                PointCharge c = charges[i];
                if (c.center.X < minX)
                {
                    minX = c.center.X;
                }
                if (c.center.Y < minY)
                {
                    minY = c.center.Y;
                }
                if (c.center.Z < minZ)
                {
                    minZ = c.center.Z;
                }

                if (c.center.X > maxX)
                {
                    maxX = c.center.X;
                }
                if (c.center.Y > maxY)
                {
                    maxY = c.center.Y;
                }
                if (c.center.Z > maxZ)
                {
                    maxZ = c.center.Z;
                }
            }

            minX -= 5;
            minY -= 5;
            minZ -= 5;

            maxX += 5;
            maxY += 5;
            maxZ += 5;
            this.SetBounds(minX, minY, minZ, maxX, maxY, maxZ);
        }
Beispiel #4
0
        protected override void LoadContent()
        {
            device = this.GraphicsDevice;

            cam = new Camera(new Vector3(-5, 5, 5), -MathHelper.PiOver4, -MathHelper.PiOver4, 0.001f, 0.1f, device);


            // Setup axes grid
            grid = new Grid(device, 0.1f);

            // Test charge
            testCharge = new PointCharge(device, 1.0f, -1.0f);
        }
Beispiel #5
0
 public void DeleteChargeWithID(int id)
 {
     for (int i = 0; i < charges.Count; i++)
     {
         PointCharge c = charges[i];
         if (c.id == id)
         {
             charges.RemoveAt(i);
             generateArrows();
             generateEFieldLines();
             return;
         }
     }
 }
Beispiel #6
0
        public static void ModifyChargeID(int id, string property, float value)
        {
            PointCharge oldC = sim.GetField().GetChargeWithID(id);

            if (oldC == null)
            {
                console.Log("Charge with ID " + id + " not defined");
                return;
            }

            float x      = oldC.center.X;
            float y      = oldC.center.Y;
            float z      = oldC.center.Z;
            float charge = oldC.charge * 1000000;

            if (property == "x")
            {
                x = value;
            }
            else if (property == "y")
            {
                y = value;
            }
            else if (property == "z")
            {
                z = value;
            }
            else if (property == "charge")
            {
                charge = value;
            }
            else
            {
                console.Log("Invalid property \"" + value + "\"");
                return;
            }

            PointCharge c = new PointCharge(sim.GraphicsDevice, new Vector3(x, y, z), charge, id);

            sim.GetField().DeleteChargeWithID(id);
            sim.GetField().AddCharge(c);

            if (shouldAutogen)
            {
                sim.GetField().generateArrows();
                sim.GetField().generateEFieldLines();
            }
        }
Beispiel #7
0
        public override Vector3 r(float x, float y, float z)
        {
            Vector3 pos = new Vector3(x, y, z);

            Vector3 sum = Vector3.Zero;

            for (int i = 0; i < charges.Count; i++)
            {
                PointCharge charge   = charges[i];
                float       distance = Vector3.Distance(pos, charge.center);

                float magnitude = -k * charge.charge / (distance * distance);

                Vector3 vec = charge.center - pos;
                vec.Normalize();
                vec *= magnitude;
                sum += vec;
            }


            return(sum);
        }
Beispiel #8
0
        public void AddCharge(Vector3 pos, float chargeInMicroCoulombs)
        {
            PointCharge pointCharge = new PointCharge(this.device, pos, chargeInMicroCoulombs, idCounter);

            AddCharge(pointCharge);
        }