From 66f368947bcbc2ff8104d382dab78229ca35de13 Mon Sep 17 00:00:00 2001
From: Erich Eckner <git@eckner.net>
Date: Tue, 23 Sep 2014 16:48:45 +0200
Subject: auch reflektierten Puls anhand von Trajektorie rekonstruieren ...
 bisher unvollständig
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 ROM.lpr | 33 ++++++++++++++++++++-------------
 1 file changed, 20 insertions(+), 13 deletions(-)

(limited to 'ROM.lpr')

diff --git a/ROM.lpr b/ROM.lpr
index 5863de0..26aebb3 100644
--- a/ROM.lpr
+++ b/ROM.lpr
@@ -12,17 +12,17 @@ uses
   { you can add units after this },
   SysUtils,ROMunit, mathunit, Math;
 
-var inPuls,refPuls,surTraj: tExtPointArray;
-    i,smooth:               longint;
-    tmax,wmax,absShift:     extended;
-    force,fourier:          boolean;
+var inPulsO,inPuls,refPuls,surTraj,cRefPuls: tExtPointArray;
+    i,smooth:                                longint;
+    tmax,wmax,absShift:                      extended;
+    force,fourier:                           boolean;
 
-const Verwendung='Verwendung: ROM ($Einfallspuls_Datei $Ausfallspuls_Datei)/(- $trace-Datei-Prefix) $output_inPuls $output_refPuls $output_Trajektorie '+
+const Verwendung='Verwendung: ROM ($Einfallspuls_Datei $Ausfallspuls_Datei)/(- $trace-Datei-Prefix) $output_inPuls $output_refPuls $output_Trajektorie $output_cRefPuls '+
    '[-s/--smooth $n] [-f/--force] [-t/--tmax $t] [-w/--wmax $w] [-F/--FFT] [-d/--dt $dt]';
 
 begin
-  if paramcount<5 then Fehler(Verwendung);
-  i:=6;
+  if paramcount<6 then Fehler(Verwendung);
+  i:=7;
   force:=false;
   smooth:=1;
   tmax:=-1;
@@ -79,23 +79,23 @@ begin
     Fehler(Verwendung);
   end;
   if paramstr(1)='-' then
-    readRawInputs(paramstr(2),inPuls,refPuls,absShift)
+    readRawInputs(paramstr(2),inPulsO,refPuls,absShift)
   else begin
-    readTextInput(paramstr(1),inPuls);
+    readTextInput(paramstr(1),inPulsO);
     readTextInput(paramstr(2),refPuls);
   end;
   write(stderr,'Input sortieren ...');
-  sort(inPuls);
+  sort(inPulsO);
   sort(refPuls);
   writeln(stderr,' fertig');
-  uniq(inPuls,false);
+  uniq(inPulsO,false);
   uniq(refPuls,false);
   write(stderr,'Input interpolieren ...');
-  interpoliere(inPuls);
+  interpoliere(inPulsO);
   interpoliere(refPuls);
   writeln(stderr,' fertig');
   flip(refPuls);
-  integrate(inPuls);
+  integrate(inPulsO,inPuls);
   integrate(refPuls);
   removeLinearOffset(inPuls);
   removeLinearOffset(refPuls);
@@ -142,9 +142,16 @@ begin
     normiere(refPuls);
     normiere(surTraj);
     writeln(stderr,' fertig');
+  end
+  else begin
+    write(stderr,'Reflektierten Puls berechnen ...');
+    berechneRefPuls(inPulsO,surTraj,cRefPuls);
+    integrate(cRefPuls);
+    writeln(stderr,' fertig');
   end;
   if paramstr(3)<>'-' then writeOutput(paramstr(3),inPuls);
   if paramstr(4)<>'-' then writeOutput(paramstr(4),refPuls);
   if paramstr(5)<>'-' then writeOutput(paramstr(5),surTraj);
+  if (paramstr(6)<>'-') and not fourier then writeOutput(paramstr(6),cRefPuls);
 end.
 
-- 
cgit v1.2.3-54-g00ecf