dolfin team mailing list archive

[osman <question124084@xxxxxxxxxxxxxxxxxxxxx>]

New question #124084 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/124084

Hello,
Just updated my binary install on lucid 64. Tried running Cahn-Hilliard python demo and it fails:
osman@osman-desktop:~/dolfin-doc/demo/pde/cahn-hilliard/python$ python demo.py
Got expression dimension = 2
Starting Newton solve.
Calling FFC just-in-time (JIT) compiler, this may take some time.
Traceback (most recent call last):
  File "demo.py", line 94, in <module>
    solver.solve(problem, u.vector())
  File "demo.py", line 29, in F
    assemble(self.L, tensor=b)
  File "/usr/lib/python2.6/dist-packages/dolfin/fem/assemble.py", line 100, in assemble
    common_cell=common_cell)
  File "/usr/lib/python2.6/dist-packages/dolfin/fem/form.py", line 34, in __init__
    (self._compiled_form, module, self.form_data) = jit(form, form_compiler_parameters, common_cell)
  File "/usr/lib/python2.6/dist-packages/dolfin/compilemodules/jit.py", line 44, in mpi_jit
    return local_jit(*args, **kwargs)
  File "/usr/lib/python2.6/dist-packages/dolfin/compilemodules/jit.py", line 103, in jit
    return jit_compile(form, parameters=p, common_cell=common_cell)
  File "/usr/lib/python2.6/dist-packages/ffc/jitcompiler.py", line 58, in jit
    return jit_form(object, parameters, common_cell)
  File "/usr/lib/python2.6/dist-packages/ffc/jitcompiler.py", line 95, in jit_form
    compile_form(preprocessed_form, prefix=jit_object.signature(), parameters=parameters)
  File "/usr/lib/python2.6/dist-packages/ffc/compiler.py", line 145, in compile_form
    oir = optimize_ir(ir, parameters)
  File "/usr/lib/python2.6/dist-packages/ffc/optimization.py", line 38, in optimize_ir
    oir_integrals = [_optimize_integral_ir(ir) for ir in ir_integrals]
  File "/usr/lib/python2.6/dist-packages/ffc/optimization.py", line 56, in _optimize_integral_ir
    oir = r.optimize_integral_ir(ir)
  File "/usr/lib/python2.6/dist-packages/ffc/tensor/tensoroptimization.py", line 42, in optimize_integral_ir
    ir["AK"][k] = (A0, GK, _optimize_tensor_contraction(A0.A0, rank))
  File "/usr/lib/python2.6/dist-packages/ffc/tensor/tensoroptimization.py", line 80, in _optimize_tensor_contraction
    return optimize(A0)
  File "/usr/lib/python2.6/dist-packages/ferari/binary.py", line 390, in optimize_action
    p = process( Adict )
  File "/usr/lib/python2.6/dist-packages/ferari/binary.py", line 188, in process
    G = get_graph( vecs , rho )
  File "/usr/lib/python2.6/dist-packages/ferari/binary.py", line 72, in get_graph
    (w,l) = dist( vecs[u] , vecs[v] )
  File "/usr/lib/python2.6/dist-packages/ferari/binary.py", line 58, in rho
    c = colinear( u , v )
  File "/usr/lib/python2.6/dist-packages/ferari/binary.py", line 44, in colinear
    uhat = util.unit_vector( u )
  File "/usr/lib/python2.6/dist-packages/ferari/util.py", line 25, in unit_vector
    raise RuntimeError, "divide by zero"
RuntimeError: divide by zero


-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.