GPIAG-Software issueshttps://git.scc.kit.edu/groups/GPIAG-Software/-/issues2017-08-23T14:07:36+02:00https://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/18Taper for stf inversion2017-08-23T14:07:36+02:00niklas.thielTaper for stf inversionUsing "TAPER_STF", the length of the taper is calculated by using the main frequency given in the source file. If the filter frequency of the seismograms is significantly lower than the main frequency, the wavelet is cut at the end. This...Using "TAPER_STF", the length of the taper is calculated by using the main frequency given in the source file. If the filter frequency of the seismograms is significantly lower than the main frequency, the wavelet is cut at the end. This leads to a high frequency component in the synthetic seismograms.
WORKAROUND:
Use the taper of the stf instead of "TAPER_STF". This taper is specified in the string of "PARA" ("irtap=...") and the length of the taper is given by hand. More information are given in the manual.
In the future, "TAPER_STF" will be enhanced (by using the filter frequency for calculating the taper length) or removed.https://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/17string length model2017-04-21T14:30:59+02:00niklas.thielstring length modelIf the string length of e.g. the given starting model is too long an error occurs during reading. The error is unspecific for the user: "BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES [...] EXIT CODE: 11 [...]"
Suggestion: Either w...If the string length of e.g. the given starting model is too long an error occurs during reading. The error is unspecific for the user: "BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES [...] EXIT CODE: 11 [...]"
Suggestion: Either we intercept too long strings in the code and provide the user with a more meaningful error message or we start a new chapter in the manual with "common errors".https://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/16Write filtered data2017-04-02T13:32:26+02:00Florian Wittkampflorian.wittkamp@kit.eduWrite filtered dataNew behaviour in branch `feature/WriteData`.
### WRITE_FILTERED_DATA==1:
- Write synthetic data directly after forward modelling: `%s_p.su.shot%d.it%d`
- Write observed data after application of time domain filter: `%s_vx.su.obs.filter...New behaviour in branch `feature/WriteData`.
### WRITE_FILTERED_DATA==1:
- Write synthetic data directly after forward modelling: `%s_p.su.shot%d.it%d`
- Write observed data after application of time domain filter: `%s_vx.su.obs.filtered.shot%d.it%d`
### WRITE_FILTERED_DATA==2:
- Write synthetic data directly after forward modelling: `%s_p.su.shot%d.it%d`
- Write observed data after application of filter, window and integration: `%s_vx.su.obs.filtered.shot%d.it%d`
- Write synthetic data after application of filter, window and integration: `%s_vx.su.syn.filtered.shot%d.it%d`
@lgass @sb5563https://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/10Pseudo Hessian2017-09-17T14:23:24+02:00Simone ButzerPseudo HessianCalculation of Pseudo Hessian for preconditioning (see IFOS2D) as alternative to the diagonal Hessian approximationCalculation of Pseudo Hessian for preconditioning (see IFOS2D) as alternative to the diagonal Hessian approximationcode maintenancehttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/9Automatic calculation of water level for Hessian preconditioning2017-09-17T14:23:24+02:00Simone ButzerAutomatic calculation of water level for Hessian preconditioningcode maintenancehttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/8Inversion of density for Hessian and L-BFGS2017-09-17T14:23:24+02:00Simone ButzerInversion of density for Hessian and L-BFGS- test inversion of density with preconditioning with the diagonal Hessian
- include inversion of density in the L-BFGS sheme and test for multi-parameter example- test inversion of density with preconditioning with the diagonal Hessian
- include inversion of density in the L-BFGS sheme and test for multi-parameter examplecode maintenancehttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/7rearange and cleanup the terminal output/.out files2017-09-17T14:23:24+02:00Tilman Steinwegrearange and cleanup the terminal output/.out filescode maintenanceTilman SteinwegTilman Steinweghttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/6seismogramm output2017-09-17T14:23:24+02:00Tilman Steinwegseismogramm output- The Seismogramm header has to be optimized in the su case,
- segy outputs needs more testing, maybe reduce su (liitle endian), binary and ascii output
would be a good idea, only ascii and su input is implemented jet anyway.
- overw...- The Seismogramm header has to be optimized in the su case,
- segy outputs needs more testing, maybe reduce su (liitle endian), binary and ascii output
would be a good idea, only ascii and su input is implemented jet anyway.
- overwriting existing su files doesn't work right it all cases
- merging the seismogram in the code would give a speedup (less data in/output) / or reactivate seismerge program to handle this issue after the inversion.code maintenanceTilman SteinwegTilman Steinweghttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/5additional misfit definitions2017-09-17T14:23:24+02:00Tilman Steinwegadditional misfit definitionspreparation of the code for field data exampleshttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/4Adding mor than one relaxatation frequency for viscoelastic modelling2017-09-17T14:23:24+02:00Tilman SteinwegAdding mor than one relaxatation frequency for viscoelastic modellingTilman SteinwegTilman Steinweghttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/3source time function inversion2017-09-17T14:23:24+02:00Tilman Steinwegsource time function inversionFor field data an inversion for the source time function
is usefull and maybe necessary.For field data an inversion for the source time function
is usefull and maybe necessary.preparation of the code for field data exampleshttps://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/2Parallel Shot calculation2017-09-17T14:23:24+02:00Tilman SteinwegParallel Shot calculationIn 3D the point where the domain decomposition
reaches a speedup not proportional to the amount of processing units
can be reached if the domains become small.
Parallel calculation of the shots will increase the number of processing u...In 3D the point where the domain decomposition
reaches a speedup not proportional to the amount of processing units
can be reached if the domains become small.
Parallel calculation of the shots will increase the number of processing units which can be used
and provide a linear speedup. But more memory has to be allocated (proportional to the amount of parallel calculated shots). https://git.scc.kit.edu/GPIAG-Software/IFOS3D/-/issues/1switch coordinate input2017-09-17T14:23:24+02:00Tilman Steinwegswitch coordinate inputAt the moment z denotes the depth in model and seismogram(header) input/output
and the input files. In the code "y" is used as depth.
y should be used therefore as depth again to be consistent.
So the z<->y transformation is no longe...At the moment z denotes the depth in model and seismogram(header) input/output
and the input files. In the code "y" is used as depth.
y should be used therefore as depth again to be consistent.
So the z<->y transformation is no longer needed.code maintenancehttps://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/14Source wavelet correction filter do not work with frequency filtering2017-04-02T13:32:26+02:00Florian Wittkampflorian.wittkamp@student.kit.eduSource wavelet correction filter do not work with frequency filteringThe calculation of the source wavelet correction filter is not working with frequency filtering.
The produced synthetic data will have full frequency content, the filtering is not working here.
We try to fix this in branch IFOS2D_ST...The calculation of the source wavelet correction filter is not working with frequency filtering.
The produced synthetic data will have full frequency content, the filtering is not working here.
We try to fix this in branch IFOS2D_STFfilter.
Notify: @lgass Florian Wittkampflorian.wittkamp@student.kit.eduFlorian Wittkampflorian.wittkamp@student.kit.eduhttps://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/13Zero-Phase filtering is not working2017-04-02T13:32:26+02:00Florian Wittkampflorian.wittkamp@student.kit.eduZero-Phase filtering is not workingThe zero-phase option within IFOS is not working probably.
The filter produces artifacts and is not zero-phase at all.
Sometimes it produces a high group delay.
Notify: @lgass The zero-phase option within IFOS is not working probably.
The filter produces artifacts and is not zero-phase at all.
Sometimes it produces a high group delay.
Notify: @lgass Florian Wittkampflorian.wittkamp@student.kit.eduFlorian Wittkampflorian.wittkamp@student.kit.eduhttps://git.scc.kit.edu/GPIAG-Software/SOFI2D/-/issues/1implement double-couple source2017-07-25T11:47:43+02:00zf8974implement double-couple sourceA double-couple source as in SOFI3D is not yet available. It could be transferred from SOFI3D/src/eqsource.c.A double-couple source as in SOFI3D is not yet available. It could be transferred from SOFI3D/src/eqsource.c.https://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/12Buffer exchange2017-07-25T18:11:08+02:00Florian Wittkampflorian.wittkamp@student.kit.eduBuffer exchangeIn branch development commit 8d859168df8b1100cab15d96f2b366ab9341f88f from @tmetz in branch master is reverted due to merge of branch from @wittkamp.
I will readd the commit ASAP.In branch development commit 8d859168df8b1100cab15d96f2b366ab9341f88f from @tmetz in branch master is reverted due to merge of branch from @wittkamp.
I will readd the commit ASAP.Florian Wittkampflorian.wittkamp@student.kit.eduFlorian Wittkampflorian.wittkamp@student.kit.eduhttps://git.scc.kit.edu/GPIAG-Software/SOFI3D/-/issues/2implemention of additional relaxation mechanisms2017-07-24T07:52:56+02:00Tilman Steinwegimplemention of additional relaxation mechanismsAt the current state of the software only L=1 is implemented.
Implementation of L>1 can be imitated from sofi2D.
Maybe an f4tensor is needed to cover all 3 dimensions and L mechanisms.
Performance has to be checked after the implement...At the current state of the software only L=1 is implemented.
Implementation of L>1 can be imitated from sofi2D.
Maybe an f4tensor is needed to cover all 3 dimensions and L mechanisms.
Performance has to be checked after the implementation (compared to old L=1).Tilman SteinwegTilman Steinweghttps://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/11Clean up MAKEFILE etc.2017-04-02T13:32:26+02:00Florian Wittkampflorian.wittkamp@student.kit.eduClean up MAKEFILE etc.Cleanup of all smaller files like Makefile and remove old stuff.Cleanup of all smaller files like Makefile and remove old stuff.Cleanup DENISEFlorian Wittkampflorian.wittkamp@student.kit.eduFlorian Wittkampflorian.wittkamp@student.kit.eduhttps://git.scc.kit.edu/GPIAG-Software/IFOS2D/-/issues/10wrong order of loops2017-04-02T13:32:26+02:00zf8974wrong order of loops- loops e.g. over gradients are ineffective- loops e.g. over gradients are ineffectiveCleanup DENISEFlorian Wittkampflorian.wittkamp@student.kit.eduFlorian Wittkampflorian.wittkamp@student.kit.edu