node_map_dof = [0, 2, 8, 6, 1, 5, 7, 3, 4] )
domain = MGridDomain( lengths = ( length, heigth, 0. ),
shape = ( 8, 4, 0 ),
adaptor = mgrid_adaptor )
# Put the tseval (time-stepper) into the spatial context of the
# discretization and specify the response tracers to evaluate there.
#
ts = TS( tse = tseval,
sdomain = domain,
bcond_list = [ BCDof( var = 'u', dof = i, value = 0. ) for i in [domain.get_bottom_left_dofs()[0, 0]] ] +
[ BCDof( var = 'u', dof = i, value = 0. ) for i in [domain.get_bottom_left_dofs()[0, 1]] ] +
[ BCDof( var = 'u', dof = i, value = 0. ) for i in [domain.get_bottom_right_dofs()[0, 1]] ] +
[ BCDof( var = 'u', dof = i, value = -1.2e-4 ) for i in [domain.get_top_middle_dofs()[0, 1]] ],
rtrace_list = [
RTraceGraph( name = 'Fi,right over u_right (iteration)' ,
var_y = 'F_int', idx_y = domain.get_bottom_left_dofs()[0, 0],
var_x = 'U_k', idx_x = domain.get_bottom_left_dofs()[0, 0],
record_on = 'update' ),
# RTraceDomainField(name = 'Flux field' ,
# var = 'eps', idx = 0,
# record_on = 'update'),
RTraceDomainField( name = 'Deformation' ,
var = 'u', idx = 1,
record_on = 'update',
warp = True ),
RTraceDomainField( name = 'Damage' ,
var = 'omega', idx = 0,
record_on = 'update',
