Квадратный корень на Assemblere

[artemavd]

Имеется алгоритм вычисления квадратного корня на Assemblere. Прошу помочь и пояснить мне, кому не трудно строчки:

	.file	"square.c"
	.section	.rodata
.LC0:
	.string	"\ninsert number: "
.LC1:
	.string	"%d"
	.align 4
.LC2:
	.string	"\n\342\214\212sqrt(%d)\342\214\213 = %d \n\nExiting...\n"
	.text
.globl main
	.type	main, @function
main:
	leal	4(%esp), %ecx
	andl	$-16, %esp
	pushl	-4(%ecx)
	pushl	%ebp
	movl	%esp, %ebp
	pushl	%ecx
	subl	$20, %esp
	subl	$12, %esp
	pushl	$.LC0
	call	puts
	addl	$16, %esp
	subl	$8, %esp
	leal	-20(%ebp), %eax
	pushl	%eax
	pushl	$.LC1
	call	scanf
	addl	$16, %esp
	movl	-20(%ebp), %eax
	movl	%eax, -16(%ebp)
	movl	$1, -12(%ebp)
	movl	$0, -8(%ebp)
	jmp	.L2
.L3:
	movl	-12(%ebp), %eax
	subl	%eax, -16(%ebp)
	incl	-8(%ebp)
	addl	$2, -12(%ebp)
.L2:
	movl	-16(%ebp), %eax
	cmpl	-12(%ebp), %eax
	jge	.L3
	movl	-20(%ebp), %eax
	subl	$4, %esp
	pushl	-8(%ebp)
	pushl	%eax
	pushl	$.LC2
	call	printf
	addl	$16, %esp
	movl	$0, %eax
	movl	-4(%ebp), %ecx
	leave
	leal	-4(%ecx), %esp
	ret
	.size	main, .-main
	.ident	"GCC: (GNU) 4.1.1 20070105 (ALT Linux, build 4.1.1-alt11)"
	.section	.note.GNU-stack,"",@progbits

Помогите плиз только комментариями. Заранее спасибо.

[alexcoder]

Алгоритм описан в педивикии: http://ru.wikipedia.org/wiki/%D0%9A%...B5%D0%BD%D1%8C
Арифметическое извлечение квадратного корня

Для квадратов чисел верны следующие равенства:

1 = 1^2
1 + 3 = 2^2
1 + 3 + 5 = 3^2

и так далее.

То есть, узнать целую часть квадратного корня числа можно, вычитая из него все нечётные числа по порядку, пока остаток не станет меньше следующего вычитаемого числа или равен нулю, и сочтя количество выполненных действий. Например, так:

9 − 1 = 8
8 − 3 = 5
5 − 5 = 0

Выполнено 3 действия, квадратный корень числа 9 равен 3.

Недостатком такого способа является то, что если извлекаемый корень не является целым числом, то можно узнать только его целую часть, но не точнее. В то же время такой способ вполне доступен детям, решающим простейшие математические задачи, требующие извлечения квадратного корня.

	movl	%eax, -16(%ebp);Число из которого нужно извлечь корень
	movl	$1, -12(%ebp);первое нечетное (1)
	movl	$0, -8(%ebp);количество вычитаний
	jmp	.L2
.L3:
	movl	-12(%ebp), %eax;текущее нечетное
	subl	%eax, -16(%ebp):вычесть из остатка текущее нечетное
	incl	-8(%ebp);увеличить количество итераций
	addl	$2, -12(%ebp);получить следующее нечетное(1,3,5...)
.L2:
	movl	-16(%ebp), %eax;остаток числа
	cmpl	-12(%ebp), %eax;сравниваем с текущим нечетным
	jge	.L3 ;если остаток больше или равен текущему нечетному, продолжаем вычисление
	а если меньше, то в [ebp-8] будет находиться значение корня, равное количеству итераций

[Mikl___]

извлечение целочисленного квадратного корня из числа Х, извиняйте, что синтаксис Intel а не AT&T. Считает на много быстрее чем алгоритм предложенный alexcoder

mov ebp,X
        inc ebp
    bsr ecx,ebp
    and ecx,-2
    mov ebx,1
    shl ebx,cl;для уменьшения количества итераций 
        xor eax,eax
@@:    lea ecx,[eax+ebx]
    shr eax,1
    cmp ecx,ebp
    sbb edx,edx
    mov edi,edx
    and edx,ecx
    sub ebp,edx
    and edi,ebx
    or eax,edi
    shr ebx,2
    jnz @b