Skip to content

GEOPY-440 #40

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

Merged
domfournier merged 3 commits into from
Apr 9, 2024
Merged

GEOPY-440 #40

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
107aab0
Add flag to Spherical directive to preserve reference angle
domfournier Apr 9, 2024
6595bf6
Change reference_angles to list
domfournier Apr 9, 2024
ab36692
Merge branch 'release/v0.19.0.dev8+geoapps.0.12.0' into GEOPY-440
domfournier Apr 9, 2024
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
9 changes: 7 additions & 2 deletions SimPEG/directives/directives.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3162,6 +3162,7 @@ class VectorInversion(InversionDirective):
def __init__(
self, simulations: list, regularizations: ComboObjectiveFunction, **kwargs
):
self.reference_angles = (False, False, False)
self.simulations = simulations
self.regularizations = regularizations

Expand Down Expand Up @@ -3216,12 +3217,16 @@ def endIter(self):
model[indices] = vec_model.flatten()

angle_map = []
for ind, reg_fun in enumerate(self.regularizations.objfcts):
for ind, (reg_fun, ref_angle) in enumerate(
zip(self.regularizations.objfcts, self.reference_angles)
):
reg_fun.model = model
reg_fun.reference_model[indices] = vec_ref

if ind > 0:
reg_fun.alpha_s = 0
if not ref_angle:
reg_fun.alpha_s = 0

reg_fun.eps_q = np.pi
reg_fun.units = "radian"
angle_map.append(reg_fun.mapping)
Expand Down