Setup, logs, etc.; Goodbye

log/2023_06_21_zVp15/RunMpp.sh

+#!/bin/sh
+### Job name
+#SBATCH --job-name=MppFWILennart
+
+### File for the output
+##vorher#SBATCH --output=Mpp%A   ???????
+#SBATCH --output=Mpp%A
+
+### Time your job needs to execute
+#SBATCH --time=280:00:00
+
+### Use more than one node for parallel jobs on distributed-memory systems
+###CR Version war 2, 4 geht, sofern genug frei.
+#SBATCH --nodes=8
+###SBATCH --nodes=1
+
+### Number of CPUS per task
+#SBATCH --cpus-per-task=1
+
+### Disable hyperthreading by setting the tasks per core to 1
+#SBATCH --ntasks-per-core=1
+
+#SBATCH --mail-type=ALL
+
+#SBATCH --mail-user=uuwyh@student.kit.edu
+
+### Number of processes per node, in our case
+#SBATCH --ntasks-per-node=16
+mpirun M++

log/2023_06_21_zVp15/testconf.conf

+# Regular Expression um aus Test_Config von Clion eine ... M++ Config zu machen
+# \{"([A-Za-z0-9_ ]+)", +"([A-Za-z0-9_ ,\.\-]+)"}.
+
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, Materialpunkt, level 1 (Ergebnisse ganz gut)
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, TwoMedium, level 1
+#Back to FastMode (2 Quellen, 2 Iterationen, Pkt.Material, level 0)
+#21.06 Receiver wieder auf 64. Bei multi noch nicht, jetzt aber bei rho rekonstruktion
+
+   #Problem= VAcousticWaveTwoMedium;//<----- Geometrie oder so Siehe: src/Problems/FWIInversionProblems.cpp
+   Problem= VAcousticWaveC2Benchmark;
+   FWISO_Steps = 15;   #Anzahl Iterationen
+   Monoparam= Vp; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   plot_tStep= -1;
+   level = 3;
+   T=0.6;
+   dt= 0.001;
+   FWISourceCount0=6;   #Anzahl Quellen
+
+   FWIReconstructRho= 0;
+   FWIReconstructV_P= 1;
+   FWIReconstructTauP= 0;
+
+   #Schnelle Configuration
+   #Problem= VAcousticWaveTwoMedium;//<----- Geometrie oder so Siehe: src/Problems/FWIInversionProblems.cpp
+   #Problem= VAcousticWaveC2Benchmark;
+   #FWISO_Steps = 2;//scheint zu klappen!
+   #Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   #plot_tStep= -1;
+   #level = 2;
+   #T=0.6;
+   #dt= 0.002;
+
+
+   ConfigVerbose=2;
+   ForwardSolver= VADGTS;
+   Distribution= RCB;
+   Model = FWIInverse_DGTS;
+   Mesh= SeismogramGrid_test;
+
+   useTaper=1;
+   taperRadiusSource=0;
+   taperRadiusReceiver=0;
+   taperExponent=-0.1;
+   taperIgnoreReceiverCells=1;
+   taperIgnoreSourceCells=1;
+   #Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   ReconstructionType=SO;#!!!
+   FWIReceiverWaypoints= 2;
+   FWIReceiverWaypoint0= 1350.0, 370.0;
+   FWIReceiverWaypoint1= 1350.0, 1630.0;
+   #FWIReceiverWaypoint0= 1600.0, 370.0;#Selbst weiter nach "rechts" geschoben
+   #FWIReceiverWaypoint1= 1600.0, 1630.0;
+   FWIReceiverIncludeLastWaypoint= 1;
+   FWIReceiverCount0= 64;//Anzahl Receiver
+   UseSparseMatrix=1;
+   Qfactor=100;
+   SignalName=RickerDerivative;
+   SignalDelay=0;
+   SignalCutoff=0.8;
+   SignalFreq=25;
+   numL=0;#!!!
+   FWISourceWaypoints=2;
+   FWISourceWaypoint0=600,370;
+   FWISourceWaypoint1=600,1620;
+   #FWISourceCount0=2;   #Anzahl Quelle
+   #FWISourceCount0=2;
+   FWIReceiverIncludeLastWaypoint=1;
+   measurement_space= point;
+   measurement_time= point;
+   compMeasure= 0;
+   source_component= 2;
+   perturbation_comp= 2;
+   dampingFlux= 1;
+   allowPlotting_M=1;
+   allowPlotting_WF=0;
+   plotfactor= 10;
+   parametrization= taup;#!!!!
+   #T=0.6;
+   #dt= 0.002;
+   deg= 1;
+   #level = 0;#Feinheit des Gitters (1 mehr -> Halbierung)//Achtung: 3 Böse, <2 Gut
+   TimeVerbose= 1;
+   rkorder= -2;
+   #plot_tStep= -1;
+   #FWISO_Steps = 2;//scheint zu klappen!
+   Steps = 2;//scheint auch nichts zu bringen :(
+
+   Kmax=150;
+   Keps= 1e-07;
+   Overlap= dG1;
+   debugAdjoint= 0;
+   debugFrechet= 1;
+   source_factor= 242574257420000.5742;
+   LinearSolver= LS;
+   LinearVerbose= -1;
+   LinearSteps=100;
+   LinearNthIterOutput= 1;
+   LinearEpsilon= 1e-12;
+   LinearReduction= 0.0001;
+   FWINonlinearSolver= FWISO;#!!!
+   FWI_DebugPlottingLevel= 1;
+   CGREGINNCanonic_Steps=3;

log/2023_06_27_IC_RhoMono_HighResolution/0testconf.conf

+# Regular Expression um aus Test_Config von Clion eine ... M++ Config zu machen
+# \{"([A-Za-z0-9_ ]+)", +"([A-Za-z0-9_ ,\.\-]+)"}.
+
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, Materialpunkt, level 1 (Ergebnisse ganz gut)
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, TwoMedium, level 1
+#Back to FastMode (2 Quellen, 2 Iterationen, Pkt.Material, level 0)
+#21.06 Receiver wieder auf 64. Bei multi noch nicht, jetzt aber bei rho rekonstruktion
+#Neuer größerer Durchlauf 26.Juni Ohne Inverse Crime für Vp
+
+   #Problem= VAcousticWaveTwoMedium;//<----- Geometrie oder so Siehe: src/Problems/FWIInversionProblems.cpp
+   Problem= VAcousticWaveC2Benchmark;
+   FWISO_Steps = 15;   #Anzahl Iterationen
+   Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   plot_tStep= -1;
+   level = 3;
+   T=0.6;
+   dt= 0.0005;
+   FWISourceCount0=6;   #Anzahl Quellen, Empfänger default 64
+   deg = 2;    #In paper CRC deg = 2. Vor 26.Juni2023 wars 1
+   FWIReconstructRho= 1;
+   FWIReconstructV_P= 0;
+   FWIReconstructTauP= 0;
+   numL=0;#!!!
+
+   #seismopath = C2NoDampingT06dt0005/C2NoDampingT06dt0005;#im Build Ordner!
+   #seismopath = readmePlease;//klappt in der ersten Iteration mit ..Please0
+   #seismopath = ../FWI/conf/C2NoDampingT06dt0005/C2NoDampingT06dt0005;//mal sehen obs klappt
+   #seismopath = /home/lionbart/fwi/FWI/conf/C2NoDampingT06dt0005/C2NoDampingT06dt0005;//OUT OF RANGE ERROR
+
+   #ReadFormat = FWI;#IFOS ist falsches Format.
+   #SpecfileNrOfSamples = 6;
+
+   SaveObservedSeismogram = 0;#Plottet "d_obs", das Vorwärtsgelöste Wellenfeld (mit echten Materialparam)
+
+   SignalName=Ricker;#RickerDerivative vorher!!! Geändert 26.Juni.2023
+
+
+
+   ConfigVerbose=2;
+   ForwardSolver= VADGTS;
+   Distribution= RCB;
+   Model = FWIInverse_DGTS;
+   Mesh= SeismogramGrid_test;
+
+   useTaper=1;
+   taperRadiusSource=0;
+   taperRadiusReceiver=0;
+   taperExponent=-0.1;
+   taperIgnoreReceiverCells=1;
+   taperIgnoreSourceCells=1;
+   #Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   ReconstructionType=SO;#!!!
+   FWIReceiverWaypoints= 2;
+   FWIReceiverWaypoint0= 1350.0, 370.0;
+   FWIReceiverWaypoint1= 1350.0, 1630.0;
+   FWIReceiverIncludeLastWaypoint= 1;
+   FWIReceiverCount0= 64;//Anzahl Receiver
+   UseSparseMatrix=1;
+   Qfactor=100;
+   #SignalName=Ricker;# RickerDerivative vorher!!! Geändert 26.Juni.2023
+   SignalDelay=0;
+   SignalCutoff=0.8;
+   SignalFreq=25;
+   #numL=0;#!!!
+   FWISourceWaypoints=2;
+   FWISourceWaypoint0=600,370;
+   FWISourceWaypoint1=600,1620;
+   #FWISourceCount0=2;   #Anzahl Quelle
+   #FWISourceCount0=2;
+   FWIReceiverIncludeLastWaypoint=1;
+   measurement_space= point;
+   measurement_time= point;
+   compMeasure= 0;
+   source_component= 2;
+   perturbation_comp= 2;
+   dampingFlux= 1;
+   allowPlotting_M=1;
+   allowPlotting_WF=0;
+   plotfactor= 10;
+   parametrization= taup;#!!!!
+   #T=0.6;
+   #dt= 0.002;
+   #deg= 1;
+   #level = 0;#Feinheit des Gitters (1 mehr -> Halbierung)//Achtung: 3 Böse, <2 Gut
+   TimeVerbose= 1;
+   rkorder= -2;
+   #plot_tStep= -1;
+   #FWISO_Steps = 2;//scheint zu klappen!
+   Steps = 2;//scheint auch nichts zu bringen :(
+
+   Kmax=150;
+   Keps= 1e-07;
+   Overlap= dG1;
+   debugAdjoint= 0;
+   debugFrechet= 1;
+   source_factor= 242574257420000.5742;
+   LinearSolver= LS;
+   LinearVerbose= -1;
+   LinearSteps=100;
+   LinearNthIterOutput= 1;
+   LinearEpsilon= 1e-12;
+   LinearReduction= 0.0001;
+   FWINonlinearSolver= FWISO;#!!!
+   FWI_DebugPlottingLevel= 1;
+   CGREGINNCanonic_Steps=3;

log/2023_06_27_noIC_VpMono_HighResolution/RunMpp.sh

+#!/bin/sh
+### Job name
+#SBATCH --job-name=MppFWILennart
+
+### File for the output
+##vorher#SBATCH --output=Mpp%A   ???????
+#SBATCH --output=Mpp%A
+
+### Time your job needs to execute
+#SBATCH --time=200:00:00
+
+### Use more than one node for parallel jobs on distributed-memory systems
+###CR Version war 2, 4 geht, sofern genug frei.
+#SBATCH --nodes=8
+###SBATCH --nodes=1
+
+### Number of CPUS per task
+#SBATCH --cpus-per-task=1
+
+### Disable hyperthreading by setting the tasks per core to 1
+#SBATCH --ntasks-per-core=1
+
+#SBATCH --mail-type=ALL
+
+#SBATCH --mail-user=uuwyh@student.kit.edu
+
+### Number of processes per node, in our case
+#SBATCH --ntasks-per-node=16
+mpirun M++

log/2023_06_27_noIC_VpMono_HighResolution/testconf.conf

+# Regular Expression um aus Test_Config von Clion eine ... M++ Config zu machen
+# \{"([A-Za-z0-9_ ]+)", +"([A-Za-z0-9_ ,\.\-]+)"}.
+
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, Materialpunkt, level 1 (Ergebnisse ganz gut)
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, TwoMedium, level 1
+#Back to FastMode (2 Quellen, 2 Iterationen, Pkt.Material, level 0)
+#21.06 Receiver wieder auf 64. Bei multi noch nicht, jetzt aber bei rho rekonstruktion
+#Neuer größerer Durchlauf 26.Juni Ohne Inverse Crime für Vp; Neu: SignalName: Ricker; deg = 2; dt = 0.0005
+
+   #Problem= VAcousticWaveTwoMedium;//<----- Geometrie oder so Siehe: src/Problems/FWIInversionProblems.cpp
+   Problem= VAcousticWaveC2Benchmark;
+   FWISO_Steps = 15;   #Anzahl Iterationen
+   Monoparam= Vp; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   plot_tStep= -1;
+   level = 3;
+   T=0.6;
+   dt= 0.0005;
+   FWISourceCount0=6;   #Anzahl Quellen 6, Empfänger default 64 unten
+   deg = 2;    #In paper CRC deg = 2. Vor 26.Juni2023 wars 1
+   FWIReconstructRho= 0;
+   FWIReconstructV_P= 1;
+   FWIReconstructTauP= 0;
+   numL=0;#!!!
+
+   seismopath = C2NoDampingT06dt0005/C2NoDampingT06dt0005;#im Build Ordner!
+   #seismopath = readmePlease;//klappt in der ersten Iteration mit ..Please0
+   #seismopath = ../FWI/conf/C2NoDampingT06dt0005/C2NoDampingT06dt0005;//mal sehen obs klappt
+   #seismopath = /home/lionbart/fwi/FWI/conf/C2NoDampingT06dt0005/C2NoDampingT06dt0005;//OUT OF RANGE ERROR
+
+   ReadFormat = FWI;#IFOS ist falsches Format.
+   SpecfileNrOfSamples = 6;
+
+   SaveObservedSeismogram = 0;#Plottet "d_obs", das Vorwärtsgelöste Wellenfeld (mit echten Materialparam)
+
+   SignalName=Ricker;#RickerDerivative vorher!!! Geändert 26.Juni.2023
+
+
+
+   ConfigVerbose=2;
+   ForwardSolver= VADGTS;
+   Distribution= RCB;
+   Model = FWIInverse_DGTS;
+   Mesh= SeismogramGrid_test;
+
+   useTaper=1;
+   taperRadiusSource=0;
+   taperRadiusReceiver=0;
+   taperExponent=-0.1;
+   taperIgnoreReceiverCells=1;
+   taperIgnoreSourceCells=1;
+   #Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   ReconstructionType=SO;#!!!
+   FWIReceiverWaypoints= 2;
+   FWIReceiverWaypoint0= 1350.0, 370.0;
+   FWIReceiverWaypoint1= 1350.0, 1630.0;
+   FWIReceiverIncludeLastWaypoint= 1;
+   FWIReceiverCount0= 64;//Anzahl Receiver
+   UseSparseMatrix=1;
+   Qfactor=100;
+   #SignalName=Ricker;# RickerDerivative vorher!!! Geändert 26.Juni.2023
+   SignalDelay=0;
+   SignalCutoff=0.8;
+   SignalFreq=25;
+   #numL=0;#!!!
+   FWISourceWaypoints=2;
+   FWISourceWaypoint0=600,370;
+   FWISourceWaypoint1=600,1620;
+   #FWISourceCount0=2;   #Anzahl Quelle
+   #FWISourceCount0=2;
+   FWIReceiverIncludeLastWaypoint=1;
+   measurement_space= point;
+   measurement_time= point;
+   compMeasure= 0;
+   source_component= 2;
+   perturbation_comp= 2;
+   dampingFlux= 1;
+   allowPlotting_M=1;
+   allowPlotting_WF=0;
+   plotfactor= 10;
+   parametrization= taup;#!!!!
+   #T=0.6;
+   #dt= 0.002;
+   #deg= 1;
+   #level = 0;#Feinheit des Gitters (1 mehr -> Halbierung)//Achtung: 3 Böse, <2 Gut
+   TimeVerbose= 1;
+   rkorder= -2;
+   #plot_tStep= -1;
+   #FWISO_Steps = 2;//scheint zu klappen!
+   Steps = 2;//scheint auch nichts zu bringen :(
+
+   Kmax=150;
+   Keps= 1e-07;
+   Overlap= dG1;
+   debugAdjoint= 0;
+   debugFrechet= 1;
+   source_factor= 242574257420000.5742;
+   LinearSolver= LS;
+   LinearVerbose= -1;
+   LinearSteps=100;
+   LinearNthIterOutput= 1;
+   LinearEpsilon= 1e-12;
+   LinearReduction= 0.0001;
+   FWINonlinearSolver= FWISO;#!!!
+   FWI_DebugPlottingLevel= 1;
+   CGREGINNCanonic_Steps=3;

log/2023_06_30_MultiHighRes_SO/RunMpp.sh

+#!/bin/sh
+### Job name
+#SBATCH --job-name=MppFWILennart
+
+### File for the output
+##vorher#SBATCH --output=Mpp%A   ???????
+#SBATCH --output=Mpp%A
+
+### Time your job needs to execute
+#SBATCH --time=200:00:00
+
+### Use more than one node for parallel jobs on distributed-memory systems
+###CR Version war 2, 4 geht, sofern genug frei.
+#SBATCH --nodes=8
+###SBATCH --nodes=1
+
+### Number of CPUS per task
+#SBATCH --cpus-per-task=1
+
+### Disable hyperthreading by setting the tasks per core to 1
+#SBATCH --ntasks-per-core=1
+
+#SBATCH --mail-type=ALL
+
+#SBATCH --mail-user=uuwyh@student.kit.edu
+
+### Number of processes per node, in our case
+#SBATCH --ntasks-per-node=16
+mpirun M++

log/2023_06_30_MultiHighRes_SO/testconf.conf

+# Regular Expression um aus Test_Config von Clion eine ... M++ Config zu machen
+# \{"([A-Za-z0-9_ ]+)", +"([A-Za-z0-9_ ,\.\-]+)"}.
+
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, Materialpunkt, level 1 (Ergebnisse ganz gut)
+#Durchlauf 2023Mai05 mit 4 Quellen, 8 Iterationen, TwoMedium, level 1
+#Back to FastMode (2 Quellen, 2 Iterationen, Pkt.Material, level 0)
+#21.06 Receiver wieder auf 64. Bei multi noch nicht, jetzt aber bei rho rekonstruktion
+#Neuer größerer Durchlauf 26.Juni Ohne Inverse Crime für Vp
+
+   #ReconstructionType=Classical;
+   #FWINonlinearSolver= FWI_mSDREGINN;
+
+
+   #Problem= VAcousticWaveTwoMedium;//<----- Geometrie oder so Siehe: src/Problems/FWIInversionProblems.cpp
+   Problem= VAcousticWaveC2Benchmark;
+   FWISO_Steps = 10;   #Anzahl Iterationen 
+   Monoparam= multi; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   plot_tStep= -1;
+   level = 3;#3 nicht auf Laptop starten, sonst stürzt ab
+   T=0.6;
+   dt= 0.0005;
+   FWISourceCount0=6;   #Anzahl Quellen, Empfänger default 64
+   deg = 2;    #In paper CRC deg = 2. Vor 26.Juni2023 wars 1
+   FWIReconstructRho= 1;
+   FWIReconstructV_P= 1;
+   FWIReconstructTauP= 0;
+
+   #seismopath = C2NoDampingT06dt0005/C2NoDampingT06dt0005;#im Build Ordner!
+   #ReadFormat = FWI;#IFOS ist falsches Format.
+   #SpecfileNrOfSamples = 6;
+
+   SaveObservedSeismogram = 0;#Plottet "d_obs", das Vorwärtsgelöste Wellenfeld (mit echten Materialparam)
+
+   SignalName=Ricker;#RickerDerivative vorher!!! Geändert 26.Juni.2023
+   ReconstructionType=SO;#!!!
+   FWINonlinearSolver= FWISO;#!!!
+   numL=0;#!!!
+
+
+
+   ConfigVerbose=2;
+   ForwardSolver= VADGTS;
+   Distribution= RCB;
+   Model = FWIInverse_DGTS;
+   Mesh= SeismogramGrid_test;
+
+   useTaper=1;
+   taperRadiusSource=0;
+   taperRadiusReceiver=0;
+   taperExponent=-0.1;
+   taperIgnoreReceiverCells=1;
+   taperIgnoreSourceCells=1;
+   #Monoparam= Rho; #can be "multi"/empty or "Vp", "Tau_p", "Rho"
+   #ReconstructionType=SO;#!!!
+   FWIReceiverWaypoints= 2;
+   FWIReceiverWaypoint0= 1350.0, 370.0;
+   FWIReceiverWaypoint1= 1350.0, 1630.0;
+   FWIReceiverIncludeLastWaypoint= 1;
+   FWIReceiverCount0= 64;//Anzahl Receiver
+   UseSparseMatrix=1;
+   Qfactor=100;
+   #SignalName=Ricker;# RickerDerivative vorher!!! Geändert 26.Juni.2023
+   SignalDelay=0;
+   SignalCutoff=0.8;
+   SignalFreq=25;
+   #numL=0;#!!!
+   FWISourceWaypoints=2;
+   FWISourceWaypoint0=600,370;
+   FWISourceWaypoint1=600,1620;
+   #FWISourceCount0=2;   #Anzahl Quelle
+   #FWISourceCount0=2;
+   FWIReceiverIncludeLastWaypoint=1;
+   measurement_space= point;
+   measurement_time= point;
+   compMeasure= 0;
+   source_component= 2;
+   perturbation_comp= 2;
+   dampingFlux= 1;
+   allowPlotting_M=1;
+   allowPlotting_WF=0;
+   plotfactor= 10;
+   parametrization= taup;#!!!!
+   #T=0.6;
+   #dt= 0.002;
+   #deg= 1;
+   #level = 0;#Feinheit des Gitters (1 mehr -> Halbierung)//Achtung: 3 Böse, <2 Gut
+   TimeVerbose= 1;
+   rkorder= -2;
+   #plot_tStep= -1;
+   #FWISO_Steps = 2;//scheint zu klappen!
+   Steps = 2;//scheint auch nichts zu bringen :(
+
+   Kmax=150;
+   Keps= 1e-07;
+   Overlap= dG1;
+   debugAdjoint= 0;
+   debugFrechet= 1;
+   source_factor= 242574257420000.5742;
+   LinearSolver= LS;#Könnte man auf GMRES (Schreibweise?) ändern..
+   LinearVerbose= -1;
+   LinearSteps=100;
+   LinearNthIterOutput= 1;
+   LinearEpsilon= 1e-12;
+   LinearReduction= 0.0001;
+   #FWINonlinearSolver= FWI_mSDREGINN;
+   FWI_DebugPlottingLevel= 1;
+   CGREGINNCanonic_Steps=3;
+

log/2023_07_06_Multi_FO_Iter10/RunMpp.sh

+#!/bin/sh
+### Job name
+#SBATCH --job-name=MppFWILennart
+
+### File for the output
+##vorher#SBATCH --output=Mpp%A   ???????
+#SBATCH --output=Mpp%A
+
+### Time your job needs to execute
+#SBATCH --time=222:00:00
+
+### Use more than one node for parallel jobs on distributed-memory systems
+###CR Version war 2, 4 geht, sofern genug frei.
+#SBATCH --nodes=8
+###SBATCH --nodes=1
+
+### Number of CPUS per task
+#SBATCH --cpus-per-task=1
+
+### Disable hyperthreading by setting the tasks per core to 1
+#SBATCH --ntasks-per-core=1
+
+#SBATCH --mail-type=ALL
+
+#SBATCH --mail-user=uuwyh@student.kit.edu
+
+### Number of processes per node, in our case
+#SBATCH --ntasks-per-node=16
+mpirun M++