dolfin team mailing list archive

[Sungick Kim <question142116@xxxxxxxxxxxxxxxxxxxxx>]

Question #142116 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/142116

Sungick Kim gave more information on the question:
https://answers.launchpad.net/dolfin/+question/102501

in the question 102501, I found similar question and it actually solved my problem.
I can write the tensor as the answers in question 102501. FEniCS doesn't complain about it.
But I think that it's not solving problem.
Here are code and the result

<my code>
f00 = Expression("(10.45-9.96)/10.45 *u1c")
f11 = Expression("(5.58-6.05)/5.58 *u1c")
f22 = Expression("(4.86-4.74)/4.86 *u1c")
z = Constant(0.0)
F2= as_matrix( (f00,z,z),(z,f11,z),(z,z,f22) )

I = Identity(v2.cell().d)  # Identity tensor
F1 = I + grad(u2c)
F = F1*F2 # Deformation gradient

<result>

(1) computing Transport problem
  Solving nonlinear variational problem
    Starting Newton solve.
      Applying boundary conditions to linear system.
      Applying boundary conditions to linear system.
      Solving linear system of size 7381 x 7381 (PETSc LU solver, umfpack).
      Applying boundary conditions to linear system.
      Newton iteration 1: r (abs) = 0.000e+00 (tol = 1.000e-10) r (rel) = -nan (tol = 1.000e-14)

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