как найти корень из числа

[$Дмитрий]

как найти корень в n-ой степени из числа

[Carbon]

Power(a,1/n)

[Квэнди]

результат:=EXP(1/<степень>*Ln(числа));

[Carbon]

Лучше юзать Power, т.к. он не вылетает при основании <=0.

[Квэнди]

Да? сравните использование памяти при использовании power.

[Carbon]

Если критична память нужно юзать asm.

[Квэнди]

function Power(const Base, Exponent: Extended): Extended;
const
Max : Double = MaxInt;
var
IntExp : Integer;
asm
fld Exponent
fld st {copy to st(1)}
fabs {abs(exp)}
fld Max
fcompp {leave exp in st(0)}
fstsw ax
sahf
jb @@RealPower {exp > MaxInt}
fld st {exp in st(0) and st(1)}
frndint {round(exp)}
fcomp {compare exp and round(exp)}
fstsw ax
sahf
jne @@RealPower
fistp IntExp
mov eax, IntExp {eax=Trunc(Exponent)}
mov ecx, eax
cdq
fld1 {Result=1}
xor eax, edx
sub eax, edx {abs(exp)}
jz @@Exit
fld Base
jmp @@Entry
@@Loop:
fmul st, st {Base * Base}
@@Entry:
shr eax, 1
jnc @@Loop
fmul st(1), st {Result * X}
jnz @@Loop
fstp st
cmp ecx, 0
jge @@Exit
fld1
fdivrp {1/Result}
jmp @@Exit
@@RealPower:
fld Base
ftst
fstsw ax
sahf
jz @@Done
fldln2
fxch
fyl2x
fxch
fmulp st(1), st
fldl2e
fmulp st(1), st
fld st(0)
frndint
fsub st(1), st
fxch st(1)
f2xm1
fld1
faddp st(1), st
fscale
@@Done:
fstp st(1)
@@Exit:
end;



function Exp(const X: Extended): Extended;
asm
{ e**x = 2**(x*log2(e)) }
FLD X
FLDL2E { y := x*log2e; }
FMUL
FLD ST(0) { i := round(y); }
FRNDINT
FSUB ST(1), ST { f := y - i; }
FXCH ST(1) { z := 2**f }
F2XM1
FLD1
FADD
FSCALE { result := z * 2**i }
FSTP ST(1)
end;

function Ln(const X: Extended): Extended;
asm
FLD X
FLDLN2
FXCH
FYL2X
FWAIT
end;


Сравнивайте. Плюс: использование Power увеличит ваш дистрибутив на размер модуля Math не забывайте об этом.

* и деление не стал приводить к атомарности, по причине того что и то и другое в исходных кодах функций используется

[zetrix]

Тема закрыта, во избежании массового наказания за флуд!
Протесты в личку!