In [2]:

    

from dolfin import *
import numpy as np
import scipy.sparse as sp
import numpy
# from scipy2Trilinos import scipy_csr_matrix2CrsMatrix
# from PyTrilinos import Epetra, ML, AztecOO

nn = 4
mesh = RectangleMesh(0, 0, 1, 1, nn, nn,'crossed')

order = 2
Magnetic = FunctionSpace(mesh, "N1curl", order)
Lagrange = FunctionSpace(mesh, "CG", order)
DG = FunctionSpace(mesh, "DG", order-1)


    

    




DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.
DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.
DEBUG:UFL:No integrals left after transformation, returning empty form.
DEBUG:FFC:Reusing form from cache.

In [ ]:

    

mesh.


    

In [34]:

    

u = Function(Lagrange)
DGu = TrialFunction(DG)
uu = grad(u)


    

In [25]:

    

N = FacetNormal(mesh)
T = as_tensor([N[1],-N[0]])


    

In [35]:

    

assemble(inner(uu,T)*DGu*dS)


    

    




Form argument must be restricted.






    




ERROR:UFL:Form argument must be restricted.






    




---------------------------------------------------------------------------
UFLException                              Traceback (most recent call last)
<ipython-input-35-8dbbadb45e57> in <module>()
----> 1 assemble(inner(uu,T)*DGu*dS)

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/dolfin/fem/assembling.pyc in assemble(form, tensor, mesh, coefficients, function_spaces, cell_domains, exterior_facet_domains, interior_facet_domains, reset_sparsity, add_values, finalize_tensor, keep_diagonal, backend, form_compiler_parameters, bcs)
    166                            coefficients=coefficients,
    167                            form_compiler_parameters=form_compiler_parameters,
--> 168                            common_cell=common_cell)
    169 
    170     # Set mesh if specified (important for functionals without a function spaces)

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/dolfin/fem/form.pyc in __init__(self, form, function_spaces, coefficients, subdomains, form_compiler_parameters, common_cell)
     54             self._compiled_form, module, self.form_data, prefix \
     55                                  = jit(form, form_compiler_parameters,
---> 56                                        common_cell)
     57 
     58         elif isinstance(form, ufc.form):

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/dolfin/compilemodules/jit.pyc in mpi_jit(*args, **kwargs)
     58         # Just call JIT compiler when running in serial
     59         if MPI.num_processes() == 1:
---> 60             return local_jit(*args, **kwargs)
     61 
     62         # Compile first on process 0

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/dolfin/compilemodules/jit.pyc in jit(form, form_compiler_parameters, common_cell)
    120         raise RuntimeError, "Form compiler must implement the jit function."
    121 
--> 122     return jit_compile(form, parameters=p, common_cell=common_cell)

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ffc/jitcompiler.pyc in jit(ufl_object, parameters, common_cell)
     76         return jit_element(ufl_object, parameters)
     77     else:
---> 78         return jit_form(ufl_object, parameters, common_cell)
     79 
     80 def _auto_select_degree(elements):

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ffc/jitcompiler.pyc in jit_form(form, parameters, common_cell)
    167     # Compute form metadata and extract preprocessed form
    168     form_data = form.compute_form_data(common_cell=common_cell,
--> 169                                        element_mapping=element_mapping)
    170 
    171     # Wrap input

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/form.pyc in compute_form_data(self, object_names, common_cell, element_mapping)
    122                                          object_names=object_names,
    123                                          common_cell=common_cell,
--> 124                                          element_mapping=element_mapping)
    125         # Always validate arguments, keeping sure that the validation works
    126         self._form_data.validate(object_names=object_names,

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/preprocess.pyc in preprocess(form, object_names, common_cell, element_mapping)
    116 
    117     # Propagate restrictions of interior facet integrals to the terminal nodes
--> 118     form = propagate_restrictions(form) # INTEGRAL, EXPR
    119 
    120     # --- BEGIN DOMAIN SPLITTING AND JOINING

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/propagate_restrictions.pyc in propagate_restrictions(expression)
    109 def propagate_restrictions(expression):
    110     "Propagate restriction nodes to wrap terminal objects directly."
--> 111     return apply_transformer(expression, RestrictionPropagator(), domain_type=Measure.INTERIOR_FACET)
    112 
    113 def check_restrictions(expression, require_restriction):

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in apply_transformer(e, transformer, domain_type)
    275     def _transform(expr):
    276         return transformer.visit(expr)
--> 277     return transform_integrands(e, _transform, domain_type)
    278 
    279 def ufl2ufl(e):

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in transform_integrands(form, transform, domain_type)
    250             integrand = itg.integrand()
    251             if domain_type is None or domain_type == itg.domain_type():
--> 252                 integrand = transform(integrand)
    253             if not isinstance(integrand, Zero):
    254                 newitg = itg.reconstruct(integrand)

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in _transform(expr)
    274     expression in form, or to form if it is an Expr."""
    275     def _transform(expr):
--> 276         return transformer.visit(expr)
    277     return transform_integrands(e, _transform, domain_type)
    278 

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in visit(self, o)
     99         if visit_children_first:
    100             # Yes, visit all children first and then call h.
--> 101             r = h(o, *map(self.visit, o.operands()))
    102         else:
    103             # No, this is a handler that handles its own children

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in visit(self, o)
    103             # No, this is a handler that handles its own children
    104             # (arguments self and o, where self is already bound)
--> 105             r = h(o)
    106 
    107         # Update stack and return

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/algorithms/propagate_restrictions.pyc in form_argument(self, o)
     67 
     68     def form_argument(self, o):
---> 69         ufl_assert(self.current_restriction is not None, "Form argument must be restricted.")
     70         #if self.current_restriction is None:
     71         #    return o

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/assertions.pyc in ufl_assert(condition, *message)
     35 def ufl_assert(condition, *message):
     36     "Assert that condition is true and otherwise issue an error with given message."
---> 37     if not condition: error(*message)
     38 

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/ufl/log.pyc in error(self, *message)
    152         "Write error message and raise an exception."
    153         self._log.error(*message)
--> 154         raise self._exception_type(self._format_raw(*message))
    155 
    156     def begin(self, *message):

UFLException: Form argument must be restricted.

In [66]:

    

n = FacetNormal(mesh)


    

In [70]:

    

di


    

    




---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-70-b0fe07840fc7> in <module>()
----> 1 vertex_to_dof_map()

/home/mwathen/Work/FEniCS/lib/python2.7/site-packages/dolfin/cpp/fem.pyc in vertex_to_dof_map(*args)
    722 
    723     """
--> 724   return _fem.vertex_to_dof_map(*args)
    725 class Equation(object):
    726     """

TypeError: vertex_to_dof_map() takes exactly one argument (0 given)

In [24]:

    

for cell in cells(mesh):
#   print "cell", cell.index(), "has edges :",
  for edge in edges(cell):
    print edge.index(),
    edge.
    print 
  pri


    

    




0
1
2

3
1
4

5
0
6

5
3
7

8
9
10

11
9
6

12
8
13

12
11
14

15
16
17

18
16
13

19
15
20

19
18
21

22
23
24

25
23
20

26
22
27

26
25
28

29
30
7

31
30
32

33
29
34

33
31
35

36
37
14

38
37
34

39
36
40

39
38
41

42
43
21

44
43
40

45
42
46

45
44
47

48
49
28

50
49
46

51
48
52

51
50
53

54
55
35

56
55
57

58
54
59

58
56
60

61
62
41

63
62
59

64
61
65

64
63
66

67
68
47

69
68
65

70
67
71

70
69
72

73
74
53

75
74
71

76
73
77

76
75
78

79
80
60

81
80
82

83
79
84

83
81
85

86
87
66

88
87
84

89
86
90

89
88
91

92
93
72

94
93
90

95
92
96

95
94
97

98
99
78

100
99
96

101
98
102

101
100
103

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    




    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.dot


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

print edge.global_index


    

In [ ]:

    

n = FacetNormal


    

In [ ]:

    




    

In [64]:

    

Magnetic.


    

    
Out[64]:





FiniteElement('Nedelec 1st kind H(curl)', Domain(Cell('triangle', 2), 'triangle_multiverse', 2, 2), 2, None)

In [ ]:

    




    

In [46]:

    

Magnetic.interpolate(DGu)


    

    




---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-46-dfc262124d54> in <module>()
----> 1 Magnetic.interpolate(grad(u))

TypeError: FunctionSpace_interpolate expected 3 arguments, got 2

In [52]:

    

b = Function(Magnetic)


    

In [55]:

    




    

    




---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
<ipython-input-55-41c132b404d7> in <module>()
----> 1 DGu.assign(grad(u))

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/functions/function.py in assign(self, rhs)
    408             cpp.dolfin_error("function.py",
    409                              "function assignment",
--> 410                              "Expects a Function or linear combinations of "\
    411                              "Functions in the same FunctionSpaces")
    412 

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/cpp/common.py in dolfin_error(*args)
   2294 
   2295     """
-> 2296   return _common.dolfin_error(*args)
   2297 
   2298 def deprecation(*args):

RuntimeError: 

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
***     fenics@fenicsproject.org
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error:   Unable to function assignment.
*** Reason:  Expects a Function or linear combinations of Functions in the same FunctionSpaces.
*** Where:   This error was encountered inside function.py.
*** Process: 0
*** 
*** DOLFIN version: 1.3.0
*** Git changeset:  
*** -------------------------------------------------------------------------

In [49]:

    

b = grad(u)


    

In [50]:

    

b


    

    
Out[50]:





$\mathbf{grad}{\left({w_h^8}\right)}$

In [51]:

    

u


    

    
Out[51]:





${w_h^8}$

In [ ]:

    




    

In [32]:

    

uu.reconstruct?


    

In [ ]:

    




    

In [ ]:

    

T = as_vector


    

In [16]:

    

T = as_vector([N[1],-N[0]])


    

In [18]:

    

u = Function(LagrangeL)


    

In [21]:

    

assemble(inner(inner(uu,T),u)*dS)


    

    




Form argument must be restricted.






    




ERROR:UFL:Form argument must be restricted.






    




---------------------------------------------------------------------------
UFLException                              Traceback (most recent call last)
<ipython-input-21-fc5978e91932> in <module>()
----> 1 assemble(inner(inner(uu,T),u)*dS)

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/fem/assembling.py in assemble(form, tensor, mesh, coefficients, function_spaces, cell_domains, exterior_facet_domains, interior_facet_domains, reset_sparsity, add_values, finalize_tensor, keep_diagonal, backend, form_compiler_parameters, bcs)
    166                            coefficients=coefficients,
    167                            form_compiler_parameters=form_compiler_parameters,
--> 168                            common_cell=common_cell)
    169 
    170     # Set mesh if specified (important for functionals without a function spaces)

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/fem/form.py in __init__(self, form, function_spaces, coefficients, subdomains, form_compiler_parameters, common_cell)
     54             self._compiled_form, module, self.form_data, prefix \
     55                                  = jit(form, form_compiler_parameters,
---> 56                                        common_cell)
     57 
     58         elif isinstance(form, ufc.form):

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/compilemodules/jit.py in mpi_jit(*args, **kwargs)
     58         # Just call JIT compiler when running in serial
     59         if MPI.num_processes() == 1:
---> 60             return local_jit(*args, **kwargs)
     61 
     62         # Compile first on process 0

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/dolfin/compilemodules/jit.py in jit(form, form_compiler_parameters, common_cell)
    120         raise RuntimeError, "Form compiler must implement the jit function."
    121 
--> 122     return jit_compile(form, parameters=p, common_cell=common_cell)

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ffc/jitcompiler.pyc in jit(ufl_object, parameters, common_cell)
     76         return jit_element(ufl_object, parameters)
     77     else:
---> 78         return jit_form(ufl_object, parameters, common_cell)
     79 
     80 def _auto_select_degree(elements):

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ffc/jitcompiler.pyc in jit_form(form, parameters, common_cell)
    167     # Compute form metadata and extract preprocessed form
    168     form_data = form.compute_form_data(common_cell=common_cell,
--> 169                                        element_mapping=element_mapping)
    170 
    171     # Wrap input

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/form.pyc in compute_form_data(self, object_names, common_cell, element_mapping)
    122                                          object_names=object_names,
    123                                          common_cell=common_cell,
--> 124                                          element_mapping=element_mapping)
    125         # Always validate arguments, keeping sure that the validation works
    126         self._form_data.validate(object_names=object_names,

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/preprocess.pyc in preprocess(form, object_names, common_cell, element_mapping)
    116 
    117     # Propagate restrictions of interior facet integrals to the terminal nodes
--> 118     form = propagate_restrictions(form) # INTEGRAL, EXPR
    119 
    120     # --- BEGIN DOMAIN SPLITTING AND JOINING

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/propagate_restrictions.pyc in propagate_restrictions(expression)
    109 def propagate_restrictions(expression):
    110     "Propagate restriction nodes to wrap terminal objects directly."
--> 111     return apply_transformer(expression, RestrictionPropagator(), domain_type=Measure.INTERIOR_FACET)
    112 
    113 def check_restrictions(expression, require_restriction):

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in apply_transformer(e, transformer, domain_type)
    275     def _transform(expr):
    276         return transformer.visit(expr)
--> 277     return transform_integrands(e, _transform, domain_type)
    278 
    279 def ufl2ufl(e):

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in transform_integrands(form, transform, domain_type)
    250             integrand = itg.integrand()
    251             if domain_type is None or domain_type == itg.domain_type():
--> 252                 integrand = transform(integrand)
    253             if not isinstance(integrand, Zero):
    254                 newitg = itg.reconstruct(integrand)

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in _transform(expr)
    274     expression in form, or to form if it is an Expr."""
    275     def _transform(expr):
--> 276         return transformer.visit(expr)
    277     return transform_integrands(e, _transform, domain_type)
    278 

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in visit(self, o)
     99         if visit_children_first:
    100             # Yes, visit all children first and then call h.
--> 101             r = h(o, *map(self.visit, o.operands()))
    102         else:
    103             # No, this is a handler that handles its own children

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/transformer.pyc in visit(self, o)
    103             # No, this is a handler that handles its own children
    104             # (arguments self and o, where self is already bound)
--> 105             r = h(o)
    106 
    107         # Update stack and return

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/algorithms/propagate_restrictions.pyc in form_argument(self, o)
     67 
     68     def form_argument(self, o):
---> 69         ufl_assert(self.current_restriction is not None, "Form argument must be restricted.")
     70         #if self.current_restriction is None:
     71         #    return o

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/assertions.pyc in ufl_assert(condition, *message)
     35 def ufl_assert(condition, *message):
     36     "Assert that condition is true and otherwise issue an error with given message."
---> 37     if not condition: error(*message)
     38 

/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/ufl/log.pyc in error(self, *message)
    152         "Write error message and raise an exception."
    153         self._log.error(*message)
--> 154         raise self._exception_type(self._format_raw(*message))
    155 
    156     def begin(self, *message):

UFLException: Form argument must be restricted.

In [24]:

    

dof = Magnetic.dofmap()


    

In [25]:

    

dof.


    

In [27]:

    

mesh.cells()


    

    
Out[27]:





array([[ 0,  1, 25],
       [ 0,  5, 25],
       [ 1,  6, 25],
       [ 5,  6, 25],
       [ 1,  2, 26],
       [ 1,  6, 26],
       [ 2,  7, 26],
       [ 6,  7, 26],
       [ 2,  3, 27],
       [ 2,  7, 27],
       [ 3,  8, 27],
       [ 7,  8, 27],
       [ 3,  4, 28],
       [ 3,  8, 28],
       [ 4,  9, 28],
       [ 8,  9, 28],
       [ 5,  6, 29],
       [ 5, 10, 29],
       [ 6, 11, 29],
       [10, 11, 29],
       [ 6,  7, 30],
       [ 6, 11, 30],
       [ 7, 12, 30],
       [11, 12, 30],
       [ 7,  8, 31],
       [ 7, 12, 31],
       [ 8, 13, 31],
       [12, 13, 31],
       [ 8,  9, 32],
       [ 8, 13, 32],
       [ 9, 14, 32],
       [13, 14, 32],
       [10, 11, 33],
       [10, 15, 33],
       [11, 16, 33],
       [15, 16, 33],
       [11, 12, 34],
       [11, 16, 34],
       [12, 17, 34],
       [16, 17, 34],
       [12, 13, 35],
       [12, 17, 35],
       [13, 18, 35],
       [17, 18, 35],
       [13, 14, 36],
       [13, 18, 36],
       [14, 19, 36],
       [18, 19, 36],
       [15, 16, 37],
       [15, 20, 37],
       [16, 21, 37],
       [20, 21, 37],
       [16, 17, 38],
       [16, 21, 38],
       [17, 22, 38],
       [21, 22, 38],
       [17, 18, 39],
       [17, 22, 39],
       [18, 23, 39],
       [22, 23, 39],
       [18, 19, 40],
       [18, 23, 40],
       [19, 24, 40],
       [23, 24, 40]], dtype=uint32)

In [ ]: