public static void Main(string[] args) { OpticalSystem.Builder systemBuilder = new OpticalSystem.Builder(); Lens.Builder lensBuilder = new Lens.Builder() .position(Vector3Pair.position_000_001) .add_surface(1 / 0.031186861, 14.934638, 4.627804137, new Abbe(Abbe.AbbeFormula.AbbeVd, 1.607170, 59.5002, 0.0)) .add_surface(0, 14.934638, 5.417429465) .add_surface(1 / -0.014065441, 12.766446, 3.728230979, new Abbe(Abbe.AbbeFormula.AbbeVd, 1.575960, 41.2999, 0.0)) .add_surface(1 / 0.034678487, 11.918098, 4.417903733) .add_stop(12.066273, 2.288913925) .add_surface(0, 12.372318, 1.499288597, new Abbe(Abbe.AbbeFormula.AbbeVd, 1.526480, 51.4000, 0.0)) .add_surface(1 / 0.035104369, 14.642815, 7.996205852, new Abbe(Abbe.AbbeFormula.AbbeVd, 1.623770, 56.8998, 0.0)) .add_surface(1 / -0.021187519, 14.642815, 85.243965130); systemBuilder.add(lensBuilder); Image.Builder imagePlaneBuilder = new Image.Builder() .position(new Vector3Pair(new Vector3(0, 0, 125.596), Vector3.vector3_001)) .curve(Flat.flat) .shape(new Rectangle(5.0 * 2)); systemBuilder.add(imagePlaneBuilder); PointSource.Builder ps = new PointSource.Builder(PointSource.SourceInfinityMode.SourceAtInfinity, Vector3.vector3_001) .add_spectral_line(SpectralLine.d) .add_spectral_line(SpectralLine.C) .add_spectral_line(SpectralLine.F); systemBuilder.add(ps); RendererSvg renderer = new RendererSvg(800, 400); OpticalSystem system = systemBuilder.build(); Console.WriteLine(system.ToString()); SystemLayout2D systemLayout2D = new SystemLayout2D(); // draw 2d system layout // system.draw_2d_fit(renderer); // system.draw_2d(renderer); systemLayout2D.layout2d(renderer, system); // RayTraceParameters parameters = new RayTraceParameters(system); // // RayTracer rayTracer = new RayTracer(); // parameters.set_default_distribution ( // new Distribution(Pattern.MeridionalDist, 10, 0.999)); // // TODO set save generated state on point source // System.out.println(parameters.sequenceToString(new StringBuilder()).toString()); // // RayTraceResults result = rayTracer.trace(system, parameters); // RayTraceRenderer.draw_2d(renderer, result, false, null); // Console.WriteLine(renderer.write(new StringBuilder()).ToString()); }
public void testBasics() { OpticalSystem.Builder sys = new OpticalSystem.Builder(); Group.Builder g1 = new Group.Builder().position(Vector3Pair.position_000_001); Group.Builder g21 = new Group.Builder().position(Vector3Pair.position_000_001); Group.Builder g22 = new Group.Builder().position(Vector3Pair.position_000_001); Surface.Builder s211 = new Surface.Builder().position(Vector3Pair.position_000_001) .curve(Flat.flat) .shape(Infinite.infinite); Surface.Builder s221 = new Surface.Builder().position(Vector3Pair.position_000_001) .curve(Flat.flat) .shape(Infinite.infinite); g22.add(s221); g1.add(g21); g1.add(g22); g21.add(s211); sys.add(g1); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(s211e.local_position().isEqual(Vector3.vector3_0, 1e-10)); Assert.IsTrue(system1.getPosition(s211e).isEqual(Vector3.vector3_0, 1e-10)); } g21.localPosition(new Vector3(1, 2, 3)); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(s211e.local_position().isEqual(Vector3.vector3_0, 1e-10)); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(1, 2, 3), 1e-10)); } g1.localPosition(new Vector3(3, 2, 1)); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(4, 4, 4), 1e-10)); } s211.localPosition(new Vector3(7, 7, 7)); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(11, 11, 11), 1e-10)); } { OpticalSystem system1 = sys.build(); system1 = sys.updatePosition(s211, new Vector3(9, 9, 9)); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(9, 9, 9), 1e-10)); Assert.IsTrue(s211e.local_position().isEqual(new Vector3(5, 5, 5), 1e-10)); } g21.rotate(0, 0, -90); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(-1, 9, 9), 1e-10)); } { OpticalSystem system1 = sys.build(); system1 = sys.updatePosition(s211, new Vector3(3, 4, 5)); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(s211e.local_position().isEqual(new Vector3(0, 1, 1), 1e-10)); } g1.rotate(0, 0, 90); { OpticalSystem system1 = sys.build(); Element s211e = system1.getGroup(0).getGroup(0).getSurface(0); Assert.IsTrue(system1.getPosition(s211e).isEqual(new Vector3(5, 2, 5), 1e-10)); } }
static void Main(string[] args) { Args arguments = Args.parseArguments(args); if (arguments.filename == null) { Console.WriteLine("Usage: --specfile inputfile [--scenario num] [--skew] [--output layout|spot] [--dump-system]"); Environment.Exit(1); } OpticalBenchDataImporter.LensSpecifications specs = new OpticalBenchDataImporter.LensSpecifications(); specs.parse_file(arguments.filename); OpticalSystem.Builder systemBuilder = OpticalBenchDataImporter.buildSystem(specs, arguments.scenario); double angleOfView = OpticalBenchDataImporter.getAngleOfViewInRadians(specs, arguments.scenario); Vector3 direction = Vector3.vector3_001; if (arguments.skewRays) { // Construct unit vector at an angle // double z1 = cos (angleOfView); // double y1 = sin (angleOfView); // unit_vector = math::Vector3 (0, y1, z1); Matrix3 r = Matrix3.get_rotation_matrix(0, angleOfView); direction = r.times(direction); } PointSource.Builder ps = new PointSource.Builder(PointSource.SourceInfinityMode.SourceAtInfinity, direction) .add_spectral_line(SpectralLine.d) .add_spectral_line(SpectralLine.C) .add_spectral_line(SpectralLine.F); systemBuilder.add(ps); OpticalSystem system = systemBuilder.build(); if (arguments.dumpSystem) { Console.WriteLine(system); } if (arguments.outputType.Equals("layout")) { RendererSvg renderer = new RendererSvg(800, 400); // draw 2d system layout SystemLayout2D systemLayout2D = new SystemLayout2D(); systemLayout2D.layout2d(renderer, system); RayTraceParameters parameters = new RayTraceParameters(system); RayTracer rayTracer = new RayTracer(); parameters.set_default_distribution( new Distribution(Pattern.MeridionalDist, 10, 0.999)); // TODO set save generated state on point source if (arguments.dumpSystem) { Console.WriteLine(parameters.sequenceToString(new StringBuilder()).ToString()); } RayTraceResults result = rayTracer.trace(system, parameters); RayTraceRenderer.draw_2d(renderer, result, false, null); Console.WriteLine(renderer.write(new StringBuilder()).ToString()); } if (arguments.outputType.Equals("spot")) { RendererSvg renderer = new RendererSvg(300, 300, Rgb.rgb_black); renderer = new RendererSvg(300, 300, Rgb.rgb_black); AnalysisSpot spot = new AnalysisSpot(system); spot.draw_diagram(renderer, true); Console.WriteLine(renderer.write(new StringBuilder()).ToString()); } }