Skip to content
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.

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

Update to use safe scaling algorithm from Reference-LAPACK PR 527 #4143

Merged
merged 8 commits into from
Sep 1, 2023

Conversation

martin-frbg
Copy link
Collaborator

fixes original topic of #4130

interface/rotg.c Outdated
*C = ONE;
*S = ZERO;
*DA = ZERO;
DA = ZERO;
Copy link
Contributor

Choose a reason for hiding this comment

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

I think this is unintentional.

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

Certainly hope so...but this is just a dump of my first draft, too much work too little brain right now

interface/rotg.c Outdated
@@ -112,7 +112,7 @@ void CNAME(FLOAT *DA, FLOAT *DB, FLOAT *C, FLOAT *S){
if (adb == ZERO) {
*C = ONE;
*S = ZERO;
DA = ZERO;
*DA = ZERO;
Copy link
Contributor

Choose a reason for hiding this comment

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

I'm not sure about this change. If the input vector is [a,0], then the rotation is the identity and the rotated vector [a,0]. Maybe your initial intention was to delete the whole line as in the other #if branch?

Copy link
Collaborator Author

Choose a reason for hiding this comment

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

yeah, looks like a cut-n-paste gone wrong. thanks for your vigilance - still haven't run any tests with this

@martin-frbg martin-frbg marked this pull request as draft July 17, 2023 21:33
@martin-frbg martin-frbg marked this pull request as ready for review September 1, 2023 07:05
@martin-frbg martin-frbg merged commit 33797c4 into OpenMathLib:develop Sep 1, 2023
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.

2 participants