Skip to content

Fix -J by correctly using orientation object and fixing read/write vec - #90

Open
MiCurry wants to merge 1 commit into
magnetotellurics:developfrom
MiCurry:fix-J
Open

MiCurry wants to merge 1 commit into
magnetotellurics:developfrom
MiCurry:fix-J

Conversation

@MiCurry

@MiCurry MiCurry commented Sep 22, 2026

Copy link
Copy Markdown
Member

This commit fixes an error in -J where NaN's would be found in b. This was because the orientation (orient) was being used incorrectly and was uninitialized. This caused NaN's to form in b, which was then detected later by the solver.

Furthermore, the sensitvity matrix was not correctly being written out as the modelParam readvec/writevec were left empty. This commit makes those routine call the binary versions when called. Later, these can be updated to write out different file types.

This issue was reported by Xiaolei (@xltu) in issue #87 and reported the fix.

This commit fixes an error in -J where NaN's would be found in b. This was
because the orientation (orient) was being used incorrectly and was
uninitialized. This caused NaN's to form in b, which was then detected later by
the solver.

Furthermore, the sensitvity matrix was not correctly being written out as the
modelParam readvec/writevec were left empty. This commit makes those routine
call the binary versions when called. Later, these can be updated to write out
different file types.

This issue was reported by Xiaolei (@xltu) in issue magnetotellurics#87 and reported the fix.
@MiCurry MiCurry changed the title Fix -J by correctly using orient and fixing read/write vec Fix -J by correctly using orientation object and fixing read/write vec Sep 22, 2026

This branch has not been deployed

No deployments
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant