Question / warning message not solved : "There is more than one boundary defined"

[rejiba-fay]

Problem description

When I execute a three shots inversion of refraction data (with a reconstructed .sgt file in attached - picking was done using Reflexw ascii export format in traveltime analysis), I have systematically the following warning message:

./core/src/meshentities.cpp:36		GIMLI::Boundary* GIMLI::findBoundary_(const std::set<Boundary*>&)  pls. check, this should not happen.  There is more than one boundary defined.2

Your environment

---

Date: Tue Oct 15 16:11:58 2024 Paris, Madrid (heure d’été)

            OS : Windows (10 10.0.19045 SP0 Multiprocessor Free)
        CPU(s) : 16
       Machine : AMD64
  Architecture : 64bit
           RAM : 63.6 GiB
   Environment : IPython

Python 3.11.9 | packaged by conda-forge | (main, Apr 19 2024, 18:27:10) [MSC
v.1938 64 bit (AMD64)]

       pygimli : 1.5.2
        pgcore : 1.5.0
         numpy : 1.26.4
    matplotlib : 3.9.2
         scipy : 1.14.1
          tqdm : 4.66.5
       IPython : 8.26.0
        tetgen : 0.6.5
       pyvista : 0.44.1

---

Steps to reproduce

import os
import pygimli as pg
import pygimli.physics.traveltime as tt
from pygimli.physics import TravelTimeManager
import matplotlib.pyplot as plt
from pygimli.physics.traveltime.importData import readTOMfile

data=tt.load('test3tirs_lot6002.txt') 

test3tirs_lot6002.txt

tt.show(data)
print(data)

fig, ax = plt.subplots()
lines = tt.drawFirstPicks(ax, data)
mgr = tt.TravelTimeManager(data)
ax, cbar = mgr.showData()
mgr.invert()
#mgr.invert(secNodes=2, paraMaxCellSize=1, paraDepth=-10, zWeight=4, vTop=300, vBottom=3000, verbose=1)
mgr.showFit(firstPicks=True)
ax, cbar = mgr.showResult(coverage=mgr.standardizedCoverage())
rays = mgr.drawRayPaths(ax=ax, color="k", lw=0.3, alpha=0.5)
ax.set_ylabel("depth (m)")
ax.set_xlabel("distance (m)")

Expected behavior

Obviously I obtain almost what is expected (despite some picking issue) ... unless the warning message concerning the boundary

Actual behavior

Data: Sensors: 26 data: 60, nonzero entries: ['err', 'g', 's', 't', 'valid']
15/10/24 - 16:16:16 - pyGIMLi - INFO - Found 1 regions.
15/10/24 - 16:16:16 - pyGIMLi - INFO - Found 1 regions.
15/10/24 - 16:16:16 - pyGIMLi - INFO - Creating forward mesh from region infos.
15/10/24 - 16:16:16 - pyGIMLi - INFO - Creating refined mesh (secnodes: 2) to solve forward task.
15/10/24 - 16:16:16 - pyGIMLi - INFO - Create gradient starting model. 500: 5000
15/10/24 - 16:16:16 - pyGIMLi - INFO - Created startmodel from forward operator: 338, min/max=0.000200/0.002000

./core/src/meshentities.cpp:36		GIMLI::Boundary* GIMLI::findBoundary_(const std::set<Boundary*>&)  pls. check, this should not happen.  There is more than one boundary defined.2
Boundary 0x24bb72a6cb0 rtti: 22 id: 52	N: 4 0  marker: -1 
Boundary 0x24bb72d0240 rtti: 22 id: 0	N: 0 1  marker: 1 

question :

I don't understand the warning message. Question : is it possible to fix it ? Does it come from my sgt file ?

Best regards,
Fayçal