@@ -46,22 +46,24 @@ function linearelasticityweakform(deltadu::Array{Float64}, w::Array{Float64})
46
46
end
47
47
48
48
# linear elasticity operators
49
- function makeoperator (basis:: TensorBasis )
50
- inputs = [
51
- OperatorField (basis, [EvaluationMode. gradient], " gradent of deformation" ),
52
- OperatorField (basis, [EvaluationMode. quadratureweights], " quadrature weights" ),
53
- ]
54
- outputs = [
55
- OperatorField (
56
- basis,
57
- [EvaluationMode. gradient],
58
- " test function gradient of deformation" ,
59
- ),
60
- ]
61
- return Operator (linearelasticityweakform, mesh, inputs, outputs)
49
+ begin
50
+ local function makeoperator (basis:: TensorBasis )
51
+ inputs = [
52
+ OperatorField (basis, [EvaluationMode. gradient], " gradent of deformation" ),
53
+ OperatorField (basis, [EvaluationMode. quadratureweights], " quadrature weights" ),
54
+ ]
55
+ outputs = [
56
+ OperatorField (
57
+ basis,
58
+ [EvaluationMode. gradient],
59
+ " test function gradient of deformation" ,
60
+ ),
61
+ ]
62
+ return Operator (linearelasticityweakform, mesh, inputs, outputs)
63
+ end
64
+ fineoperator = makeoperator (finebasis)
65
+ coarseoperator = makeoperator (coarsebasis)
62
66
end
63
- fineoperator = makeoperator (finebasis)
64
- coarseoperator = makeoperator (coarsebasis)
65
67
66
68
# Chebyshev smoother
67
69
chebyshev = Chebyshev (fineoperator)
0 commit comments