Skip to content

snes solver: Pseudo Transient Continuation method #3181

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 11 commits into
base: next
Choose a base branch
from
+1,302 −627
Open

snes solver: Pseudo Transient Continuation method #3181

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
11 commits
Select commit Hold shift + click to select a range
5ded72d
snes solver: Pseudo Transient Continuation method
bendudson Oct 12, 2025
e9b8d12
snes: Remove pseudo_time switch, use equation_form=pseudo_transient
bendudson Oct 13, 2025
a024d58
snes: Refactoring, split rhs_function from snes_function
bendudson Oct 13, 2025
37187e3
snes: Add pseudo_strategy option
bendudson Oct 13, 2025
07c2882
snes: Replace CHKERRQ with PetscCall macro
bendudson Oct 13, 2025
cc58791
snes: Split PTC implementation into separate functions
bendudson Oct 13, 2025
dcbdc5d
snes: Refactoring Jacobian coloring and pruning
bendudson Oct 14, 2025
48c0b18
snes: Expand manual and minor tidying
bendudson Oct 14, 2025
ca04bc8
snes: Add missing returns, clang tidying
bendudson Oct 14, 2025
9b1eb95
snes: pid_consider_failures and asinh_vars
bendudson Oct 15, 2025
ed642ee
Merge branch 'next' into snes-pseudo-time
bendudson Oct 15, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion include/bout/mesh.hxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -757,7 +757,7 @@ public:
void addRegionPerp(const std::string& region_name, const Region<IndPerp>& region);

/// Converts an Ind2D to an Ind3D using calculation
Ind3D ind2Dto3D(const Ind2D& ind2D, int jz = 0) {
Ind3D ind2Dto3D(const Ind2D& ind2D, int jz = 0) const {
return {ind2D.ind * LocalNz + jz, LocalNy, LocalNz};
}

Expand Down
Loading
Loading