luah 1.1.0

Unofficial lua impl by HydroRoll-Team.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44

-- The Computer Language Benchmarks Game
-- http://benchmarksgame.alioth.debian.org/
-- contributed by Mike Pall

local function A(i, j)
    local ij = i+j-1
    return 1.0 / (ij * (ij-1) * 0.5 + i)
  end
  
  local function Av(x, y, N)
    for i=1,N do
      local a = 0
      for j=1,N do a = a + x[j] * A(i, j) end
      y[i] = a
    end
  end
  
  local function Atv(x, y, N)
    for i=1,N do
      local a = 0
      for j=1,N do a = a + x[j] * A(j, i) end
      y[i] = a
    end
  end
  
  local function AtAv(x, y, t, N)
    Av(x, t, N)
    Atv(t, y, N)
  end
  
  local N = 100 -- input argument
  local u, v, t = {}, {}, {}
  for i=1,N do u[i] = 1 end
  
  for i=1,10 do AtAv(u, v, t, N) AtAv(v, u, t, N) end
  
  local vBv, vv = 0, 0
  for i=1,N do
    local ui, vi = u[i], v[i]
    vBv = vBv + ui*vi
    vv = vv + vi*vi
  end
  print(vBv / vv)
  --print(math.sqrt(vBv / vv))