numeric

C++ library with numerical algorithms
git clone git://src.adamsgaard.dk/numeric # fast
git clone https://src.adamsgaard.dk/numeric.git # slow
Log | Files | Refs | README | LICENSE Back to index

lua-linarr.lua (341B)

      1 #!/usr/bin/env lua
      2 
      3 N = tonumber(arg[1])
      4 if N == nil then error("Sorry, I need matrix width as command line argument") end
      5 
      6 A = {}
      7 B = {}
      8 for i=1,(N*N) do
      9   A[i] = 2.0
     10   B[i] = i-1
     11 end
     12 
     13 C = {}
     14 for i=1,N do
     15   for j=1,N do
     16     sum = 0.0
     17     for k=1,N do
     18       sum = sum + A[(k-1)*N+i] * B[(j-1)*N+k]
     19     end
     20     C[(i-1)*N + j] = sum
     21   end
     22 end