I need to use fenics to do some processing on my solution field. I solve my PDE with Abaqus and extract the displacement for all the nodes and I can put it in vector D. I am using the same mesh for both Abaqus and Fenics. so is there any way that I can assign the vector of my solution which is an array of number to U so I can have the solution function ?
from __future__ import division
from dolfin import *
import numpy
mesh=UnitSquareMesh(10, 10)
#D is the displacment obtained from Abaqus
V = VectorFunctionSpace(mesh, "Lagrange", 1)
U = Function(V)
Thanks.