.data


nazwa:  	.asciiz "kolo_na_trojkacie.bmp" #alokuje pamiec
zapytanie_x1:	.asciiz "podaj wspolrzedna x pierwszego wierzcholka \n"
zapytanie_y1:	.asciiz "podaj wspolrzedna y pierwszej wspolrzednej trojkata \n"
zapytanie_x2:	.asciiz "podaj wspolrzedna x drugiej wspolrzednej trojkata \n"
zapytanie_y2:	.asciiz "podaj wspolrzedna y drugiej wspolrzednej trojkata\n"
zapytanie_x3:	.asciiz "podaj wspolrzedna x trzeciej wspolrzednej trojkata \n"
zapytanie_y3:	.asciiz "podaj wspolrzedna y trzeciej wspolrzednej trojkata\n"
#----------------- PLIK BMP -------------
buf_zawartosc:
header:
file_size:	.word 1960054 #800x600*4 +14 +40
reserved:	.word 0
offset_to_bitmap: .word 0
dib_header:
header_size:	.word 40 #
width: 		.word 800
height:		.word 600
planes:		.half  1
bits_per_pixel:	.half 32
compresion:	.word 0
image_size:	.word 1920000 # wielkosc bitmapy bez naglowka
x_meter:	.word 0
y_meter:	.word 0
colors:		.word 0
important:	.word 0


bitmap:		.word 0
#-------------------------------------

tmp:		.ascii "BM"

.text
	#allocate mem (sbrk)
	li	$a0, 1920000
	li	$v0, 9
	syscall
	move	$s7, $v0


	la $a0, zapytanie_x1
	li $v0, 4
	syscall
	la $v0, 5
	syscall
	move $s0, $v0 #zapisuje x pierwszego wierzcholka, s0


	la $a0, zapytanie_y1
	li $v0, 4
	syscall
	la $v0, 5
	syscall
	move $s1, $v0 #zapisuje y pierwszego wierzcholka, s1


	la $a0, zapytanie_x2
	li $v0, 4
	syscall
	li $v0, 5
	syscall
	move $s2, $v0 #zapisuje x drugiego wierzcholka, s3



	la $a0, zapytanie_y2
	li $v0, 4
	syscall
	la $v0, 5
	syscall
	move $s3, $v0 #zapisuje y drugiego wierzcholka, s4

        la $a0, zapytanie_x3
	li $v0, 4
	syscall
	li $v0, 5
	syscall
	move $s4, $v0 #zapisuje x trzeciego wierzcholka, s4

	la $a0, zapytanie_y3
	li $v0, 4
	syscall
	li $v0, 5
	syscall
	move $s5, $v0 #zapisuje y trzeciego wierzcholka, s5



	move $a0, $s0
	move $a1, $s1
	move $a2, $s2
	move $a3, $s3
	jal	plotline


	move $a0, $s2
	move $a1, $s3
	move $a2, $s4
	move $a3, $s5
	jal	plotline

	move $a0, $s4
	move $a1, $s5
	move $a2, $s0
	move $a3, $s1
	jal	plotline


 	jal 	srodek

	 j open_file



plotline: #metoda bresenhama
		li $t9, 0x00ffffff
plotline_x: # a0=x1, a1=y1, a2=x2, a3=y2  t0=sx t1=sy, t3=dx t4=dy t5=error, t6=e2
		bge $a2, $a0, x2_wiekszy	# druga wspolrzedna
		li $t0, -1		#sx=-1
		sub $t3, $a0, $a2 	#dx=x1-x2
		j plotline_y
x2_wiekszy:
		li $t0, 1		#sx=1
		sub $t3, $a2, $a0	#dx=x2-x1
plotline_y:
		bge $a3, $a1, y2_wiekszy	#
		li $t1, -1		#sx=-1
		sub $t4, $a3, $a1	#dx=-(y1-y2)=y2-y1
		j error
y2_wiekszy:
		li $t1, 1
		sub $t4, $a1, $a3	#dy=-(y2-y1)=y1-y2
error:
		add $t5, $t3, $t4 	#error= dx+dy

for:		#(szerokosc*y1+x1)*4=pozycja bitu  wstawianie pixela
set_pixel:	mul $t7, $a1, 800	#szerokosc* y1
		addu $t7, $t7, $a0	#(szerokosc*y1+x1)
		sll $t7, $t7, 2	#przemnozenie przez 4
		addu $t7, $s7, $t7
		sw $t9, ($t7)	# zapis bialgo do t7

x1_equal_x2:	bne $a0, $a2, e2_2_err  #x1!=x2
		beq $a1, $a3, end  #y1=y2


e2_2_err:	sll $t6, $t5, 1 # e2=2*error
		blt  $t6, $t4, e2_blt_dy #  if (e2>=dy) spwdzam czy e2<=dy jesli tak to omijam warunek
		addu $t5, $t5, $t4   	# error+=dy
		addu $a0, $a0, $t0	#x0=x0+sx

e2_blt_dy: # jesli e2<dy  to przechodze do if (e2<dx)
		bgt $t6, $t3, for	# if (e2>dx) to wracam do na poczatek petli for
		addu $t5, $t5, $t3  	#err=err+dx
		addu $a1, $a1, $t1	#y0= y0+sy
		j for

end:
		jr $ra

srodek:

	## wyliczam tu dlugosc odcinka 1-2
		sub $t0, $s0, $s2 #t0= x1-x2
		mul $t0, $t0, $t0 #t0=t0*t0
		sub $t1, $s1, $s3 # t1=y1-y2
		mul $t1, $t1, $t1 # t1=t1*t1
		addu $t0, $t0, $t1 # t0=t0+t1
		mtc1 $t0, $f0
		cvt.s.w $f0, $f0
		sqrt.s $f0, $f0
		cvt.w.s $f0,$f0
		mfc1 $t7, $f0
	## wyliczam tu dlugosc odcinka 1-3
		sub $t0, $s0, $s4 #t0= x1-x2
		mul $t0, $t0, $t0 #t0=t0*t0
		sub $t1, $s1, $s5 # t1=y1-y2
		mul $t1, $t1, $t1 # t1=t1*t1
		addu $t0, $t0, $t1 # t0=t0+t1
		mtc1 $t0, $f0
		cvt.s.w $f0, $f0
		sqrt.s $f0, $f0
		cvt.w.s $f0,$f0
		mfc1 $t8, $f0
	## wyliczam tu dlugosc odcinka 2-3
		sub $t0, $s2, $s4 #t0= x1-x2
		mul $t0, $t0, $t0 #t0=t0*t0
		sub $t1, $s3, $s5 # t1=y1-y2
		mul $t1, $t1, $t1 # t1=t1*t1
		addu $t0, $t0, $t1 # t0=t0+t1
		mtc1 $t0, $f0
		cvt.s.w $f0, $f0
		sqrt.s $f0, $f0
		cvt.w.s $f0,$f0
		mfc1 $t9, $f0

	#### licze srodek okregu wpisanego

		## licze Ox
		add $t0, $t9, $t8
		add $t0, $t0, $t7
		mul $t1, $s0, $t9
		mul $t2, $s2, $t8
		mul $t3, $s4, $t7
		add $t1, $t1, $t2
		add $t1, $t1, $t3
		div $t5, $t1, $t0
		## licze Oy
		add $t0, $t9, $t8
		add $t0, $t0, $t7
		mul $t1, $s1, $t9
		mul $t2, $s3, $t8
		mul $t3, $s5, $t7
		add $t1, $t1, $t2
		add $t1, $t1, $t3
		div $t4, $t1, $t0

		mul $t7, $t4, 800	#szerokosc* y1
		addu $t7, $t7, $t5	#(szerokosc*y1+x1)
		sll $t7, $t7, 2	#przemnozenie przez 4
		addu $t7, $s7, $t7
		sw $t9, ($t7)	# zapis bialgo do t7

		## licze promien korzystajac z tego ze r = 2P/abc  P = 1/2 | [ (x1-x2)(y3-y1) - (y2-y1)(x3-x1) ] |
			sub $t0, $s0, $s2 # x1-x2
			sub $t1, $s5, $s1 #y3-y1
			mul $t0, $t0, $t1  # (x1-x2)(y3-y1)
			sub $t1, $s3, $s1 #y2-y1
			sub $t2, $s4, $s0 #x3-x1
			mul $t1, $t1, $t2 # (y2-y1)(x3-x1)













open_file:
	la $a0, nazwa	  # adres to nazwy pliku
	li $a1, 1	 # 0 odpowiada za czytanie 1 za pisanie
	li $a2, 0	#sposob
	li $v0, 13 	 #odpowiada za otwieranie
	syscall

	move $s6, $v0    #identyfikator pliku uchwyt pliku

write_to_file:
	move $a0, $s6	#pobranie identyfikatora
	la $a1,	tmp	#adress bufora
	li $a2, 2    # liczba bajtow do zapisania
	li $v0, 15
	syscall

	move $a0, $s6
	la $a1,	buf_zawartosc	#adress bufora
	li $a2, 52    # liczba bajtow do zapisania zapisuje wartosc
	li $v0, 15
	syscall

	move $a0, $s6
	move $a1, $s7	#adress bufora
	li $a2, 1920000    # liczba bajtow do zapisania zapisuje wartosc
	li $v0, 15
	syscall

close_file:
	move $a0, $s6  #zamyka file descriptor
	li $v0, 16
	syscall