2026.1.2.8 List<DetectorResultObject> detectorResults = new List<DetectorResultObject>();␍ ␍ ␍ if (WavefrontData == null) {␍ throw new ArgumentException("No wavefront data given.");␍ }␍ if (WavefrontData.Count != 1) {␍ throw new ArgumentException("Only single input field supported.");␍ }␍ ␍ DataArrayGridLess daGridlessWavefrontData = WavefrontData[0] as DataArrayGridLess;␍ ␍ if (daGridlessWavefrontData == null) {␍ throw new ArgumentException("Only gridless input data supported.");␍ }␍ ␍ //calculate radial wavefront information␍ DataArray1D daRadialWavefrontInformation = calculateRadialWavefrontData(daGridlessWavefrontData);␍ ␍ // sample detector output for documents␍ detectorResults.Add(VL_Detectors.CreateDetectorResult(daRadialWavefrontInformation, "Radial Wavefront Information"));␍ ␍ return detectorResults;␍ ␍ /// <summary>␍ /// Defines how the view settings shall be changed after the resulting document windows ␍ /// are shown.␍ /// </summary>␍ /// <param name="formDocument">The document window for which the view settings are adapted.</param>␍ public static void ChangeViewSettingsAfterShow(FormDocument formDocument) {␍ }␍ ␍ ␍ //support method to initialize the cache of the snippet␍ private DataArray1D calculateRadialWavefrontData(DataArrayGridLess inputWavefrontData) {␍ ␍ //Globals.DataDisplay.LogMessage("init");␍ //generate list of positions and phase values for interpolation␍ List<VectorD> listPositionsForInterpolation = new List<VectorD>();␍ List<double> listPhaseValuesForInterpolation = new List<double>();␍ ␍ ␍ double minX = double.MaxValue;␍ double minY = double.MaxValue;␍ double maxX = double.MinValue;␍ double maxY = double.MinValue;␍ ␍ ␍ //get current value␍ double currValue_0 = inputWavefrontData.Data[0][0].Re;␍ ␍ //loop over all point of the input wavefront data to compensate␍ for (int runPoints = 0; runPoints < inputWavefrontData.NoOfDataPoints; runPoints++) {␍ //get position␍ VectorD currPosition = new VectorD(inputWavefrontData.Coordinates[0, runPoints].Re,␍ inputWavefrontData.Coordinates[1, runPoints].Re);␍ //get current value␍ double currValue = inputWavefrontData.Data[0][runPoints].Re;␍ ␍ if (!double.IsNaN(currValue) &&␍ !currPosition.IsUndefined) {␍ listPositionsForInterpolation.Add(currPosition);␍ listPhaseValuesForInterpolation.Add(currValue - currValue_0);␍ ␍ minX = Math.Min(minX, currPosition.X);␍ minY = Math.Min(minY, currPosition.Y);␍ maxX = Math.Max(maxX, currPosition.X);␍ maxY = Math.Max(maxY, currPosition.Y);␍ }␍ }␍ ␍ ␍ ␍ //generate complex field to ␍ ComplexField cfSplineInputParameter = new ComplexField(new Vector(listPositionsForInterpolation.Count, 3), false, PrecisionMode.Double);␍ for (int runLists = 0; runLists < listPositionsForInterpolation.Count; runLists++) {␍ cfSplineInputParameter[runLists, 0] = listPositionsForInterpolation[runLists].X;␍ cfSplineInputParameter[runLists, 1] = listPositionsForInterpolation[runLists].Y;␍ cfSplineInputParameter[runLists, 2] = listPhaseValuesForInterpolation[runLists];␍ }␍ //generate spline interpolation object for initial input parameter␍ SplineInterpolationGenerator splineGenerator = new SplineInterpolationGenerator(cfSplineInputParameter, Vector3D.UndefinedVector, 3, 3, false);␍ SplineInterpolation splineInterpolationInputWavefront = splineGenerator.GenerateSplineInterpolationObject();␍ ␍ double samplingDistanceRadial = 100e-9;␍ int numberPoints = (int)(maxX / samplingDistanceRadial);␍ if (numberPoints % 2 == 0) {␍ numberPoints++;␍ }␍ ␍ CFieldDerivative1DReal cfRadialPhaseValue = new CFieldDerivative1DReal(numberPoints, double.NaN);␍ CFieldDerivative1DReal cfRadialPhaseGradient = new CFieldDerivative1DReal(numberPoints, double.NaN);␍ ␍ //loop over all points␍ for (int runPoints = 0; runPoints < numberPoints; runPoints++) {␍ //define radial position␍ double currR = runPoints * samplingDistanceRadial;␍ //evaluate phase ␍ double evaluatedRadialPhaseValue = splineInterpolationInputWavefront.PointInterpolate(new VectorD(currR, 0));␍ double evaluatedRadialPhaseGradient = splineInterpolationInputWavefront.PointInterpolationGradient(new VectorD(currR, 0)).X;␍ //store values in data matrix␍ cfRadialPhaseValue[runPoints] = evaluatedRadialPhaseValue;␍ cfRadialPhaseGradient[runPoints] = evaluatedRadialPhaseGradient;␍ }␍ ␍ //generate 1D data array ␍ DataArray1D daRadialPhaseResponseInput = new DataArray1D(new ComplexField1DArray(new ComplexField[] { cfRadialPhaseValue, cfRadialPhaseGradient }),␍ new PhysicalProperty[] { PhysicalProperty.AngleRad, PhysicalProperty.NoUnit },␍ new string[] { "Interpolated Input Phase", "Interpolated Input Phase Gradient" },␍ samplingDistanceRadial,␍ 0,␍ PhysicalProperty.Length,␍ "Radius");␍ ␍ return daRadialPhaseResponseInput;␍ ␍ }␍ {}␀Generate Radial Wavefront Phase Info4/7/2026 5:26:55 PMAttribution␣1false␀1