kaufmanno / GSDMA

GNU General Public License v3.0
0 stars 1 forks source link

Improve the "Example of reading borehole data from electrodes string files inserted into a sqlite database" notebook for documentation #14

Open kaufmanno opened 3 years ago

kaufmanno commented 3 years ago

I did it.

kaufmanno commented 3 years ago

No it is not fixed yet

YanickNdepo commented 3 years ago

Here the description of the task in the documentation is not clear. It’s written ‘building of a lexicon’ but the code don’t really show that. For enhancement, try to describe every task in documentation by specifying the right code used

Joseph-Wandji commented 3 years ago

This code cannot work to my notebook. Perhaps, the repertory is not really specify for the main code. Code work but we need to classified this ones and their main functions in the same folder.

kaufmanno commented 3 years ago

This code cannot work to my notebook. Perhaps, the repertory is not really specify for the main code. Code work but we need to classified this ones and their main functions in the same folder.

@Joseph-Wandji You should create a new issue and report the problem in details.

Here is an example of how to submit a such an issue report:

Problem description

I am trying to invert an ert dataset with geostatistical regularization. I keep getting an error message and I can not figure out where it comes from. I am actually not sure that I am using the ERTManager in the right way. Any help is appreciated ;-)

Olivier

Your environment

Operating system: Ubuntu 18.04LTS Python version: 3.7.6 pyGIMLi version: 1.1.0+0.g47c8d810 (with local changes) Way of installation: Conda package

Steps to reproduce

Here is a short script that produces the encountered error.

import pygimli as pg

data=pg.DataContainerERT('Eugies_ERT.dat') ert = pg.physics.ERTManager(sr=False, useBert=True, verbose=True, debug=False) ert.applyMesh(ert.createMesh(data=data, paraDX=0.5, maxCellArea=100.0), ignoreRegionManager=False) C = pg.matrix.GeostatisticConstraintsMatrix(mesh=ert.fop.mesh(), I=[25, 5], dip=-25) ert.fop.setConstraints(C) ert.invert(data=data, lam=10) ert.showResultAndFit(cMap='jet')

Expected behavior

I was expecting that ert.invert would run without returning this error message...

Actual behavior

an error message is thrown by ert.invert :


RuntimeError Traceback (most recent call last)

in ----> 1 ert.invert(data=data, lam=10) 2 ert.showResultAndFit(cMap='jet') ~/miniconda3/envs/pg/lib/python3.7/site-packages/pygimli/frameworks/methodManager.py in invert(self, data, mesh, zWeight, startModel, **kwargs) 695 errorVals = self._ensureError(self.fop.data, dataVals) 696 --> 697 if self.fop.mesh() is None: 698 pg.critical('Please provide a mesh') 699 # inversion will call this itsself as default behaviour ~/miniconda3/envs/pg/lib/python3.7/site-packages/pygimli/frameworks/modelling.py in mesh(self) 520 """""" 521 if self._regionManagerInUse and self._regionChanged: --> 522 self.createFwdMesh_() 523 524 return super(Modelling, self).mesh() ~/miniconda3/envs/pg/lib/python3.7/site-packages/pygimli/frameworks/modelling.py in createFwdMesh_(self) 515 self.setMeshPost(m) 516 self._regionChanged = False --> 517 super(Modelling, self).setMesh(m, ignoreRegionManager=True) 518 519 def mesh(self): RuntimeError: /home/wagner/miniconda3/conda-bld/pygimli_1589132014332/work/gimli/gimli/core/src/matrix.h: 181 virtual void GIMLI::MatrixBase::clear() not yet implemented libgimli-v1.1.0
Joseph-Wandji commented 3 years ago

@YanickNdepo has resolved this problem. Pleasure to interact with him

kaufmanno commented 3 years ago

Changed the notebook in commit 596367a9b5da3188ff2fd98a66d37004d4b42f17 to improve some of the comments. An .rst file was not created yet for update of the documentation...