MD5/Implementation: Difference between revisions

Line 4,255:

=={{header|Phix}}==

Non-optimised. Originally written by Davi Tassinari de Figueiredo.

<lang Phix>-- ~~demo\rosetta\md5.exw~~

-- demo\rosetta\md5.exw

function uxor(atom data1,atom data2)

with javascript_semantics

atom result = xor_bits(data1,data2)

function uxor(atom data1,atom data2)

if result<0 then result += #100000000 end if

atom result = xor_bits(data1,data2)

return result

if result<0 then result += #100000000 end if

end function

return result

end function

function uor(atom data1,atom data2)

atom result = or_bits(data1,data2)

function uor(atom data1,atom data2)

if result<0 then result += #100000000 end if

atom result = or_bits(data1,data2)

return result

if result<0 then result += #100000000 end if

end function

return result

end function

function r32(atom a)

return remainder(a,#100000000)

function r32(atom a)

end function

return remainder(a,#100000000)

end function

function rol(atom word,integer bits)

-- left rotate the bits of a 32-bit number by the specified number of bits

function rol(atom word,integer bits)

return r32(word*power(2,bits))+floor(word/power(2,32-bits))

-- left rotate the bits of a 32-bit number by the specified number of bits

end function

return r32(word*power(2,bits))+floor(word/power(2,32-bits))

end function

constant K =

{#d76aa478, #e8c7b756, #242070db, #c1bdceee, #f57c0faf, #4787c62a, #a8304613, #fd469501,

constant K =

#698098d8, #8b44f7af, #ffff5bb1, #895cd7be, #6b901122, #fd987193, #a679438e, #49b40821,

{#d76aa478, #e8c7b756, #242070db, #c1bdceee, #f57c0faf, #4787c62a, #a8304613, #fd469501,

#f61e2562, #c040b340, #265e5a51, #e9b6c7aa, #d62f105d, #02441453, #d8a1e681, #e7d3fbc8,

#698098d8, #8b44f7af, #ffff5bb1, #895cd7be, #6b901122, #fd987193, #a679438e, #49b40821,

#21e1cde6, #c33707d6, #f4d50d87, #455a14ed, #a9e3e905, #fcefa3f8, #676f02d9, #8d2a4c8a,

#f61e2562, #c040b340, #265e5a51, #e9b6c7aa, #d62f105d, #02441453, #d8a1e681, #e7d3fbc8,

#fffa3942, #8771f681, #6d9d6122, #fde5380c, #a4beea44, #4bdecfa9, #f6bb4b60, #bebfbc70,

#21e1cde6, #c33707d6, #f4d50d87, #455a14ed, #a9e3e905, #fcefa3f8, #676f02d9, #8d2a4c8a,

#289b7ec6, #eaa127fa, #d4ef3085, #04881d05, #d9d4d039, #e6db99e5, #1fa27cf8, #c4ac5665,

#fffa3942, #8771f681, #6d9d6122, #fde5380c, #a4beea44, #4bdecfa9, #f6bb4b60, #bebfbc70,

#f4292244, #432aff97, #ab9423a7, #fc93a039, #655b59c3, #8f0ccc92, #ffeff47d, #85845dd1,

#289b7ec6, #eaa127fa, #d4ef3085, #04881d05, #d9d4d039, #e6db99e5, #1fa27cf8, #c4ac5665,

#6fa87e4f, #fe2ce6e0, #a3014314, #4e0811a1, #f7537e82, #bd3af235, #2ad7d2bb, #eb86d391}

#f4292244, #432aff97, #ab9423a7, #fc93a039, #655b59c3, #8f0ccc92, #ffeff47d, #85845dd1,

#6fa87e4f, #fe2ce6e0, #a3014314, #4e0811a1, #f7537e82, #bd3af235, #2ad7d2bb, #eb86d391}

constant m_block = {1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,

2, 7,12, 1, 6,11,16, 5,10,15, 4, 9,14, 3, 8,13,

constant m_block = {1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,

6, 9,12,15, 2, 5, 8,11,14, 1, 4, 7,10,13,16, 3,

2, 7,12, 1, 6,11,16, 5,10,15, 4, 9,14, 3, 8,13,

1, 8,15, 6,13, 4,11, 2, 9,16, 7,14, 5,12, 3,10}

6, 9,12,15, 2, 5, 8,11,14, 1, 4, 7,10,13,16, 3,

1, 8,15, 6,13, 4,11, 2, 9,16, 7,14, 5,12, 3,10}

constant c_words = {#67452301,#efcdab89,#98badcfe,#10325476}

constant c_words = {#67452301,#efcdab89,#98badcfe,#10325476}

sequence words

sequence words

function divide_in_words(sequence message)

-- Divides the string into words (32-bit numbers)

function divide_in_words(sequence message)

sequence res

-- Divides the string into words (32-bit numbers)

res = repeat(0,length(message)/4)

integer l = floor(length(message)/4)

for word=1 to length(message)/4 do

sequence res = repeat(0,l)

res[word] = bytes_to_int(message[word*4-3..word*4])

for word=1 to l do

end for

res[word] = bytes_to_int(message[word*4-3..word*4])

return res

end for

end function

return res

end function

procedure process_block(sequence block)

-- Updates the words according to the contents of the block

procedure process_block(sequence block)

atom a,b,c,d

-- Updates the words according to the contents of the block

atom a,b,c,d

block = divide_in_words(block)

block = divide_in_words(block)

a = words[1]

b = words[2]

a = words[1]

c = words[3]

b = words[2]

d = words[4]

c = words[3]

d = words[4]

-- Round 1

for step=1 to 16 by 4 do

-- Round 1

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uor(and_bits(b,c),and_bits(not_bits(b),d))), 7))

for step=1 to 16 by 4 do

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uor(and_bits(a,b),and_bits(not_bits(a),c))),12))

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uor(and_bits(b,c),and_bits(not_bits(b),d))), 7))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uor(and_bits(d,a),and_bits(not_bits(d),b))),17))

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uor(and_bits(a,b),and_bits(not_bits(a),c))),12))

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uor(and_bits(c,d),and_bits(not_bits(c),a))),22))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uor(and_bits(d,a),and_bits(not_bits(d),b))),17))

end for

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uor(and_bits(c,d),and_bits(not_bits(c),a))),22))

end for

-- Round 2

for step=17 to 32 by 4 do

-- Round 2

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uor(and_bits(b,d),and_bits(c,not_bits(d)))), 5))

for step=17 to 32 by 4 do

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uor(and_bits(a,c),and_bits(b,not_bits(c)))), 9))

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uor(and_bits(b,d),and_bits(c,not_bits(d)))), 5))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uor(and_bits(d,b),and_bits(a,not_bits(b)))),14))

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uor(and_bits(a,c),and_bits(b,not_bits(c)))), 9))

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uor(and_bits(c,a),and_bits(d,not_bits(a)))),20))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uor(and_bits(d,b),and_bits(a,not_bits(b)))),14))

end for

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uor(and_bits(c,a),and_bits(d,not_bits(a)))),20))

end for

-- Round 3

for step=33 to 48 by 4 do

-- Round 3

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uxor(b,xor_bits(c,d))), 4))

for step=33 to 48 by 4 do

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uxor(a,xor_bits(b,c))),11))

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uxor(b,xor_bits(c,d))), 4))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uxor(d,xor_bits(a,b))),16))

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uxor(a,xor_bits(b,c))),11))

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uxor(c,xor_bits(d,a))),23))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uxor(d,xor_bits(a,b))),16))

end for

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uxor(c,xor_bits(d,a))),23))

end for

-- Round 4

for step=49 to 64 by 4 do

-- Round 4

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uxor(c,or_bits(b,not_bits(d)))), 6))

for step=49 to 64 by 4 do

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uxor(b,or_bits(a,not_bits(c)))),10))

a = r32(b+rol(r32(a+block[m_block[step ]]+K[step ]+uxor(c,or_bits(b,not_bits(d)))), 6))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uxor(a,or_bits(d,not_bits(b)))),15))

d = r32(a+rol(r32(d+block[m_block[step+1]]+K[step+1]+uxor(b,or_bits(a,not_bits(c)))),10))

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uxor(d,or_bits(c,not_bits(a)))),21))

c = r32(d+rol(r32(c+block[m_block[step+2]]+K[step+2]+uxor(a,or_bits(d,not_bits(b)))),15))

end for

b = r32(c+rol(r32(b+block[m_block[step+3]]+K[step+3]+uxor(d,or_bits(c,not_bits(a)))),21))

end for

-- Update the words

words[1] = r32(words[1]+a)

-- Update the words

words[2] = r32(words[2]+b)

words = deep_copy(words)

words[3] = r32(words[3]+c)

words[1] = r32(words[1]+a)

words[4] = r32(words[4]+d)

words[2] = r32(words[2]+b)

end procedure

words[3] = r32(words[3]+c)

words[4] = r32(words[4]+d)

function pad_message(sequence message)

end procedure

-- Add bytes to the end of the message so it can be divided

-- in an exact number of 64-byte blocks.

function pad_message(sequence message)

integer bytes_to_add

-- Add bytes to the end of the message so it can be divided

bytes_to_add = 64-remainder(length(message)+9,64)

-- in an exact number of 64-byte blocks.

if bytes_to_add=64 then bytes_to_add = 0 end if

integer l = length(message),

bytes_to_add = 64-remainder(l+9,64)

message = message&#80&repeat(0,bytes_to_add)&

if bytes_to_add=64 then bytes_to_add = 0 end if

int_to_bytes(length(message)*8)&{0,0,0,0}

message = deep_copy(message)

return message

message &= #80

end function

message &= repeat(0,bytes_to_add)

message &= int_to_bytes(l*8)&{0,0,0,0}

global function md5(sequence message)

return message

-- Given a string, returns a 16-byte hash of it.

end function

words = c_words -- Initialize the H words

global function md5(sequence message)

message = pad_message(message) -- Add bytes to the message

-- Given a string, returns a 16-byte hash of it.

-- Process each 64-byte block

words = c_words -- Initialize the H words

for block=1 to length(message) by 64 do

process_block(message[block..block+63])

message = pad_message(message) -- Add bytes to the message

end for

-- ~~Convert~~ ~~hash~~ ~~into~~ ~~bytes~~

-- Process each 64-byte block

for block=1 to length(message) by 64 do

return int_to_bytes(words[1])& -- Return the hash

process_block(message[block..block+63])

int_to_bytes(words[2])&

end for

int_to_bytes(words[3])&

int_to_bytes(words[4])

-- Convert hash into bytes

return int_to_bytes(words[1])& -- Return the hash

end function

int_to_bytes(words[2])&

int_to_bytes(words[3])&

if include_file()=1 then

int_to_bytes(words[4])

string fmt = "0x%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X\n"

end function

printf(1,fmt,md5(""))

printf(1,fmt,md5("a"))

if platform()=JS or include_file()=1 then

printf(1,fmt,md5("abc"))

string fmt = "0x%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X%02X\n"

printf(1,fmt,md5("message digest"))

printf(1,fmt,md5("abcdefghijklmnopqrstuvwxyz"))

printf(1,fmt,md5(""))

printf(1,fmt,md5("ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789"))

printf(1,fmt,md5("a"))

printf(1,fmt,md5("12345678901234567890123456789012345678901234567890123456789012345678901234567890"))

printf(1,fmt,md5("abc"))

end if</lang>

printf(1,fmt,md5("message digest"))

printf(1,fmt,md5("abcdefghijklmnopqrstuvwxyz"))

printf(1,fmt,md5("ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789"))

printf(1,fmt,md5("12345678901234567890123456789012345678901234567890123456789012345678901234567890"))

end if

<pre>