OCD cosmetics.

[dagnn.git] / test-dagnn.lua

diff --git a/test-dagnn.lua b/test-dagnn.lua
index 8f92ccf..b390a29 100755 (executable)
--- a/test-dagnn.lua
+++ b/test-dagnn.lua
@@ -21,16 +21,18 @@
 
 require 'torch'
 require 'nn'
+
+-- require 'cunn'
+
 require 'dagnn'
 
--- torch.setnumthreads(params.nbThreads)
 torch.setdefaulttensortype('torch.DoubleTensor')
-torch.manualSeed(2)
+torch.manualSeed(1)
 
-function checkGrad(model, criterion, input, target)
+function checkGrad(model, criterion, input, target, epsilon)
    local params, gradParams = model:getParameters()
 
-   local epsilon = 1e-5
+   local epsilon = epsilon or 1e-5
 
    local output = model:forward(input)
    local loss = criterion:forward(output, target)
@@ -39,6 +41,8 @@ function checkGrad(model, criterion, input, target)
    model:backward(input, gradOutput)
    local analyticalGradParam = gradParams:clone()
 
+   local err = 0
+
    for i = 1, params:size(1) do
       local x = params[i]
 
@@ -54,23 +58,13 @@ function checkGrad(model, criterion, input, target)
 
       local ana = analyticalGradParam[i]
       local num = (loss1 - loss0) / (2 * epsilon)
-      local err
 
-      if num == ana then
-         err = 0
-      else
-         err = torch.abs(num - ana) / torch.abs(num)
+      if num ~= ana then
+         err = math.max(err, math.abs(num - ana) / math.max(epsilon, math.abs(num)))
       end
-
-      print(
-         'CHECK '
-            .. err
-            .. ' checkGrad ' .. i
-            .. ' analytical ' .. ana
-            .. ' numerical ' .. num
-      )
    end
 
+   return err
 end
 
 function printTensorTable(t)
@@ -84,39 +78,66 @@ function printTensorTable(t)
    end
 end
 
---                     +--> c ----> e --+
---                    /            /     \
---                   /            /       \
---  input --> a --> b ---> d ----+         g --> output
---                          \             /
---                           \           /
---                            +--> f ---+
+--               +-- Linear(10, 10) --> ReLU --> d -->
+--              /                               /
+--             /                               /
+--  --> a --> b -----------> c ---------------+
+--                            \
+--                             \
+--                              +--------------- e -->
+
+dag = nn.DAG()
 
 a = nn.Linear(50, 10)
 b = nn.ReLU()
 c = nn.Linear(10, 15)
-d = nn.Linear(10, 15)
-e = nn.CMulTable()
-f = nn.Linear(15, 15)
-g = nn.CAddTable()
-
-model = nn.DAG()
-
-model:addEdge(a, b)
-model:addEdge(b, nn.Linear(10, 5), nn.ReLU(), nn.Linear(5, 10), c)
-model:addEdge(b, d)
-model:addEdge(c, e)
-model:addEdge(d, e)
-model:addEdge(d, f)
-model:addEdge(e, g)
-model:addEdge(f, nn.Mul(-1), g)
-
-model:setInput(a)
-model:setOutput(g)
+d = nn.CMulTable()
+e = nn.Mul(-1)
+
+dag:connect(a, b, c)
+dag:connect(b, nn.Linear(10, 15), nn.ReLU(), d)
+dag:connect(c, d)
+dag:connect(c, e)
+
+dag:setLabel(a, 'first module')
+
+dag:setInput(a)
+dag:setOutput({ d, e })
+
+-- Check the output of the dot file. Generate a pdf with:
+--
+-- dot ./graph.dot -Lg -T pdf -o ./graph.pdf
+--
+print('Writing ./graph.dot')
+dag:saveDot('./graph.dot')
+
+-- Let's make a model where the dag is inside another nn.Container.
+model = nn.Sequential()
+   :add(nn.Linear(50, 50))
+   :add(dag)
+   :add(nn.CAddTable())
+
+criterion = nn.MSECriterion()
+
+if cunn then
+   print("Using CUDA")
+   model:cuda()
+   criterion:cuda()
+   torch.setdefaulttensortype('torch.CudaTensor')
+   epsilon = 1e-3
+end
 
 local input = torch.Tensor(30, 50):uniform()
 local output = model:updateOutput(input):clone()
-
 output:uniform()
 
-checkGrad(model, nn.MSECriterion(), input, output)
+-- Check that DAG:accGradParameters and friends work okay
+print('Gradient estimate error ' .. checkGrad(model, criterion, input, output, epsilon))
+
+-- Check that we can save and reload the model
+model:clearState()
+torch.save('./test.t7', model)
+local otherModel = torch.load('./test.t7')
+print('Gradient estimate error ' .. checkGrad(otherModel, criterion, input, output, epsilon))
+
+dag:print()