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

Eltwise unary reciprocal and rsqrt require Newton iteration on ASIMD #774

Open
hfp opened this issue Jun 5, 2023 · 2 comments
Open

Eltwise unary reciprocal and rsqrt require Newton iteration on ASIMD #774

hfp opened this issue Jun 5, 2023 · 2 comments
Labels

Comments

@hfp
Copy link
Collaborator

hfp commented Jun 5, 2023

Eltwise unary reciprocal and rsqrt require Newton iteration to meet reasonable accuracy with ASIMD. Accuracy is reasonable for the SVE based code-path (Newton iteration).

@hfp hfp changed the title Eltwise unary reciprocal and rsqrt require Newton iteration to meet accuracy with ASIMD Eltwise unary reciprocal and rsqrt require Newton iteration on ASIMD Jun 5, 2023
@hfp
Copy link
Collaborator Author

hfp commented Jun 12, 2023

This is a low-priority issue if an issue at all. The unwritten policy for numerical accuracy can be "architectural accuracy" being "embraced" by default (even if exposing big differences betwee uArchs), i.e., LIBXSMM_ULP_PRECISION may be the only way to trigger higher precision (or Newton iterations on top of "native" precision).

Remaining work (need to keep issue open?): check if LIBXSMM_ULP_PRECISION can be made a vehicle to close the gap between architectural differences. For example, on x86 native sqrt, native rcp, and native rsqrt exist depending on type (FP32 vs F64), but even that exposes differences due to type (and/or vendor's implementation). On ARM64, precision and ISA coverage is apparently the weakest.

@hfp hfp added the revisit label Jun 15, 2023
@hfp hfp closed this as completed Jun 15, 2023
@hfp hfp reopened this Jun 20, 2023
@hfp
Copy link
Collaborator Author

hfp commented Aug 17, 2023

To confirm: SVE performs Newton iteration(s) as opposed to ASIMD.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant