Borrar filtros
Borrar filtros

Is rsqrt the same as Fast inverse square root?

26 visualizaciones (últimos 30 días)
Hongbo Zhu
Hongbo Zhu el 13 de Sept. de 2021
Editada: Walter Roberson el 31 de Jul. de 2022
Hello,
I would like to us sqrt() and run it on a embeded controller. I was wondering sqrt() is the same as Fast inverse square root ( https://en.wikipedia.org/wiki/Fast_inverse_square_root ). if not, is it efficent as Fast inverse square root. Does it support code generation? Can I use it in embeded control software?
Thanks,
Hongbo
  4 comentarios
Rik
Rik el 13 de Sept. de 2021
It might be, although it doesn't have to be. Implementing this C code as a mex function might speed up what you want to do, or might slow it down.
If you don't know what you're working with well enough to know this code is not Matlab but C, I am not sure you would be able to write something faster than the engineers working for Mathworks. This algorithm is famous, I would expect them to know it. It would highly surprise me if they haven't tried implementing it to see the relative performance.
What is the reason for your concern about performance? If performance is so important, why aren't you doing everything in C anyway (or in a language even closer to your bare metal of choice)? I personally have the view that Matlab (and Simulink) are great tools for rapid prototyping and doing complex things. However, they tend to be the slower choice if performance is absolutely critical (it often isn't).
Hongbo Zhu
Hongbo Zhu el 13 de Sept. de 2021
@Rik Thank you very much for your quick response and explanation. It will be great if someone from MathWork can confirm rsqrt() is efficient for embeded control system.

Iniciar sesión para comentar.

Respuesta aceptada

stozaki
stozaki el 13 de Sept. de 2021
Hi,
The Simulink sqrt block supports rsqrt.
The algorithm is Newton-Raphson, so it can be considered the same as the one described on the wiki.
With that setting, C code generation is possible.

Más respuestas (1)

Jan
Jan el 13 de Sept. de 2021
i = * ( long * ) &y
This is equivalent to:
y = single(pi);
i = typecast(y, 'int32');
The shown code of Q_rsqrt is an approximation of 1/sqrt(x) for single precision floating point values. It might have a fair speed on CPUs without a floating point unit, but it is less accurate then the direct calculation.
  1 comentario
Hongbo Zhu
Hongbo Zhu el 13 de Sept. de 2021
Hi @Jan, Thank you very much. I noticed there is also a simulink block named Float Typecast. Does it supported by code generation for embeded controller?

Iniciar sesión para comentar.

Categorías

Más información sobre Simulink en Help Center y File Exchange.

Productos


Versión

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by