|
99 | 99 |
|
100 | 100 |
|
101 | 101 | ke=Kinematics(Electro, Solid) |
102 | | - km=Kinematics(Mechano,Solid) |
| 102 | + km=Kinematics(Mechano, Solid) |
103 | 103 | k=(km,ke) |
104 | 104 | function jach(uh, φh) |
105 | | - jac((du, dφ), (v, vφ)) = jacobian(modelelectro,k, (uh, φh), (du, dφ), (v, vφ), dΩ) |
| 105 | + jac((du, dφ), (v, vφ)) = jacobian(modelelectro, k, (uh, φh), (du, dφ), (v, vφ), dΩ) |
106 | 106 | end |
107 | 107 | function jac_mech(uh, φh) |
108 | | - jac(du, v) = jacobian(modelelectro, Mechano, k,(uh, φh), du, v, dΩ) |
| 108 | + jac(du, v) = jacobian(modelelectro, Mechano, k, (uh, φh), du, v, dΩ) |
109 | 109 | end |
110 | 110 | function jac_elech(uh, φh) |
111 | | - jac(dφ, vφ) = jacobian(modelelectro, Electro, k,(uh, φh), dφ, vφ, dΩ) |
| 111 | + jac(dφ, vφ) = jacobian(modelelectro, Electro, k, (uh, φh), dφ, vφ, dΩ) |
112 | 112 | end |
113 | 113 |
|
114 | 114 | jac_ = assemble_matrix(jach(uh, φh), V, V) |
115 | 115 | jac_m = assemble_matrix(jac_mech(uh, φh), Vu, Vu) |
116 | 116 | jac_e = assemble_matrix(jac_elech(uh, φh), Vφ, Vφ) |
117 | | - ≈ |
| 117 | + |
118 | 118 | @test norm(jac_) ≈ 18.934585248125135 |
119 | 119 | @test jac_[1] ≈ 0.7777777777777775 |
120 | 120 | @test jac_[end] ≈ -1.3333333333333326 |
|
0 commit comments