Skip to content

fix: Fix bug in L2 norm calculation for isothermal#3954

Open
dkachuma wants to merge 1 commit intodevelopfrom
dkachuma/fix-ls-residual-check
Open

fix: Fix bug in L2 norm calculation for isothermal#3954
dkachuma wants to merge 1 commit intodevelopfrom
dkachuma/fix-ls-residual-check

Conversation

@dkachuma
Copy link
Contributor

For isothermal models, the normaliser for the L2 normal calculation was not being populated leading to a divide by zero down the line. This PR simply sets this to the minimum value to bypass the floating point exception.

@dkachuma dkachuma self-assigned this Jan 30, 2026
@dkachuma dkachuma requested a review from rrsettgast as a code owner January 30, 2026 22:07
@dkachuma dkachuma added the type: bug Something isn't working label Jan 30, 2026
@dkachuma dkachuma added ci: run CUDA builds Allows to triggers (costly) CUDA jobs flag: no rebaseline Does not require rebaseline ci: run code coverage enables running of the code coverage CI jobs labels Jan 30, 2026
@codecov
Copy link

codecov bot commented Jan 30, 2026

Codecov Report

✅ All modified and coverable lines are covered by tests.
✅ Project coverage is 58.87%. Comparing base (2acb3ef) to head (c11a5d7).
⚠️ Report is 1 commits behind head on develop.

Additional details and impacted files
@@           Coverage Diff            @@
##           develop    #3954   +/-   ##
========================================
  Coverage    58.87%   58.87%           
========================================
  Files         1353     1353           
  Lines       112953   112954    +1     
========================================
+ Hits         66504    66505    +1     
  Misses       46449    46449           

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:
  • ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should defaults in well also be set ?

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Well looks to be handled differently. There's no volume balance and it seems it's all sized correctly (1 for isothermal and 2 for thermal). CompositionalMultiphaseWell.cpp:1379-1388

@dkachuma dkachuma added flag: ready for review ci: run integrated tests Allows to run the integrated tests in GEOS CI labels Feb 3, 2026
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

ci: run code coverage enables running of the code coverage CI jobs ci: run CUDA builds Allows to triggers (costly) CUDA jobs ci: run integrated tests Allows to run the integrated tests in GEOS CI flag: no rebaseline Does not require rebaseline flag: ready for review type: bug Something isn't working

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants