geopsy.org

Hi,
Today I downloaded the new release (3.3.2) interested in verifying that SPAC & spac2disp worked (since this didn't happen with the previous release). I have verified with different databases that:
1) SPAC allows the setting but then does not start. The program stop, but does not crash and you have to force to close it.
2) spac2disp opens but after loading the target file a warning appears that says: "kmin is less than 1e-5". If you accept anyway and press enter, the autocorrelation curves appear, but the dispersion graph is empty.
3) gpfksimulator crash when you ask to save the results (cross section - after writing the output file)
4) the bug report cannot be sent yet because the usual warning appears saying that I have to upgrade to the latest release ...
This is the bug report
CrashReport>
<Error>
<severity>--FATAL ERROR- </severity>
<message>ASSERT: &quot;_ncol==m._nrow&quot; in file c:/Geopsy/distribs/src/3.3/geopsypack/QGpCoreMath/src/Matrix.h, line 136</message>
<time>2020-08-01 12:04:03.001</time>
<threads>main,layerPainter1,layerPainter2,layerPainter3,layerPainter4,layerPainter5,layerPainter6,layerPainter7,layerPainter8</threads>
<currentThread></currentThread>
<stack>---- thread main
double FKArrayMap::HRFKRayleigh(const QGpCoreMath::Point&amp;) const:848
void Simulator::on_crossSectionBut_clicked():407
---- thread layerPainter1
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter2
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter3
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter4
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter5
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter6
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter7
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter8
virtual void SciFigs::LayerPainterThread::run():59
</stack>
<bugStack>---- thread main
---- thread layerPainter1
---- thread layerPainter2
---- thread layerPainter3
---- thread layerPainter4
---- thread layerPainter5
---- thread layerPainter6
---- thread layerPainter7
---- thread layerPainter8
---- thread main
---- thread layerPainter1
---- thread layerPainter2
---- thread layerPainter3
---- thread layerPainter4
---- thread layerPainter5
---- thread layerPainter6
---- thread layerPainter7
---- thread layerPainter8
</bugStack>
</Error>
<application>gpfksimulator</application>
<version> Qt 5.12.3
QGpCoreTools 2.3.1 from 3.3.2
(commit b7d0c8ab6bb4726c71937b15816dcacae2ee5ea9)
QGpCoreMath 1.3.0 from 3.3.2
(commit 27776e1a7638081a41ef3d2815e998f8024e34e2)
DinverCore 1.2.1 from 3.3.2
(commit 7966e349ea0c51cc3b3fc3805c9faf24718b7d26)
QGpCoreWave 2.1.1 from 3.3.2
(commit 1818f01257daa38705071ccc310003c78a58aa22)
QGpCompatibility 1.1.0 from 3.3.2
(commit 2a1a1e452a5f36d25299613473280f699fc01e4f)
DinverDCCore 2.1.1 from 3.3.2
(commit e506d28eca9e687e4a4fe913816b19a2bc26c76b)
GeopsyCore 4.3.1 from 3.3.2
(commit 49e9ff90660eab0d1dbe7448d42476887e7fe319)
QGpCoreStat 1.1.0 from 3.3.2
(commit 48e42460af0d9650451b48ef167146c0a0bcccaf)
ArrayCore 4.0.2 from 3.3.2
(commit 4698505cfdde25441ce47254793337429dee9119)
QGpGuiTools 2.3.0 from 3.3.2
(commit e33f0db59e1efd148abc0d16c31d93bb5bf4f314)
SciFigs 3.2.0 from 3.3.2
(commit 44d8371e06a3e408858d30f6e430aea0e0e7e5a4)
QGpGuiMath 1.2.0 from 3.3.2
(commit 598890f255a708b64139acb2093a5557dd6214fc)
gpfksimulator 1.1.0 from 3.3.2
(commit a2e2a6a632b3758bd40b31153eebbd27509728fe)
geopsyfigs 1.0.1 from 3.3.2
(commit 9f7e51239f43119cc3f955538683224881c80599)
</version>
<system>windows 10</system>
</CrashReport>

Best regards
Luigi

Regards
Luigi

Hi Luigi,

Thanks for your reports.
1) I've just checked SPAC on 3.3.2 and it runs nicely. Is there any message in the log? Can you send a snapshot of the parameters (all tabs including rings)? There is maybe something linked to particular value.
2) If you increase kmin to something positive, the next time you restart spac2disp it will keep this value and it will not complain. If you want to get the dispersion plot without restarting, just modify a little bit the minimum velocity.
3) This tools is for experimental use only. It work fine for processing "FK Vertical". Cross-section works only for processing "HRFK Rayleigh". The reason is only a lack of time to implement it correctly for all other processing options. I disabled it when HRFKRayleigh is not selected. to avoid useless crashes.
4( Yes I know I have to repair this bug facility

Best regards,

Marc

Hi Marc,
Thanks for answer.
About the first point I tried to run SPAC using your database Lep_ring01.
The tool worked nicely with V data, but crash immediately with 3C data. I attach the screenshots for V and 3C parameters setting (I used exactly your rings dimensions) and, for the last one, the crash report too.
<CrashReport>
<Error>
<severity>--FATAL ERROR- </severity>
<message>ASSERT: &quot;_nrow==m._nrow &amp;&amp; _ncol==m._ncol&quot; in file c:/Geopsy/distribs/src/3.3/geopsypack/QGpCoreMath/src/Matrix.h, line 114</message>
<time>2020-08-04 11:34:00.377</time>
<threads>main,layerPainter1,layerPainter2,layerPainter3,layerPainter4,layerPainter5,layerPainter6,layerPainter7,layerPainter8,LoopWorker_0,LoopWorker_1,LoopWorker_2,LoopWorker_3</threads>
<currentThread></currentThread>
<stack>---- thread main
virtual void QGpCoreTools::ParallelLoop::start(int, int, bool):150
void QGpCoreTools::TaskManager::start(bool):75
virtual void GeopsyCore::AbstractTool::start():93
virtual void ArrayCore::AbstractArrayTool::start():203
virtual void SPACToolWidget::start():210
---- thread layerPainter1
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter2
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter3
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter4
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter5
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter6
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter7
virtual void SciFigs::LayerPainterThread::run():59
---- thread layerPainter8
virtual void SciFigs::LayerPainterThread::run():59
---- thread LoopWorker_0
bool ArrayCore::SPACCrossSpectrum::calculate(const QVector&lt;int&gt;&amp;):73
virtual void QGpCoreTools::LoopWorker::run():47
---- thread LoopWorker_1
bool QGpCoreTools::Cache::makeAvailable(const QGpCoreTools::CacheItem*):165
bool GeopsyCore::DoubleSignal::setValue(double):740
bool GeopsyCore::DoubleSignal::subtractValue(double):172
bool ArrayCore::SPACCrossSpectrum::calculate(const QVector&lt;int&gt;&amp;):73
virtual void QGpCoreTools::LoopWorker::run():47
---- thread LoopWorker_2
void GeopsyCore::TimeRangeList::add(double, const GeopsyCore::WindowingParameters&amp;, const GeopsyCore::SparseKeepSignal&amp;, const GeopsyCore::TimeRange&amp;):125
void GeopsyCore::TimeRangeList::add(double, const GeopsyCore::WindowingParameters&amp;, const GeopsyCore::SparseKeepSignal&amp;, const GeopsyCore::TimeRange&amp;, const GeopsyCore::SeismicEventTable*):294
virtual void QGpCoreTools::LoopWorker::run():47
---- thread LoopWorker_3
</stack>
<bugStack>---- thread main
---- thread layerPainter1
---- thread layerPainter2
---- thread layerPainter3
---- thread layerPainter4
---- thread layerPainter5
---- thread layerPainter6
---- thread layerPainter7
---- thread layerPainter8
---- thread LoopWorker_0
---- thread LoopWorker_1
---- thread LoopWorker_2
---- thread LoopWorker_3
---- thread main
---- thread layerPainter1
---- thread layerPainter2
---- thread layerPainter3
---- thread layerPainter4
---- thread layerPainter5
---- thread layerPainter6
---- thread layerPainter7
---- thread layerPainter8
---- thread LoopWorker_0
---- thread LoopWorker_1
---- thread LoopWorker_2
---- thread LoopWorker_3
</bugStack>
</Error>
<application>geopsy</application>
<version> Qt 5.12.3
QGpCoreTools 2.3.1 from 3.3.2
(commit b7d0c8ab6bb4726c71937b15816dcacae2ee5ea9)
QGpCoreMath 1.3.0 from 3.3.2
(commit 27776e1a7638081a41ef3d2815e998f8024e34e2)
GeopsyCore 4.3.1 from 3.3.2
(commit 49e9ff90660eab0d1dbe7448d42476887e7fe319)
QGpGuiTools 2.3.0 from 3.3.2
(commit e33f0db59e1efd148abc0d16c31d93bb5bf4f314)
SciFigs 3.2.0 from 3.3.2
(commit 44d8371e06a3e408858d30f6e430aea0e0e7e5a4)
QGpGuiMath 1.2.0 from 3.3.2
(commit 598890f255a708b64139acb2093a5557dd6214fc)
GeopsyGui 3.3.1 from 3.3.2
(commit 39ad97613d99f84d1b4bf7a9ece468668c6aff0b)
GeopsySLink 2.0.2 from 3.3.2
(commit 55c95aa444ec44fa21c3a0b298c4cea7044b0e10)
GeopsySLinkGui 2.0.2 from 3.3.2
(commit 1d7ec3e61b3b8ca7bfc6ca264d479ead63444996)
geopsy 3.3.1 from 3.3.2
(commit d9a7e21d6a9fc68e2364f63a91cc9c3ceecbe1a9)
geopsyfigs 1.0.1 from 3.3.2
(commit 9f7e51239f43119cc3f955538683224881c80599)
campbelltob3 1.0.1 from 3.3.2
(commit d27183a55c6eac7d8421708245705b326bb6e8d1)
DinverCore 1.2.1 from 3.3.2
(commit 7966e349ea0c51cc3b3fc3805c9faf24718b7d26)
QGpCoreWave 2.1.1 from 3.3.2
(commit 1818f01257daa38705071ccc310003c78a58aa22)
QGpCompatibility 1.1.0 from 3.3.2
(commit 2a1a1e452a5f36d25299613473280f699fc01e4f)
DinverDCCore 2.1.1 from 3.3.2
(commit e506d28eca9e687e4a4fe913816b19a2bc26c76b)
QGpCoreStat 1.1.0 from 3.3.2
(commit 48e42460af0d9650451b48ef167146c0a0bcccaf)
ArrayCore 4.0.2 from 3.3.2
(commit 4698505cfdde25441ce47254793337429dee9119)
QGpGuiWave 1.1.0 from 3.3.2
(commit d583719a3637376c7c7768a28327f0eb41314a6a)
ArrayGui 1.2.1 from 3.3.2
(commit 6c77f5a2a0f6ac8e78b16116ecca9ac3c2927291)
geopsyarray 3.2.1 from 3.3.2
(commit d438abdae142857b27ee9f1732fec9a0c7d58b1c)
geopsydamping 2.1.1 from 3.3.2
(commit 69d1c7368e67ed6cae9ad5cdb77e91efc45e6b81)
HVCore 1.0.1 from 3.3.2
(commit c0b656aa841a043c00bdbaf70b4775d0b2821159)
HVGui 1.0.0 from 3.3.2
(commit 2ea23fc57892395d61463951b8ffcba4e4720f03)
geopsyhv 4.0.0 from 3.3.2
(commit 1a2e6cfd3b85039a5be2f6f56e2f908d9d80490c)
geopsynr 1.1.0 from 3.3.2
(commit 8f0be8657f805827cd451123a099f4cd032e8275)
geopsyptmotion 2.0.1 from 3.3.2
(commit e919e37b91b1663739c89f4a241e0212906b8e85)
geopsyrefra 2.0.1 from 3.3.2
(commit 5d83035f77cfe994395e415c09a4c1df5c917741)
geopsytfa 2.1.0 from 3.3.2
(commit 9f0da99aad097a5627720caedf9ec28ff9d4cb80)
hvtfa 2.1.1 from 3.3.2
(commit 4266298abee589a926de3f233c5cf66bd47a8f56)
matfiles 1.0.2 from 3.3.2
(commit 1f654122774812cbb48563ee3ab70ccb93863a10)
</version>
<system>windows 10</system>
</CrashReport>

About the second point you are rigth: both kmin and kmax were set to 0.0 (default values because that occurred at first run of this release). Turning to realistic values, the dispersion graph shows the densities grid.

Thanks a lot and best regards
Luigi

OK I could reproduce the error while processing 3C SPAC.
It was a wrong initialization of the horizontal cross-spectrum matrix.

The processing of 3C SPAC targets is still experimental. You must start spac2disp with option '-VRT'.
If not, you will get an error complaining that there are 3 curves per ring.