Changeset 0260d3 for test/unit_test/library/force_sinus.cpp
- Timestamp:
- Jul 5, 2012, 11:54:39 AM (13 years ago)
- Children:
- 185795
- Parents:
- ef94e7
- File:
-
- 1 edited
-
test/unit_test/library/force_sinus.cpp (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
test/unit_test/library/force_sinus.cpp
ref94e7 r0260d3 131 131 132 132 } 133 134 pos.X() = 0.25;135 pos.Y() = 0.25;136 for (pos.Z() = 0.25; pos.Z()<=0.75; pos.Z() += 0.001) {137 p.Pos() = pos;138 139 const Index index = (p.Pos()-sol.Extent().Begin())/sol.Extent().MeshWidth()-sol.Global().LocalBegin()+sol.Local().Begin();140 141 ip.ComputeCoefficients(sol, index);142 ip.Evaluate(p);143 144 Vector ref = Vector(std::cos(sine_factor*p.Pos().X())*std::sin(sine_factor*p.Pos().Y())*std::sin(sine_factor*p.Pos().Z()),145 std::sin(sine_factor*p.Pos().X())*std::cos(sine_factor*p.Pos().Y())*std::sin(sine_factor*p.Pos().Z()),146 std::sin(sine_factor*p.Pos().X())*std::sin(sine_factor*p.Pos().Y())*std::cos(sine_factor*p.Pos().Z()));147 ref *= -1.0 * sine_factor;148 std::cout << pos.Z() << " " << std::abs((ref.Z() - p.Field().Z()) / ref.Z()) << std::endl;149 }150 133 } 151 134
Note:
See TracChangeset
for help on using the changeset viewer.
