diff --git a/Data/Mesh/nonlinearmembrane_4elts_rectangle.mesh b/Data/Mesh/nonlinearmembrane_4elts_rectangle.mesh new file mode 100644 index 0000000000000000000000000000000000000000..fb2b81475214097b89d6c25ac52294da690555fb --- /dev/null +++ b/Data/Mesh/nonlinearmembrane_4elts_rectangle.mesh @@ -0,0 +1,31 @@ +MeshVersionFormatted 2 + +Dimension +3 + +Vertices +6 +0 0 0 1 +1 0 0 2 +2 0 0 3 +0 2 0 4 +1 2 0 5 +2 2 0 6 + +Edges +6 +1 2 11 +2 3 12 +3 6 13 +6 5 14 +5 4 15 +4 1 16 + +Triangles +4 +1 2 4 21 +5 4 2 22 +2 3 5 23 +6 5 3 24 + +End \ No newline at end of file diff --git a/Data/Mesh/nonlinearmembrane_4elts_rectangle_random_basis.mesh b/Data/Mesh/nonlinearmembrane_4elts_rectangle_random_basis.mesh new file mode 100644 index 0000000000000000000000000000000000000000..c7e43da0933498a951bec94734b1429682213051 --- /dev/null +++ b/Data/Mesh/nonlinearmembrane_4elts_rectangle_random_basis.mesh @@ -0,0 +1,31 @@ +MeshVersionFormatted 2 + +Dimension +3 + +Vertices +6 +0 0 0 1 +1 0 -0.5 2 +2 0 0 3 +0 2 0 4 +1 2 -0.5 5 +2 2 0 6 + +Edges +6 +1 2 11 +2 3 12 +3 6 13 +6 5 14 +5 4 15 +4 1 16 + +Triangles +4 +4 1 2 21 +2 5 4 22 +5 2 3 23 +3 6 5 24 + +End \ No newline at end of file diff --git a/Data/Mesh/nonlinearmembrane_quartercylinder.mesh b/Data/Mesh/nonlinearmembrane_quartercylinder.mesh new file mode 100644 index 0000000000000000000000000000000000000000..fa1d8044ebdc5f08635b202e467cfee4b2c4d89d --- /dev/null +++ b/Data/Mesh/nonlinearmembrane_quartercylinder.mesh @@ -0,0 +1,3225 @@ +MeshVersionFormatted 2 + +Dimension +3 +Vertices +1071 +1.0 0.0 -1 1 +0.999506560366 0.0314107590781 -1 0 +0.998026728428 0.0627905195293 -1 0 +0.995561964603 0.0941083133185 -1 0 +0.992114701314 0.125333233564 -1 0 +0.987688340595 0.15643446504 -1 0 +0.982287250729 0.187381314586 -1 0 +0.975916761939 0.218143241397 -1 0 +0.968583161129 0.248689887165 -1 0 +0.960293685677 0.278991106039 -1 0 +0.951056516295 0.309016994375 -1 0 +0.940880768954 0.338737920245 -1 0 +0.929776485888 0.368124552685 -1 0 +0.917754625684 0.397147890635 -1 0 +0.904827052466 0.425779291565 -1 0 +0.891006524188 0.45399049974 -1 0 +0.876306680044 0.481753674102 -1 0 +0.860742027004 0.50904141575 -1 0 +0.844327925502 0.535826794979 -1 0 +0.827080574275 0.562083377852 -1 0 +0.809016994375 0.587785252292 -1 0 +0.790155012376 0.612907053653 -1 0 +0.770513242776 0.637423989749 -1 0 +0.75011106963 0.661311865324 -1 0 +0.728968627421 0.684547105929 -1 0 +0.707106781187 0.707106781187 -1 0 +0.684547105929 0.728968627421 -1 0 +0.661311865324 0.75011106963 -1 0 +0.637423989749 0.770513242776 -1 0 +0.612907053653 0.790155012376 -1 0 +0.587785252292 0.809016994375 -1 0 +0.562083377852 0.827080574275 -1 0 +0.535826794979 0.844327925502 -1 0 +0.50904141575 0.860742027004 -1 0 +0.481753674102 0.876306680044 -1 0 +0.45399049974 0.891006524188 -1 0 +0.425779291565 0.904827052466 -1 0 +0.397147890635 0.917754625684 -1 0 +0.368124552685 0.929776485888 -1 0 +0.338737920245 0.940880768954 -1 0 +0.309016994375 0.951056516295 -1 0 +0.278991106039 0.960293685677 -1 0 +0.248689887165 0.968583161129 -1 0 +0.218143241397 0.975916761939 -1 0 +0.187381314586 0.982287250729 -1 0 +0.15643446504 0.987688340595 -1 0 +0.125333233564 0.992114701314 -1 0 +0.0941083133185 0.995561964603 -1 0 +0.0627905195293 0.998026728428 -1 0 +0.0314107590781 0.999506560366 -1 0 +6.12323399574e-17 1.0 -1 2 +1.0 0.0 -0.9 1 +0.999506560366 0.0314107590781 -0.9 0 +0.998026728428 0.0627905195293 -0.9 0 +0.995561964603 0.0941083133185 -0.9 0 +0.992114701314 0.125333233564 -0.9 0 +0.987688340595 0.15643446504 -0.9 0 +0.982287250729 0.187381314586 -0.9 0 +0.975916761939 0.218143241397 -0.9 0 +0.968583161129 0.248689887165 -0.9 0 +0.960293685677 0.278991106039 -0.9 0 +0.951056516295 0.309016994375 -0.9 0 +0.940880768954 0.338737920245 -0.9 0 +0.929776485888 0.368124552685 -0.9 0 +0.917754625684 0.397147890635 -0.9 0 +0.904827052466 0.425779291565 -0.9 0 +0.891006524188 0.45399049974 -0.9 0 +0.876306680044 0.481753674102 -0.9 0 +0.860742027004 0.50904141575 -0.9 0 +0.844327925502 0.535826794979 -0.9 0 +0.827080574275 0.562083377852 -0.9 0 +0.809016994375 0.587785252292 -0.9 0 +0.790155012376 0.612907053653 -0.9 0 +0.770513242776 0.637423989749 -0.9 0 +0.75011106963 0.661311865324 -0.9 0 +0.728968627421 0.684547105929 -0.9 0 +0.707106781187 0.707106781187 -0.9 0 +0.684547105929 0.728968627421 -0.9 0 +0.661311865324 0.75011106963 -0.9 0 +0.637423989749 0.770513242776 -0.9 0 +0.612907053653 0.790155012376 -0.9 0 +0.587785252292 0.809016994375 -0.9 0 +0.562083377852 0.827080574275 -0.9 0 +0.535826794979 0.844327925502 -0.9 0 +0.50904141575 0.860742027004 -0.9 0 +0.481753674102 0.876306680044 -0.9 0 +0.45399049974 0.891006524188 -0.9 0 +0.425779291565 0.904827052466 -0.9 0 +0.397147890635 0.917754625684 -0.9 0 +0.368124552685 0.929776485888 -0.9 0 +0.338737920245 0.940880768954 -0.9 0 +0.309016994375 0.951056516295 -0.9 0 +0.278991106039 0.960293685677 -0.9 0 +0.248689887165 0.968583161129 -0.9 0 +0.218143241397 0.975916761939 -0.9 0 +0.187381314586 0.982287250729 -0.9 0 +0.15643446504 0.987688340595 -0.9 0 +0.125333233564 0.992114701314 -0.9 0 +0.0941083133185 0.995561964603 -0.9 0 +0.0627905195293 0.998026728428 -0.9 0 +0.0314107590781 0.999506560366 -0.9 0 +6.12323399574e-17 1.0 -0.9 2 +1.0 0.0 -0.8 1 +0.999506560366 0.0314107590781 -0.8 0 +0.998026728428 0.0627905195293 -0.8 0 +0.995561964603 0.0941083133185 -0.8 0 +0.992114701314 0.125333233564 -0.8 0 +0.987688340595 0.15643446504 -0.8 0 +0.982287250729 0.187381314586 -0.8 0 +0.975916761939 0.218143241397 -0.8 0 +0.968583161129 0.248689887165 -0.8 0 +0.960293685677 0.278991106039 -0.8 0 +0.951056516295 0.309016994375 -0.8 0 +0.940880768954 0.338737920245 -0.8 0 +0.929776485888 0.368124552685 -0.8 0 +0.917754625684 0.397147890635 -0.8 0 +0.904827052466 0.425779291565 -0.8 0 +0.891006524188 0.45399049974 -0.8 0 +0.876306680044 0.481753674102 -0.8 0 +0.860742027004 0.50904141575 -0.8 0 +0.844327925502 0.535826794979 -0.8 0 +0.827080574275 0.562083377852 -0.8 0 +0.809016994375 0.587785252292 -0.8 0 +0.790155012376 0.612907053653 -0.8 0 +0.770513242776 0.637423989749 -0.8 0 +0.75011106963 0.661311865324 -0.8 0 +0.728968627421 0.684547105929 -0.8 0 +0.707106781187 0.707106781187 -0.8 0 +0.684547105929 0.728968627421 -0.8 0 +0.661311865324 0.75011106963 -0.8 0 +0.637423989749 0.770513242776 -0.8 0 +0.612907053653 0.790155012376 -0.8 0 +0.587785252292 0.809016994375 -0.8 0 +0.562083377852 0.827080574275 -0.8 0 +0.535826794979 0.844327925502 -0.8 0 +0.50904141575 0.860742027004 -0.8 0 +0.481753674102 0.876306680044 -0.8 0 +0.45399049974 0.891006524188 -0.8 0 +0.425779291565 0.904827052466 -0.8 0 +0.397147890635 0.917754625684 -0.8 0 +0.368124552685 0.929776485888 -0.8 0 +0.338737920245 0.940880768954 -0.8 0 +0.309016994375 0.951056516295 -0.8 0 +0.278991106039 0.960293685677 -0.8 0 +0.248689887165 0.968583161129 -0.8 0 +0.218143241397 0.975916761939 -0.8 0 +0.187381314586 0.982287250729 -0.8 0 +0.15643446504 0.987688340595 -0.8 0 +0.125333233564 0.992114701314 -0.8 0 +0.0941083133185 0.995561964603 -0.8 0 +0.0627905195293 0.998026728428 -0.8 0 +0.0314107590781 0.999506560366 -0.8 0 +6.12323399574e-17 1.0 -0.8 2 +1.0 0.0 -0.7 1 +0.999506560366 0.0314107590781 -0.7 0 +0.998026728428 0.0627905195293 -0.7 0 +0.995561964603 0.0941083133185 -0.7 0 +0.992114701314 0.125333233564 -0.7 0 +0.987688340595 0.15643446504 -0.7 0 +0.982287250729 0.187381314586 -0.7 0 +0.975916761939 0.218143241397 -0.7 0 +0.968583161129 0.248689887165 -0.7 0 +0.960293685677 0.278991106039 -0.7 0 +0.951056516295 0.309016994375 -0.7 0 +0.940880768954 0.338737920245 -0.7 0 +0.929776485888 0.368124552685 -0.7 0 +0.917754625684 0.397147890635 -0.7 0 +0.904827052466 0.425779291565 -0.7 0 +0.891006524188 0.45399049974 -0.7 0 +0.876306680044 0.481753674102 -0.7 0 +0.860742027004 0.50904141575 -0.7 0 +0.844327925502 0.535826794979 -0.7 0 +0.827080574275 0.562083377852 -0.7 0 +0.809016994375 0.587785252292 -0.7 0 +0.790155012376 0.612907053653 -0.7 0 +0.770513242776 0.637423989749 -0.7 0 +0.75011106963 0.661311865324 -0.7 0 +0.728968627421 0.684547105929 -0.7 0 +0.707106781187 0.707106781187 -0.7 0 +0.684547105929 0.728968627421 -0.7 0 +0.661311865324 0.75011106963 -0.7 0 +0.637423989749 0.770513242776 -0.7 0 +0.612907053653 0.790155012376 -0.7 0 +0.587785252292 0.809016994375 -0.7 0 +0.562083377852 0.827080574275 -0.7 0 +0.535826794979 0.844327925502 -0.7 0 +0.50904141575 0.860742027004 -0.7 0 +0.481753674102 0.876306680044 -0.7 0 +0.45399049974 0.891006524188 -0.7 0 +0.425779291565 0.904827052466 -0.7 0 +0.397147890635 0.917754625684 -0.7 0 +0.368124552685 0.929776485888 -0.7 0 +0.338737920245 0.940880768954 -0.7 0 +0.309016994375 0.951056516295 -0.7 0 +0.278991106039 0.960293685677 -0.7 0 +0.248689887165 0.968583161129 -0.7 0 +0.218143241397 0.975916761939 -0.7 0 +0.187381314586 0.982287250729 -0.7 0 +0.15643446504 0.987688340595 -0.7 0 +0.125333233564 0.992114701314 -0.7 0 +0.0941083133185 0.995561964603 -0.7 0 +0.0627905195293 0.998026728428 -0.7 0 +0.0314107590781 0.999506560366 -0.7 0 +6.12323399574e-17 1.0 -0.7 2 +1.0 0.0 -0.6 1 +0.999506560366 0.0314107590781 -0.6 0 +0.998026728428 0.0627905195293 -0.6 0 +0.995561964603 0.0941083133185 -0.6 0 +0.992114701314 0.125333233564 -0.6 0 +0.987688340595 0.15643446504 -0.6 0 +0.982287250729 0.187381314586 -0.6 0 +0.975916761939 0.218143241397 -0.6 0 +0.968583161129 0.248689887165 -0.6 0 +0.960293685677 0.278991106039 -0.6 0 +0.951056516295 0.309016994375 -0.6 0 +0.940880768954 0.338737920245 -0.6 0 +0.929776485888 0.368124552685 -0.6 0 +0.917754625684 0.397147890635 -0.6 0 +0.904827052466 0.425779291565 -0.6 0 +0.891006524188 0.45399049974 -0.6 0 +0.876306680044 0.481753674102 -0.6 0 +0.860742027004 0.50904141575 -0.6 0 +0.844327925502 0.535826794979 -0.6 0 +0.827080574275 0.562083377852 -0.6 0 +0.809016994375 0.587785252292 -0.6 0 +0.790155012376 0.612907053653 -0.6 0 +0.770513242776 0.637423989749 -0.6 0 +0.75011106963 0.661311865324 -0.6 0 +0.728968627421 0.684547105929 -0.6 0 +0.707106781187 0.707106781187 -0.6 0 +0.684547105929 0.728968627421 -0.6 0 +0.661311865324 0.75011106963 -0.6 0 +0.637423989749 0.770513242776 -0.6 0 +0.612907053653 0.790155012376 -0.6 0 +0.587785252292 0.809016994375 -0.6 0 +0.562083377852 0.827080574275 -0.6 0 +0.535826794979 0.844327925502 -0.6 0 +0.50904141575 0.860742027004 -0.6 0 +0.481753674102 0.876306680044 -0.6 0 +0.45399049974 0.891006524188 -0.6 0 +0.425779291565 0.904827052466 -0.6 0 +0.397147890635 0.917754625684 -0.6 0 +0.368124552685 0.929776485888 -0.6 0 +0.338737920245 0.940880768954 -0.6 0 +0.309016994375 0.951056516295 -0.6 0 +0.278991106039 0.960293685677 -0.6 0 +0.248689887165 0.968583161129 -0.6 0 +0.218143241397 0.975916761939 -0.6 0 +0.187381314586 0.982287250729 -0.6 0 +0.15643446504 0.987688340595 -0.6 0 +0.125333233564 0.992114701314 -0.6 0 +0.0941083133185 0.995561964603 -0.6 0 +0.0627905195293 0.998026728428 -0.6 0 +0.0314107590781 0.999506560366 -0.6 0 +6.12323399574e-17 1.0 -0.6 2 +1.0 0.0 -0.5 1 +0.999506560366 0.0314107590781 -0.5 0 +0.998026728428 0.0627905195293 -0.5 0 +0.995561964603 0.0941083133185 -0.5 0 +0.992114701314 0.125333233564 -0.5 0 +0.987688340595 0.15643446504 -0.5 0 +0.982287250729 0.187381314586 -0.5 0 +0.975916761939 0.218143241397 -0.5 0 +0.968583161129 0.248689887165 -0.5 0 +0.960293685677 0.278991106039 -0.5 0 +0.951056516295 0.309016994375 -0.5 0 +0.940880768954 0.338737920245 -0.5 0 +0.929776485888 0.368124552685 -0.5 0 +0.917754625684 0.397147890635 -0.5 0 +0.904827052466 0.425779291565 -0.5 0 +0.891006524188 0.45399049974 -0.5 0 +0.876306680044 0.481753674102 -0.5 0 +0.860742027004 0.50904141575 -0.5 0 +0.844327925502 0.535826794979 -0.5 0 +0.827080574275 0.562083377852 -0.5 0 +0.809016994375 0.587785252292 -0.5 0 +0.790155012376 0.612907053653 -0.5 0 +0.770513242776 0.637423989749 -0.5 0 +0.75011106963 0.661311865324 -0.5 0 +0.728968627421 0.684547105929 -0.5 0 +0.707106781187 0.707106781187 -0.5 0 +0.684547105929 0.728968627421 -0.5 0 +0.661311865324 0.75011106963 -0.5 0 +0.637423989749 0.770513242776 -0.5 0 +0.612907053653 0.790155012376 -0.5 0 +0.587785252292 0.809016994375 -0.5 0 +0.562083377852 0.827080574275 -0.5 0 +0.535826794979 0.844327925502 -0.5 0 +0.50904141575 0.860742027004 -0.5 0 +0.481753674102 0.876306680044 -0.5 0 +0.45399049974 0.891006524188 -0.5 0 +0.425779291565 0.904827052466 -0.5 0 +0.397147890635 0.917754625684 -0.5 0 +0.368124552685 0.929776485888 -0.5 0 +0.338737920245 0.940880768954 -0.5 0 +0.309016994375 0.951056516295 -0.5 0 +0.278991106039 0.960293685677 -0.5 0 +0.248689887165 0.968583161129 -0.5 0 +0.218143241397 0.975916761939 -0.5 0 +0.187381314586 0.982287250729 -0.5 0 +0.15643446504 0.987688340595 -0.5 0 +0.125333233564 0.992114701314 -0.5 0 +0.0941083133185 0.995561964603 -0.5 0 +0.0627905195293 0.998026728428 -0.5 0 +0.0314107590781 0.999506560366 -0.5 0 +6.12323399574e-17 1.0 -0.5 2 +1.0 0.0 -0.4 1 +0.999506560366 0.0314107590781 -0.4 0 +0.998026728428 0.0627905195293 -0.4 0 +0.995561964603 0.0941083133185 -0.4 0 +0.992114701314 0.125333233564 -0.4 0 +0.987688340595 0.15643446504 -0.4 0 +0.982287250729 0.187381314586 -0.4 0 +0.975916761939 0.218143241397 -0.4 0 +0.968583161129 0.248689887165 -0.4 0 +0.960293685677 0.278991106039 -0.4 0 +0.951056516295 0.309016994375 -0.4 0 +0.940880768954 0.338737920245 -0.4 0 +0.929776485888 0.368124552685 -0.4 0 +0.917754625684 0.397147890635 -0.4 0 +0.904827052466 0.425779291565 -0.4 0 +0.891006524188 0.45399049974 -0.4 0 +0.876306680044 0.481753674102 -0.4 0 +0.860742027004 0.50904141575 -0.4 0 +0.844327925502 0.535826794979 -0.4 0 +0.827080574275 0.562083377852 -0.4 0 +0.809016994375 0.587785252292 -0.4 0 +0.790155012376 0.612907053653 -0.4 0 +0.770513242776 0.637423989749 -0.4 0 +0.75011106963 0.661311865324 -0.4 0 +0.728968627421 0.684547105929 -0.4 0 +0.707106781187 0.707106781187 -0.4 0 +0.684547105929 0.728968627421 -0.4 0 +0.661311865324 0.75011106963 -0.4 0 +0.637423989749 0.770513242776 -0.4 0 +0.612907053653 0.790155012376 -0.4 0 +0.587785252292 0.809016994375 -0.4 0 +0.562083377852 0.827080574275 -0.4 0 +0.535826794979 0.844327925502 -0.4 0 +0.50904141575 0.860742027004 -0.4 0 +0.481753674102 0.876306680044 -0.4 0 +0.45399049974 0.891006524188 -0.4 0 +0.425779291565 0.904827052466 -0.4 0 +0.397147890635 0.917754625684 -0.4 0 +0.368124552685 0.929776485888 -0.4 0 +0.338737920245 0.940880768954 -0.4 0 +0.309016994375 0.951056516295 -0.4 0 +0.278991106039 0.960293685677 -0.4 0 +0.248689887165 0.968583161129 -0.4 0 +0.218143241397 0.975916761939 -0.4 0 +0.187381314586 0.982287250729 -0.4 0 +0.15643446504 0.987688340595 -0.4 0 +0.125333233564 0.992114701314 -0.4 0 +0.0941083133185 0.995561964603 -0.4 0 +0.0627905195293 0.998026728428 -0.4 0 +0.0314107590781 0.999506560366 -0.4 0 +6.12323399574e-17 1.0 -0.4 2 +1.0 0.0 -0.3 1 +0.999506560366 0.0314107590781 -0.3 0 +0.998026728428 0.0627905195293 -0.3 0 +0.995561964603 0.0941083133185 -0.3 0 +0.992114701314 0.125333233564 -0.3 0 +0.987688340595 0.15643446504 -0.3 0 +0.982287250729 0.187381314586 -0.3 0 +0.975916761939 0.218143241397 -0.3 0 +0.968583161129 0.248689887165 -0.3 0 +0.960293685677 0.278991106039 -0.3 0 +0.951056516295 0.309016994375 -0.3 0 +0.940880768954 0.338737920245 -0.3 0 +0.929776485888 0.368124552685 -0.3 0 +0.917754625684 0.397147890635 -0.3 0 +0.904827052466 0.425779291565 -0.3 0 +0.891006524188 0.45399049974 -0.3 0 +0.876306680044 0.481753674102 -0.3 0 +0.860742027004 0.50904141575 -0.3 0 +0.844327925502 0.535826794979 -0.3 0 +0.827080574275 0.562083377852 -0.3 0 +0.809016994375 0.587785252292 -0.3 0 +0.790155012376 0.612907053653 -0.3 0 +0.770513242776 0.637423989749 -0.3 0 +0.75011106963 0.661311865324 -0.3 0 +0.728968627421 0.684547105929 -0.3 0 +0.707106781187 0.707106781187 -0.3 0 +0.684547105929 0.728968627421 -0.3 0 +0.661311865324 0.75011106963 -0.3 0 +0.637423989749 0.770513242776 -0.3 0 +0.612907053653 0.790155012376 -0.3 0 +0.587785252292 0.809016994375 -0.3 0 +0.562083377852 0.827080574275 -0.3 0 +0.535826794979 0.844327925502 -0.3 0 +0.50904141575 0.860742027004 -0.3 0 +0.481753674102 0.876306680044 -0.3 0 +0.45399049974 0.891006524188 -0.3 0 +0.425779291565 0.904827052466 -0.3 0 +0.397147890635 0.917754625684 -0.3 0 +0.368124552685 0.929776485888 -0.3 0 +0.338737920245 0.940880768954 -0.3 0 +0.309016994375 0.951056516295 -0.3 0 +0.278991106039 0.960293685677 -0.3 0 +0.248689887165 0.968583161129 -0.3 0 +0.218143241397 0.975916761939 -0.3 0 +0.187381314586 0.982287250729 -0.3 0 +0.15643446504 0.987688340595 -0.3 0 +0.125333233564 0.992114701314 -0.3 0 +0.0941083133185 0.995561964603 -0.3 0 +0.0627905195293 0.998026728428 -0.3 0 +0.0314107590781 0.999506560366 -0.3 0 +6.12323399574e-17 1.0 -0.3 2 +1.0 0.0 -0.2 1 +0.999506560366 0.0314107590781 -0.2 0 +0.998026728428 0.0627905195293 -0.2 0 +0.995561964603 0.0941083133185 -0.2 0 +0.992114701314 0.125333233564 -0.2 0 +0.987688340595 0.15643446504 -0.2 0 +0.982287250729 0.187381314586 -0.2 0 +0.975916761939 0.218143241397 -0.2 0 +0.968583161129 0.248689887165 -0.2 0 +0.960293685677 0.278991106039 -0.2 0 +0.951056516295 0.309016994375 -0.2 0 +0.940880768954 0.338737920245 -0.2 0 +0.929776485888 0.368124552685 -0.2 0 +0.917754625684 0.397147890635 -0.2 0 +0.904827052466 0.425779291565 -0.2 0 +0.891006524188 0.45399049974 -0.2 0 +0.876306680044 0.481753674102 -0.2 0 +0.860742027004 0.50904141575 -0.2 0 +0.844327925502 0.535826794979 -0.2 0 +0.827080574275 0.562083377852 -0.2 0 +0.809016994375 0.587785252292 -0.2 0 +0.790155012376 0.612907053653 -0.2 0 +0.770513242776 0.637423989749 -0.2 0 +0.75011106963 0.661311865324 -0.2 0 +0.728968627421 0.684547105929 -0.2 0 +0.707106781187 0.707106781187 -0.2 0 +0.684547105929 0.728968627421 -0.2 0 +0.661311865324 0.75011106963 -0.2 0 +0.637423989749 0.770513242776 -0.2 0 +0.612907053653 0.790155012376 -0.2 0 +0.587785252292 0.809016994375 -0.2 0 +0.562083377852 0.827080574275 -0.2 0 +0.535826794979 0.844327925502 -0.2 0 +0.50904141575 0.860742027004 -0.2 0 +0.481753674102 0.876306680044 -0.2 0 +0.45399049974 0.891006524188 -0.2 0 +0.425779291565 0.904827052466 -0.2 0 +0.397147890635 0.917754625684 -0.2 0 +0.368124552685 0.929776485888 -0.2 0 +0.338737920245 0.940880768954 -0.2 0 +0.309016994375 0.951056516295 -0.2 0 +0.278991106039 0.960293685677 -0.2 0 +0.248689887165 0.968583161129 -0.2 0 +0.218143241397 0.975916761939 -0.2 0 +0.187381314586 0.982287250729 -0.2 0 +0.15643446504 0.987688340595 -0.2 0 +0.125333233564 0.992114701314 -0.2 0 +0.0941083133185 0.995561964603 -0.2 0 +0.0627905195293 0.998026728428 -0.2 0 +0.0314107590781 0.999506560366 -0.2 0 +6.12323399574e-17 1.0 -0.2 2 +1.0 0.0 -0.1 1 +0.999506560366 0.0314107590781 -0.1 0 +0.998026728428 0.0627905195293 -0.1 0 +0.995561964603 0.0941083133185 -0.1 0 +0.992114701314 0.125333233564 -0.1 0 +0.987688340595 0.15643446504 -0.1 0 +0.982287250729 0.187381314586 -0.1 0 +0.975916761939 0.218143241397 -0.1 0 +0.968583161129 0.248689887165 -0.1 0 +0.960293685677 0.278991106039 -0.1 0 +0.951056516295 0.309016994375 -0.1 0 +0.940880768954 0.338737920245 -0.1 0 +0.929776485888 0.368124552685 -0.1 0 +0.917754625684 0.397147890635 -0.1 0 +0.904827052466 0.425779291565 -0.1 0 +0.891006524188 0.45399049974 -0.1 0 +0.876306680044 0.481753674102 -0.1 0 +0.860742027004 0.50904141575 -0.1 0 +0.844327925502 0.535826794979 -0.1 0 +0.827080574275 0.562083377852 -0.1 0 +0.809016994375 0.587785252292 -0.1 0 +0.790155012376 0.612907053653 -0.1 0 +0.770513242776 0.637423989749 -0.1 0 +0.75011106963 0.661311865324 -0.1 0 +0.728968627421 0.684547105929 -0.1 0 +0.707106781187 0.707106781187 -0.1 0 +0.684547105929 0.728968627421 -0.1 0 +0.661311865324 0.75011106963 -0.1 0 +0.637423989749 0.770513242776 -0.1 0 +0.612907053653 0.790155012376 -0.1 0 +0.587785252292 0.809016994375 -0.1 0 +0.562083377852 0.827080574275 -0.1 0 +0.535826794979 0.844327925502 -0.1 0 +0.50904141575 0.860742027004 -0.1 0 +0.481753674102 0.876306680044 -0.1 0 +0.45399049974 0.891006524188 -0.1 0 +0.425779291565 0.904827052466 -0.1 0 +0.397147890635 0.917754625684 -0.1 0 +0.368124552685 0.929776485888 -0.1 0 +0.338737920245 0.940880768954 -0.1 0 +0.309016994375 0.951056516295 -0.1 0 +0.278991106039 0.960293685677 -0.1 0 +0.248689887165 0.968583161129 -0.1 0 +0.218143241397 0.975916761939 -0.1 0 +0.187381314586 0.982287250729 -0.1 0 +0.15643446504 0.987688340595 -0.1 0 +0.125333233564 0.992114701314 -0.1 0 +0.0941083133185 0.995561964603 -0.1 0 +0.0627905195293 0.998026728428 -0.1 0 +0.0314107590781 0.999506560366 -0.1 0 +6.12323399574e-17 1.0 -0.1 2 +1.0 0.0 0 1 +0.999506560366 0.0314107590781 0 0 +0.998026728428 0.0627905195293 0 0 +0.995561964603 0.0941083133185 0 0 +0.992114701314 0.125333233564 0 0 +0.987688340595 0.15643446504 0 0 +0.982287250729 0.187381314586 0 0 +0.975916761939 0.218143241397 0 0 +0.968583161129 0.248689887165 0 0 +0.960293685677 0.278991106039 0 0 +0.951056516295 0.309016994375 0 0 +0.940880768954 0.338737920245 0 0 +0.929776485888 0.368124552685 0 0 +0.917754625684 0.397147890635 0 0 +0.904827052466 0.425779291565 0 0 +0.891006524188 0.45399049974 0 0 +0.876306680044 0.481753674102 0 0 +0.860742027004 0.50904141575 0 0 +0.844327925502 0.535826794979 0 0 +0.827080574275 0.562083377852 0 0 +0.809016994375 0.587785252292 0 0 +0.790155012376 0.612907053653 0 0 +0.770513242776 0.637423989749 0 0 +0.75011106963 0.661311865324 0 0 +0.728968627421 0.684547105929 0 0 +0.707106781187 0.707106781187 0 0 +0.684547105929 0.728968627421 0 0 +0.661311865324 0.75011106963 0 0 +0.637423989749 0.770513242776 0 0 +0.612907053653 0.790155012376 0 0 +0.587785252292 0.809016994375 0 0 +0.562083377852 0.827080574275 0 0 +0.535826794979 0.844327925502 0 0 +0.50904141575 0.860742027004 0 0 +0.481753674102 0.876306680044 0 0 +0.45399049974 0.891006524188 0 0 +0.425779291565 0.904827052466 0 0 +0.397147890635 0.917754625684 0 0 +0.368124552685 0.929776485888 0 0 +0.338737920245 0.940880768954 0 0 +0.309016994375 0.951056516295 0 0 +0.278991106039 0.960293685677 0 0 +0.248689887165 0.968583161129 0 0 +0.218143241397 0.975916761939 0 0 +0.187381314586 0.982287250729 0 0 +0.15643446504 0.987688340595 0 0 +0.125333233564 0.992114701314 0 0 +0.0941083133185 0.995561964603 0 0 +0.0627905195293 0.998026728428 0 0 +0.0314107590781 0.999506560366 0 0 +6.12323399574e-17 1.0 0 2 +1.0 0.0 0.1 1 +0.999506560366 0.0314107590781 0.1 0 +0.998026728428 0.0627905195293 0.1 0 +0.995561964603 0.0941083133185 0.1 0 +0.992114701314 0.125333233564 0.1 0 +0.987688340595 0.15643446504 0.1 0 +0.982287250729 0.187381314586 0.1 0 +0.975916761939 0.218143241397 0.1 0 +0.968583161129 0.248689887165 0.1 0 +0.960293685677 0.278991106039 0.1 0 +0.951056516295 0.309016994375 0.1 0 +0.940880768954 0.338737920245 0.1 0 +0.929776485888 0.368124552685 0.1 0 +0.917754625684 0.397147890635 0.1 0 +0.904827052466 0.425779291565 0.1 0 +0.891006524188 0.45399049974 0.1 0 +0.876306680044 0.481753674102 0.1 0 +0.860742027004 0.50904141575 0.1 0 +0.844327925502 0.535826794979 0.1 0 +0.827080574275 0.562083377852 0.1 0 +0.809016994375 0.587785252292 0.1 0 +0.790155012376 0.612907053653 0.1 0 +0.770513242776 0.637423989749 0.1 0 +0.75011106963 0.661311865324 0.1 0 +0.728968627421 0.684547105929 0.1 0 +0.707106781187 0.707106781187 0.1 0 +0.684547105929 0.728968627421 0.1 0 +0.661311865324 0.75011106963 0.1 0 +0.637423989749 0.770513242776 0.1 0 +0.612907053653 0.790155012376 0.1 0 +0.587785252292 0.809016994375 0.1 0 +0.562083377852 0.827080574275 0.1 0 +0.535826794979 0.844327925502 0.1 0 +0.50904141575 0.860742027004 0.1 0 +0.481753674102 0.876306680044 0.1 0 +0.45399049974 0.891006524188 0.1 0 +0.425779291565 0.904827052466 0.1 0 +0.397147890635 0.917754625684 0.1 0 +0.368124552685 0.929776485888 0.1 0 +0.338737920245 0.940880768954 0.1 0 +0.309016994375 0.951056516295 0.1 0 +0.278991106039 0.960293685677 0.1 0 +0.248689887165 0.968583161129 0.1 0 +0.218143241397 0.975916761939 0.1 0 +0.187381314586 0.982287250729 0.1 0 +0.15643446504 0.987688340595 0.1 0 +0.125333233564 0.992114701314 0.1 0 +0.0941083133185 0.995561964603 0.1 0 +0.0627905195293 0.998026728428 0.1 0 +0.0314107590781 0.999506560366 0.1 0 +6.12323399574e-17 1.0 0.1 2 +1.0 0.0 0.2 1 +0.999506560366 0.0314107590781 0.2 0 +0.998026728428 0.0627905195293 0.2 0 +0.995561964603 0.0941083133185 0.2 0 +0.992114701314 0.125333233564 0.2 0 +0.987688340595 0.15643446504 0.2 0 +0.982287250729 0.187381314586 0.2 0 +0.975916761939 0.218143241397 0.2 0 +0.968583161129 0.248689887165 0.2 0 +0.960293685677 0.278991106039 0.2 0 +0.951056516295 0.309016994375 0.2 0 +0.940880768954 0.338737920245 0.2 0 +0.929776485888 0.368124552685 0.2 0 +0.917754625684 0.397147890635 0.2 0 +0.904827052466 0.425779291565 0.2 0 +0.891006524188 0.45399049974 0.2 0 +0.876306680044 0.481753674102 0.2 0 +0.860742027004 0.50904141575 0.2 0 +0.844327925502 0.535826794979 0.2 0 +0.827080574275 0.562083377852 0.2 0 +0.809016994375 0.587785252292 0.2 0 +0.790155012376 0.612907053653 0.2 0 +0.770513242776 0.637423989749 0.2 0 +0.75011106963 0.661311865324 0.2 0 +0.728968627421 0.684547105929 0.2 0 +0.707106781187 0.707106781187 0.2 0 +0.684547105929 0.728968627421 0.2 0 +0.661311865324 0.75011106963 0.2 0 +0.637423989749 0.770513242776 0.2 0 +0.612907053653 0.790155012376 0.2 0 +0.587785252292 0.809016994375 0.2 0 +0.562083377852 0.827080574275 0.2 0 +0.535826794979 0.844327925502 0.2 0 +0.50904141575 0.860742027004 0.2 0 +0.481753674102 0.876306680044 0.2 0 +0.45399049974 0.891006524188 0.2 0 +0.425779291565 0.904827052466 0.2 0 +0.397147890635 0.917754625684 0.2 0 +0.368124552685 0.929776485888 0.2 0 +0.338737920245 0.940880768954 0.2 0 +0.309016994375 0.951056516295 0.2 0 +0.278991106039 0.960293685677 0.2 0 +0.248689887165 0.968583161129 0.2 0 +0.218143241397 0.975916761939 0.2 0 +0.187381314586 0.982287250729 0.2 0 +0.15643446504 0.987688340595 0.2 0 +0.125333233564 0.992114701314 0.2 0 +0.0941083133185 0.995561964603 0.2 0 +0.0627905195293 0.998026728428 0.2 0 +0.0314107590781 0.999506560366 0.2 0 +6.12323399574e-17 1.0 0.2 2 +1.0 0.0 0.3 1 +0.999506560366 0.0314107590781 0.3 0 +0.998026728428 0.0627905195293 0.3 0 +0.995561964603 0.0941083133185 0.3 0 +0.992114701314 0.125333233564 0.3 0 +0.987688340595 0.15643446504 0.3 0 +0.982287250729 0.187381314586 0.3 0 +0.975916761939 0.218143241397 0.3 0 +0.968583161129 0.248689887165 0.3 0 +0.960293685677 0.278991106039 0.3 0 +0.951056516295 0.309016994375 0.3 0 +0.940880768954 0.338737920245 0.3 0 +0.929776485888 0.368124552685 0.3 0 +0.917754625684 0.397147890635 0.3 0 +0.904827052466 0.425779291565 0.3 0 +0.891006524188 0.45399049974 0.3 0 +0.876306680044 0.481753674102 0.3 0 +0.860742027004 0.50904141575 0.3 0 +0.844327925502 0.535826794979 0.3 0 +0.827080574275 0.562083377852 0.3 0 +0.809016994375 0.587785252292 0.3 0 +0.790155012376 0.612907053653 0.3 0 +0.770513242776 0.637423989749 0.3 0 +0.75011106963 0.661311865324 0.3 0 +0.728968627421 0.684547105929 0.3 0 +0.707106781187 0.707106781187 0.3 0 +0.684547105929 0.728968627421 0.3 0 +0.661311865324 0.75011106963 0.3 0 +0.637423989749 0.770513242776 0.3 0 +0.612907053653 0.790155012376 0.3 0 +0.587785252292 0.809016994375 0.3 0 +0.562083377852 0.827080574275 0.3 0 +0.535826794979 0.844327925502 0.3 0 +0.50904141575 0.860742027004 0.3 0 +0.481753674102 0.876306680044 0.3 0 +0.45399049974 0.891006524188 0.3 0 +0.425779291565 0.904827052466 0.3 0 +0.397147890635 0.917754625684 0.3 0 +0.368124552685 0.929776485888 0.3 0 +0.338737920245 0.940880768954 0.3 0 +0.309016994375 0.951056516295 0.3 0 +0.278991106039 0.960293685677 0.3 0 +0.248689887165 0.968583161129 0.3 0 +0.218143241397 0.975916761939 0.3 0 +0.187381314586 0.982287250729 0.3 0 +0.15643446504 0.987688340595 0.3 0 +0.125333233564 0.992114701314 0.3 0 +0.0941083133185 0.995561964603 0.3 0 +0.0627905195293 0.998026728428 0.3 0 +0.0314107590781 0.999506560366 0.3 0 +6.12323399574e-17 1.0 0.3 2 +1.0 0.0 0.4 1 +0.999506560366 0.0314107590781 0.4 0 +0.998026728428 0.0627905195293 0.4 0 +0.995561964603 0.0941083133185 0.4 0 +0.992114701314 0.125333233564 0.4 0 +0.987688340595 0.15643446504 0.4 0 +0.982287250729 0.187381314586 0.4 0 +0.975916761939 0.218143241397 0.4 0 +0.968583161129 0.248689887165 0.4 0 +0.960293685677 0.278991106039 0.4 0 +0.951056516295 0.309016994375 0.4 0 +0.940880768954 0.338737920245 0.4 0 +0.929776485888 0.368124552685 0.4 0 +0.917754625684 0.397147890635 0.4 0 +0.904827052466 0.425779291565 0.4 0 +0.891006524188 0.45399049974 0.4 0 +0.876306680044 0.481753674102 0.4 0 +0.860742027004 0.50904141575 0.4 0 +0.844327925502 0.535826794979 0.4 0 +0.827080574275 0.562083377852 0.4 0 +0.809016994375 0.587785252292 0.4 0 +0.790155012376 0.612907053653 0.4 0 +0.770513242776 0.637423989749 0.4 0 +0.75011106963 0.661311865324 0.4 0 +0.728968627421 0.684547105929 0.4 0 +0.707106781187 0.707106781187 0.4 0 +0.684547105929 0.728968627421 0.4 0 +0.661311865324 0.75011106963 0.4 0 +0.637423989749 0.770513242776 0.4 0 +0.612907053653 0.790155012376 0.4 0 +0.587785252292 0.809016994375 0.4 0 +0.562083377852 0.827080574275 0.4 0 +0.535826794979 0.844327925502 0.4 0 +0.50904141575 0.860742027004 0.4 0 +0.481753674102 0.876306680044 0.4 0 +0.45399049974 0.891006524188 0.4 0 +0.425779291565 0.904827052466 0.4 0 +0.397147890635 0.917754625684 0.4 0 +0.368124552685 0.929776485888 0.4 0 +0.338737920245 0.940880768954 0.4 0 +0.309016994375 0.951056516295 0.4 0 +0.278991106039 0.960293685677 0.4 0 +0.248689887165 0.968583161129 0.4 0 +0.218143241397 0.975916761939 0.4 0 +0.187381314586 0.982287250729 0.4 0 +0.15643446504 0.987688340595 0.4 0 +0.125333233564 0.992114701314 0.4 0 +0.0941083133185 0.995561964603 0.4 0 +0.0627905195293 0.998026728428 0.4 0 +0.0314107590781 0.999506560366 0.4 0 +6.12323399574e-17 1.0 0.4 2 +1.0 0.0 0.5 1 +0.999506560366 0.0314107590781 0.5 0 +0.998026728428 0.0627905195293 0.5 0 +0.995561964603 0.0941083133185 0.5 0 +0.992114701314 0.125333233564 0.5 0 +0.987688340595 0.15643446504 0.5 0 +0.982287250729 0.187381314586 0.5 0 +0.975916761939 0.218143241397 0.5 0 +0.968583161129 0.248689887165 0.5 0 +0.960293685677 0.278991106039 0.5 0 +0.951056516295 0.309016994375 0.5 0 +0.940880768954 0.338737920245 0.5 0 +0.929776485888 0.368124552685 0.5 0 +0.917754625684 0.397147890635 0.5 0 +0.904827052466 0.425779291565 0.5 0 +0.891006524188 0.45399049974 0.5 0 +0.876306680044 0.481753674102 0.5 0 +0.860742027004 0.50904141575 0.5 0 +0.844327925502 0.535826794979 0.5 0 +0.827080574275 0.562083377852 0.5 0 +0.809016994375 0.587785252292 0.5 0 +0.790155012376 0.612907053653 0.5 0 +0.770513242776 0.637423989749 0.5 0 +0.75011106963 0.661311865324 0.5 0 +0.728968627421 0.684547105929 0.5 0 +0.707106781187 0.707106781187 0.5 0 +0.684547105929 0.728968627421 0.5 0 +0.661311865324 0.75011106963 0.5 0 +0.637423989749 0.770513242776 0.5 0 +0.612907053653 0.790155012376 0.5 0 +0.587785252292 0.809016994375 0.5 0 +0.562083377852 0.827080574275 0.5 0 +0.535826794979 0.844327925502 0.5 0 +0.50904141575 0.860742027004 0.5 0 +0.481753674102 0.876306680044 0.5 0 +0.45399049974 0.891006524188 0.5 0 +0.425779291565 0.904827052466 0.5 0 +0.397147890635 0.917754625684 0.5 0 +0.368124552685 0.929776485888 0.5 0 +0.338737920245 0.940880768954 0.5 0 +0.309016994375 0.951056516295 0.5 0 +0.278991106039 0.960293685677 0.5 0 +0.248689887165 0.968583161129 0.5 0 +0.218143241397 0.975916761939 0.5 0 +0.187381314586 0.982287250729 0.5 0 +0.15643446504 0.987688340595 0.5 0 +0.125333233564 0.992114701314 0.5 0 +0.0941083133185 0.995561964603 0.5 0 +0.0627905195293 0.998026728428 0.5 0 +0.0314107590781 0.999506560366 0.5 0 +6.12323399574e-17 1.0 0.5 2 +1.0 0.0 0.6 1 +0.999506560366 0.0314107590781 0.6 0 +0.998026728428 0.0627905195293 0.6 0 +0.995561964603 0.0941083133185 0.6 0 +0.992114701314 0.125333233564 0.6 0 +0.987688340595 0.15643446504 0.6 0 +0.982287250729 0.187381314586 0.6 0 +0.975916761939 0.218143241397 0.6 0 +0.968583161129 0.248689887165 0.6 0 +0.960293685677 0.278991106039 0.6 0 +0.951056516295 0.309016994375 0.6 0 +0.940880768954 0.338737920245 0.6 0 +0.929776485888 0.368124552685 0.6 0 +0.917754625684 0.397147890635 0.6 0 +0.904827052466 0.425779291565 0.6 0 +0.891006524188 0.45399049974 0.6 0 +0.876306680044 0.481753674102 0.6 0 +0.860742027004 0.50904141575 0.6 0 +0.844327925502 0.535826794979 0.6 0 +0.827080574275 0.562083377852 0.6 0 +0.809016994375 0.587785252292 0.6 0 +0.790155012376 0.612907053653 0.6 0 +0.770513242776 0.637423989749 0.6 0 +0.75011106963 0.661311865324 0.6 0 +0.728968627421 0.684547105929 0.6 0 +0.707106781187 0.707106781187 0.6 0 +0.684547105929 0.728968627421 0.6 0 +0.661311865324 0.75011106963 0.6 0 +0.637423989749 0.770513242776 0.6 0 +0.612907053653 0.790155012376 0.6 0 +0.587785252292 0.809016994375 0.6 0 +0.562083377852 0.827080574275 0.6 0 +0.535826794979 0.844327925502 0.6 0 +0.50904141575 0.860742027004 0.6 0 +0.481753674102 0.876306680044 0.6 0 +0.45399049974 0.891006524188 0.6 0 +0.425779291565 0.904827052466 0.6 0 +0.397147890635 0.917754625684 0.6 0 +0.368124552685 0.929776485888 0.6 0 +0.338737920245 0.940880768954 0.6 0 +0.309016994375 0.951056516295 0.6 0 +0.278991106039 0.960293685677 0.6 0 +0.248689887165 0.968583161129 0.6 0 +0.218143241397 0.975916761939 0.6 0 +0.187381314586 0.982287250729 0.6 0 +0.15643446504 0.987688340595 0.6 0 +0.125333233564 0.992114701314 0.6 0 +0.0941083133185 0.995561964603 0.6 0 +0.0627905195293 0.998026728428 0.6 0 +0.0314107590781 0.999506560366 0.6 0 +6.12323399574e-17 1.0 0.6 2 +1.0 0.0 0.7 1 +0.999506560366 0.0314107590781 0.7 0 +0.998026728428 0.0627905195293 0.7 0 +0.995561964603 0.0941083133185 0.7 0 +0.992114701314 0.125333233564 0.7 0 +0.987688340595 0.15643446504 0.7 0 +0.982287250729 0.187381314586 0.7 0 +0.975916761939 0.218143241397 0.7 0 +0.968583161129 0.248689887165 0.7 0 +0.960293685677 0.278991106039 0.7 0 +0.951056516295 0.309016994375 0.7 0 +0.940880768954 0.338737920245 0.7 0 +0.929776485888 0.368124552685 0.7 0 +0.917754625684 0.397147890635 0.7 0 +0.904827052466 0.425779291565 0.7 0 +0.891006524188 0.45399049974 0.7 0 +0.876306680044 0.481753674102 0.7 0 +0.860742027004 0.50904141575 0.7 0 +0.844327925502 0.535826794979 0.7 0 +0.827080574275 0.562083377852 0.7 0 +0.809016994375 0.587785252292 0.7 0 +0.790155012376 0.612907053653 0.7 0 +0.770513242776 0.637423989749 0.7 0 +0.75011106963 0.661311865324 0.7 0 +0.728968627421 0.684547105929 0.7 0 +0.707106781187 0.707106781187 0.7 0 +0.684547105929 0.728968627421 0.7 0 +0.661311865324 0.75011106963 0.7 0 +0.637423989749 0.770513242776 0.7 0 +0.612907053653 0.790155012376 0.7 0 +0.587785252292 0.809016994375 0.7 0 +0.562083377852 0.827080574275 0.7 0 +0.535826794979 0.844327925502 0.7 0 +0.50904141575 0.860742027004 0.7 0 +0.481753674102 0.876306680044 0.7 0 +0.45399049974 0.891006524188 0.7 0 +0.425779291565 0.904827052466 0.7 0 +0.397147890635 0.917754625684 0.7 0 +0.368124552685 0.929776485888 0.7 0 +0.338737920245 0.940880768954 0.7 0 +0.309016994375 0.951056516295 0.7 0 +0.278991106039 0.960293685677 0.7 0 +0.248689887165 0.968583161129 0.7 0 +0.218143241397 0.975916761939 0.7 0 +0.187381314586 0.982287250729 0.7 0 +0.15643446504 0.987688340595 0.7 0 +0.125333233564 0.992114701314 0.7 0 +0.0941083133185 0.995561964603 0.7 0 +0.0627905195293 0.998026728428 0.7 0 +0.0314107590781 0.999506560366 0.7 0 +6.12323399574e-17 1.0 0.7 2 +1.0 0.0 0.8 1 +0.999506560366 0.0314107590781 0.8 0 +0.998026728428 0.0627905195293 0.8 0 +0.995561964603 0.0941083133185 0.8 0 +0.992114701314 0.125333233564 0.8 0 +0.987688340595 0.15643446504 0.8 0 +0.982287250729 0.187381314586 0.8 0 +0.975916761939 0.218143241397 0.8 0 +0.968583161129 0.248689887165 0.8 0 +0.960293685677 0.278991106039 0.8 0 +0.951056516295 0.309016994375 0.8 0 +0.940880768954 0.338737920245 0.8 0 +0.929776485888 0.368124552685 0.8 0 +0.917754625684 0.397147890635 0.8 0 +0.904827052466 0.425779291565 0.8 0 +0.891006524188 0.45399049974 0.8 0 +0.876306680044 0.481753674102 0.8 0 +0.860742027004 0.50904141575 0.8 0 +0.844327925502 0.535826794979 0.8 0 +0.827080574275 0.562083377852 0.8 0 +0.809016994375 0.587785252292 0.8 0 +0.790155012376 0.612907053653 0.8 0 +0.770513242776 0.637423989749 0.8 0 +0.75011106963 0.661311865324 0.8 0 +0.728968627421 0.684547105929 0.8 0 +0.707106781187 0.707106781187 0.8 0 +0.684547105929 0.728968627421 0.8 0 +0.661311865324 0.75011106963 0.8 0 +0.637423989749 0.770513242776 0.8 0 +0.612907053653 0.790155012376 0.8 0 +0.587785252292 0.809016994375 0.8 0 +0.562083377852 0.827080574275 0.8 0 +0.535826794979 0.844327925502 0.8 0 +0.50904141575 0.860742027004 0.8 0 +0.481753674102 0.876306680044 0.8 0 +0.45399049974 0.891006524188 0.8 0 +0.425779291565 0.904827052466 0.8 0 +0.397147890635 0.917754625684 0.8 0 +0.368124552685 0.929776485888 0.8 0 +0.338737920245 0.940880768954 0.8 0 +0.309016994375 0.951056516295 0.8 0 +0.278991106039 0.960293685677 0.8 0 +0.248689887165 0.968583161129 0.8 0 +0.218143241397 0.975916761939 0.8 0 +0.187381314586 0.982287250729 0.8 0 +0.15643446504 0.987688340595 0.8 0 +0.125333233564 0.992114701314 0.8 0 +0.0941083133185 0.995561964603 0.8 0 +0.0627905195293 0.998026728428 0.8 0 +0.0314107590781 0.999506560366 0.8 0 +6.12323399574e-17 1.0 0.8 2 +1.0 0.0 0.9 1 +0.999506560366 0.0314107590781 0.9 0 +0.998026728428 0.0627905195293 0.9 0 +0.995561964603 0.0941083133185 0.9 0 +0.992114701314 0.125333233564 0.9 0 +0.987688340595 0.15643446504 0.9 0 +0.982287250729 0.187381314586 0.9 0 +0.975916761939 0.218143241397 0.9 0 +0.968583161129 0.248689887165 0.9 0 +0.960293685677 0.278991106039 0.9 0 +0.951056516295 0.309016994375 0.9 0 +0.940880768954 0.338737920245 0.9 0 +0.929776485888 0.368124552685 0.9 0 +0.917754625684 0.397147890635 0.9 0 +0.904827052466 0.425779291565 0.9 0 +0.891006524188 0.45399049974 0.9 0 +0.876306680044 0.481753674102 0.9 0 +0.860742027004 0.50904141575 0.9 0 +0.844327925502 0.535826794979 0.9 0 +0.827080574275 0.562083377852 0.9 0 +0.809016994375 0.587785252292 0.9 0 +0.790155012376 0.612907053653 0.9 0 +0.770513242776 0.637423989749 0.9 0 +0.75011106963 0.661311865324 0.9 0 +0.728968627421 0.684547105929 0.9 0 +0.707106781187 0.707106781187 0.9 0 +0.684547105929 0.728968627421 0.9 0 +0.661311865324 0.75011106963 0.9 0 +0.637423989749 0.770513242776 0.9 0 +0.612907053653 0.790155012376 0.9 0 +0.587785252292 0.809016994375 0.9 0 +0.562083377852 0.827080574275 0.9 0 +0.535826794979 0.844327925502 0.9 0 +0.50904141575 0.860742027004 0.9 0 +0.481753674102 0.876306680044 0.9 0 +0.45399049974 0.891006524188 0.9 0 +0.425779291565 0.904827052466 0.9 0 +0.397147890635 0.917754625684 0.9 0 +0.368124552685 0.929776485888 0.9 0 +0.338737920245 0.940880768954 0.9 0 +0.309016994375 0.951056516295 0.9 0 +0.278991106039 0.960293685677 0.9 0 +0.248689887165 0.968583161129 0.9 0 +0.218143241397 0.975916761939 0.9 0 +0.187381314586 0.982287250729 0.9 0 +0.15643446504 0.987688340595 0.9 0 +0.125333233564 0.992114701314 0.9 0 +0.0941083133185 0.995561964603 0.9 0 +0.0627905195293 0.998026728428 0.9 0 +0.0314107590781 0.999506560366 0.9 0 +6.12323399574e-17 1.0 0.9 2 +1.0 0.0 1 1 +0.999506560366 0.0314107590781 1 0 +0.998026728428 0.0627905195293 1 0 +0.995561964603 0.0941083133185 1 0 +0.992114701314 0.125333233564 1 0 +0.987688340595 0.15643446504 1 0 +0.982287250729 0.187381314586 1 0 +0.975916761939 0.218143241397 1 0 +0.968583161129 0.248689887165 1 0 +0.960293685677 0.278991106039 1 0 +0.951056516295 0.309016994375 1 0 +0.940880768954 0.338737920245 1 0 +0.929776485888 0.368124552685 1 0 +0.917754625684 0.397147890635 1 0 +0.904827052466 0.425779291565 1 0 +0.891006524188 0.45399049974 1 0 +0.876306680044 0.481753674102 1 0 +0.860742027004 0.50904141575 1 0 +0.844327925502 0.535826794979 1 0 +0.827080574275 0.562083377852 1 0 +0.809016994375 0.587785252292 1 0 +0.790155012376 0.612907053653 1 0 +0.770513242776 0.637423989749 1 0 +0.75011106963 0.661311865324 1 0 +0.728968627421 0.684547105929 1 0 +0.707106781187 0.707106781187 1 0 +0.684547105929 0.728968627421 1 0 +0.661311865324 0.75011106963 1 0 +0.637423989749 0.770513242776 1 0 +0.612907053653 0.790155012376 1 0 +0.587785252292 0.809016994375 1 0 +0.562083377852 0.827080574275 1 0 +0.535826794979 0.844327925502 1 0 +0.50904141575 0.860742027004 1 0 +0.481753674102 0.876306680044 1 0 +0.45399049974 0.891006524188 1 0 +0.425779291565 0.904827052466 1 0 +0.397147890635 0.917754625684 1 0 +0.368124552685 0.929776485888 1 0 +0.338737920245 0.940880768954 1 0 +0.309016994375 0.951056516295 1 0 +0.278991106039 0.960293685677 1 0 +0.248689887165 0.968583161129 1 0 +0.218143241397 0.975916761939 1 0 +0.187381314586 0.982287250729 1 0 +0.15643446504 0.987688340595 1 0 +0.125333233564 0.992114701314 1 0 +0.0941083133185 0.995561964603 1 0 +0.0627905195293 0.998026728428 1 0 +0.0314107590781 0.999506560366 1 0 +6.12323399574e-17 1.0 1 0 + +Edges +140 +1 2 3 +2 3 3 +3 4 3 +4 5 3 +5 6 3 +6 7 3 +7 8 3 +8 9 3 +9 10 3 +10 11 3 +11 12 3 +12 13 3 +13 14 3 +14 15 3 +15 16 3 +16 17 3 +17 18 3 +18 19 3 +19 20 3 +20 21 3 +21 22 3 +22 23 3 +23 24 3 +24 25 3 +25 26 3 +26 27 3 +27 28 3 +28 29 3 +29 30 3 +30 31 3 +31 32 3 +32 33 3 +33 34 3 +34 35 3 +35 36 3 +36 37 3 +37 38 3 +38 39 3 +39 40 3 +40 41 3 +41 42 3 +42 43 3 +43 44 3 +44 45 3 +45 46 3 +46 47 3 +47 48 3 +48 49 3 +49 50 3 +50 51 3 +51 102 2 +102 153 2 +153 204 2 +204 255 2 +255 306 2 +306 357 2 +357 408 2 +408 459 2 +459 510 2 +510 561 2 +561 612 2 +612 663 2 +663 714 2 +714 765 2 +765 816 2 +816 867 2 +867 918 2 +918 969 2 +969 1020 2 +1020 1071 2 +1022 1021 4 +1023 1022 4 +1024 1023 4 +1025 1024 4 +1026 1025 4 +1027 1026 4 +1028 1027 4 +1029 1028 4 +1030 1029 4 +1031 1030 4 +1032 1031 4 +1033 1032 4 +1034 1033 4 +1035 1034 4 +1036 1035 4 +1037 1036 4 +1038 1037 4 +1039 1038 4 +1040 1039 4 +1041 1040 4 +1042 1041 4 +1043 1042 4 +1044 1043 4 +1045 1044 4 +1046 1045 4 +1047 1046 4 +1048 1047 4 +1049 1048 4 +1050 1049 4 +1051 1050 4 +1052 1051 4 +1053 1052 4 +1054 1053 4 +1055 1054 4 +1056 1055 4 +1057 1056 4 +1058 1057 4 +1059 1058 4 +1060 1059 4 +1061 1060 4 +1062 1061 4 +1063 1062 4 +1064 1063 4 +1065 1064 4 +1066 1065 4 +1067 1066 4 +1068 1067 4 +1069 1068 4 +1070 1069 4 +1071 1070 4 +52 1 1 +103 52 1 +154 103 1 +205 154 1 +256 205 1 +307 256 1 +358 307 1 +409 358 1 +460 409 1 +511 460 1 +562 511 1 +613 562 1 +664 613 1 +715 664 1 +766 715 1 +817 766 1 +868 817 1 +919 868 1 +970 919 1 +1021 970 1 + +Triangles +2000 +1 2 53 0 +1 53 52 0 +2 3 54 0 +2 54 53 0 +3 4 55 0 +3 55 54 0 +4 5 56 0 +4 56 55 0 +5 6 57 0 +5 57 56 0 +6 7 58 0 +6 58 57 0 +7 8 59 0 +7 59 58 0 +8 9 60 0 +8 60 59 0 +9 10 61 0 +9 61 60 0 +10 11 62 0 +10 62 61 0 +11 12 63 0 +11 63 62 0 +12 13 64 0 +12 64 63 0 +13 14 65 0 +13 65 64 0 +14 15 66 0 +14 66 65 0 +15 16 67 0 +15 67 66 0 +16 17 68 0 +16 68 67 0 +17 18 69 0 +17 69 68 0 +18 19 70 0 +18 70 69 0 +19 20 71 0 +19 71 70 0 +20 21 72 0 +20 72 71 0 +21 22 73 0 +21 73 72 0 +22 23 74 0 +22 74 73 0 +23 24 75 0 +23 75 74 0 +24 25 76 0 +24 76 75 0 +25 26 77 0 +25 77 76 0 +26 27 78 0 +26 78 77 0 +27 28 79 0 +27 79 78 0 +28 29 80 0 +28 80 79 0 +29 30 81 0 +29 81 80 0 +30 31 82 0 +30 82 81 0 +31 32 83 0 +31 83 82 0 +32 33 84 0 +32 84 83 0 +33 34 85 0 +33 85 84 0 +34 35 86 0 +34 86 85 0 +35 36 87 0 +35 87 86 0 +36 37 88 0 +36 88 87 0 +37 38 89 0 +37 89 88 0 +38 39 90 0 +38 90 89 0 +39 40 91 0 +39 91 90 0 +40 41 92 0 +40 92 91 0 +41 42 93 0 +41 93 92 0 +42 43 94 0 +42 94 93 0 +43 44 95 0 +43 95 94 0 +44 45 96 0 +44 96 95 0 +45 46 97 0 +45 97 96 0 +46 47 98 0 +46 98 97 0 +47 48 99 0 +47 99 98 0 +48 49 100 0 +48 100 99 0 +49 50 101 0 +49 101 100 0 +50 51 102 0 +50 102 101 0 +52 53 104 0 +52 104 103 0 +53 54 105 0 +53 105 104 0 +54 55 106 0 +54 106 105 0 +55 56 107 0 +55 107 106 0 +56 57 108 0 +56 108 107 0 +57 58 109 0 +57 109 108 0 +58 59 110 0 +58 110 109 0 +59 60 111 0 +59 111 110 0 +60 61 112 0 +60 112 111 0 +61 62 113 0 +61 113 112 0 +62 63 114 0 +62 114 113 0 +63 64 115 0 +63 115 114 0 +64 65 116 0 +64 116 115 0 +65 66 117 0 +65 117 116 0 +66 67 118 0 +66 118 117 0 +67 68 119 0 +67 119 118 0 +68 69 120 0 +68 120 119 0 +69 70 121 0 +69 121 120 0 +70 71 122 0 +70 122 121 0 +71 72 123 0 +71 123 122 0 +72 73 124 0 +72 124 123 0 +73 74 125 0 +73 125 124 0 +74 75 126 0 +74 126 125 0 +75 76 127 0 +75 127 126 0 +76 77 128 0 +76 128 127 0 +77 78 129 0 +77 129 128 0 +78 79 130 0 +78 130 129 0 +79 80 131 0 +79 131 130 0 +80 81 132 0 +80 132 131 0 +81 82 133 0 +81 133 132 0 +82 83 134 0 +82 134 133 0 +83 84 135 0 +83 135 134 0 +84 85 136 0 +84 136 135 0 +85 86 137 0 +85 137 136 0 +86 87 138 0 +86 138 137 0 +87 88 139 0 +87 139 138 0 +88 89 140 0 +88 140 139 0 +89 90 141 0 +89 141 140 0 +90 91 142 0 +90 142 141 0 +91 92 143 0 +91 143 142 0 +92 93 144 0 +92 144 143 0 +93 94 145 0 +93 145 144 0 +94 95 146 0 +94 146 145 0 +95 96 147 0 +95 147 146 0 +96 97 148 0 +96 148 147 0 +97 98 149 0 +97 149 148 0 +98 99 150 0 +98 150 149 0 +99 100 151 0 +99 151 150 0 +100 101 152 0 +100 152 151 0 +101 102 153 0 +101 153 152 0 +103 104 155 0 +103 155 154 0 +104 105 156 0 +104 156 155 0 +105 106 157 0 +105 157 156 0 +106 107 158 0 +106 158 157 0 +107 108 159 0 +107 159 158 0 +108 109 160 0 +108 160 159 0 +109 110 161 0 +109 161 160 0 +110 111 162 0 +110 162 161 0 +111 112 163 0 +111 163 162 0 +112 113 164 0 +112 164 163 0 +113 114 165 0 +113 165 164 0 +114 115 166 0 +114 166 165 0 +115 116 167 0 +115 167 166 0 +116 117 168 0 +116 168 167 0 +117 118 169 0 +117 169 168 0 +118 119 170 0 +118 170 169 0 +119 120 171 0 +119 171 170 0 +120 121 172 0 +120 172 171 0 +121 122 173 0 +121 173 172 0 +122 123 174 0 +122 174 173 0 +123 124 175 0 +123 175 174 0 +124 125 176 0 +124 176 175 0 +125 126 177 0 +125 177 176 0 +126 127 178 0 +126 178 177 0 +127 128 179 0 +127 179 178 0 +128 129 180 0 +128 180 179 0 +129 130 181 0 +129 181 180 0 +130 131 182 0 +130 182 181 0 +131 132 183 0 +131 183 182 0 +132 133 184 0 +132 184 183 0 +133 134 185 0 +133 185 184 0 +134 135 186 0 +134 186 185 0 +135 136 187 0 +135 187 186 0 +136 137 188 0 +136 188 187 0 +137 138 189 0 +137 189 188 0 +138 139 190 0 +138 190 189 0 +139 140 191 0 +139 191 190 0 +140 141 192 0 +140 192 191 0 +141 142 193 0 +141 193 192 0 +142 143 194 0 +142 194 193 0 +143 144 195 0 +143 195 194 0 +144 145 196 0 +144 196 195 0 +145 146 197 0 +145 197 196 0 +146 147 198 0 +146 198 197 0 +147 148 199 0 +147 199 198 0 +148 149 200 0 +148 200 199 0 +149 150 201 0 +149 201 200 0 +150 151 202 0 +150 202 201 0 +151 152 203 0 +151 203 202 0 +152 153 204 0 +152 204 203 0 +154 155 206 0 +154 206 205 0 +155 156 207 0 +155 207 206 0 +156 157 208 0 +156 208 207 0 +157 158 209 0 +157 209 208 0 +158 159 210 0 +158 210 209 0 +159 160 211 0 +159 211 210 0 +160 161 212 0 +160 212 211 0 +161 162 213 0 +161 213 212 0 +162 163 214 0 +162 214 213 0 +163 164 215 0 +163 215 214 0 +164 165 216 0 +164 216 215 0 +165 166 217 0 +165 217 216 0 +166 167 218 0 +166 218 217 0 +167 168 219 0 +167 219 218 0 +168 169 220 0 +168 220 219 0 +169 170 221 0 +169 221 220 0 +170 171 222 0 +170 222 221 0 +171 172 223 0 +171 223 222 0 +172 173 224 0 +172 224 223 0 +173 174 225 0 +173 225 224 0 +174 175 226 0 +174 226 225 0 +175 176 227 0 +175 227 226 0 +176 177 228 0 +176 228 227 0 +177 178 229 0 +177 229 228 0 +178 179 230 0 +178 230 229 0 +179 180 231 0 +179 231 230 0 +180 181 232 0 +180 232 231 0 +181 182 233 0 +181 233 232 0 +182 183 234 0 +182 234 233 0 +183 184 235 0 +183 235 234 0 +184 185 236 0 +184 236 235 0 +185 186 237 0 +185 237 236 0 +186 187 238 0 +186 238 237 0 +187 188 239 0 +187 239 238 0 +188 189 240 0 +188 240 239 0 +189 190 241 0 +189 241 240 0 +190 191 242 0 +190 242 241 0 +191 192 243 0 +191 243 242 0 +192 193 244 0 +192 244 243 0 +193 194 245 0 +193 245 244 0 +194 195 246 0 +194 246 245 0 +195 196 247 0 +195 247 246 0 +196 197 248 0 +196 248 247 0 +197 198 249 0 +197 249 248 0 +198 199 250 0 +198 250 249 0 +199 200 251 0 +199 251 250 0 +200 201 252 0 +200 252 251 0 +201 202 253 0 +201 253 252 0 +202 203 254 0 +202 254 253 0 +203 204 255 0 +203 255 254 0 +205 206 257 0 +205 257 256 0 +206 207 258 0 +206 258 257 0 +207 208 259 0 +207 259 258 0 +208 209 260 0 +208 260 259 0 +209 210 261 0 +209 261 260 0 +210 211 262 0 +210 262 261 0 +211 212 263 0 +211 263 262 0 +212 213 264 0 +212 264 263 0 +213 214 265 0 +213 265 264 0 +214 215 266 0 +214 266 265 0 +215 216 267 0 +215 267 266 0 +216 217 268 0 +216 268 267 0 +217 218 269 0 +217 269 268 0 +218 219 270 0 +218 270 269 0 +219 220 271 0 +219 271 270 0 +220 221 272 0 +220 272 271 0 +221 222 273 0 +221 273 272 0 +222 223 274 0 +222 274 273 0 +223 224 275 0 +223 275 274 0 +224 225 276 0 +224 276 275 0 +225 226 277 0 +225 277 276 0 +226 227 278 0 +226 278 277 0 +227 228 279 0 +227 279 278 0 +228 229 280 0 +228 280 279 0 +229 230 281 0 +229 281 280 0 +230 231 282 0 +230 282 281 0 +231 232 283 0 +231 283 282 0 +232 233 284 0 +232 284 283 0 +233 234 285 0 +233 285 284 0 +234 235 286 0 +234 286 285 0 +235 236 287 0 +235 287 286 0 +236 237 288 0 +236 288 287 0 +237 238 289 0 +237 289 288 0 +238 239 290 0 +238 290 289 0 +239 240 291 0 +239 291 290 0 +240 241 292 0 +240 292 291 0 +241 242 293 0 +241 293 292 0 +242 243 294 0 +242 294 293 0 +243 244 295 0 +243 295 294 0 +244 245 296 0 +244 296 295 0 +245 246 297 0 +245 297 296 0 +246 247 298 0 +246 298 297 0 +247 248 299 0 +247 299 298 0 +248 249 300 0 +248 300 299 0 +249 250 301 0 +249 301 300 0 +250 251 302 0 +250 302 301 0 +251 252 303 0 +251 303 302 0 +252 253 304 0 +252 304 303 0 +253 254 305 0 +253 305 304 0 +254 255 306 0 +254 306 305 0 +256 257 308 0 +256 308 307 0 +257 258 309 0 +257 309 308 0 +258 259 310 0 +258 310 309 0 +259 260 311 0 +259 311 310 0 +260 261 312 0 +260 312 311 0 +261 262 313 0 +261 313 312 0 +262 263 314 0 +262 314 313 0 +263 264 315 0 +263 315 314 0 +264 265 316 0 +264 316 315 0 +265 266 317 0 +265 317 316 0 +266 267 318 0 +266 318 317 0 +267 268 319 0 +267 319 318 0 +268 269 320 0 +268 320 319 0 +269 270 321 0 +269 321 320 0 +270 271 322 0 +270 322 321 0 +271 272 323 0 +271 323 322 0 +272 273 324 0 +272 324 323 0 +273 274 325 0 +273 325 324 0 +274 275 326 0 +274 326 325 0 +275 276 327 0 +275 327 326 0 +276 277 328 0 +276 328 327 0 +277 278 329 0 +277 329 328 0 +278 279 330 0 +278 330 329 0 +279 280 331 0 +279 331 330 0 +280 281 332 0 +280 332 331 0 +281 282 333 0 +281 333 332 0 +282 283 334 0 +282 334 333 0 +283 284 335 0 +283 335 334 0 +284 285 336 0 +284 336 335 0 +285 286 337 0 +285 337 336 0 +286 287 338 0 +286 338 337 0 +287 288 339 0 +287 339 338 0 +288 289 340 0 +288 340 339 0 +289 290 341 0 +289 341 340 0 +290 291 342 0 +290 342 341 0 +291 292 343 0 +291 343 342 0 +292 293 344 0 +292 344 343 0 +293 294 345 0 +293 345 344 0 +294 295 346 0 +294 346 345 0 +295 296 347 0 +295 347 346 0 +296 297 348 0 +296 348 347 0 +297 298 349 0 +297 349 348 0 +298 299 350 0 +298 350 349 0 +299 300 351 0 +299 351 350 0 +300 301 352 0 +300 352 351 0 +301 302 353 0 +301 353 352 0 +302 303 354 0 +302 354 353 0 +303 304 355 0 +303 355 354 0 +304 305 356 0 +304 356 355 0 +305 306 357 0 +305 357 356 0 +307 308 359 0 +307 359 358 0 +308 309 360 0 +308 360 359 0 +309 310 361 0 +309 361 360 0 +310 311 362 0 +310 362 361 0 +311 312 363 0 +311 363 362 0 +312 313 364 0 +312 364 363 0 +313 314 365 0 +313 365 364 0 +314 315 366 0 +314 366 365 0 +315 316 367 0 +315 367 366 0 +316 317 368 0 +316 368 367 0 +317 318 369 0 +317 369 368 0 +318 319 370 0 +318 370 369 0 +319 320 371 0 +319 371 370 0 +320 321 372 0 +320 372 371 0 +321 322 373 0 +321 373 372 0 +322 323 374 0 +322 374 373 0 +323 324 375 0 +323 375 374 0 +324 325 376 0 +324 376 375 0 +325 326 377 0 +325 377 376 0 +326 327 378 0 +326 378 377 0 +327 328 379 0 +327 379 378 0 +328 329 380 0 +328 380 379 0 +329 330 381 0 +329 381 380 0 +330 331 382 0 +330 382 381 0 +331 332 383 0 +331 383 382 0 +332 333 384 0 +332 384 383 0 +333 334 385 0 +333 385 384 0 +334 335 386 0 +334 386 385 0 +335 336 387 0 +335 387 386 0 +336 337 388 0 +336 388 387 0 +337 338 389 0 +337 389 388 0 +338 339 390 0 +338 390 389 0 +339 340 391 0 +339 391 390 0 +340 341 392 0 +340 392 391 0 +341 342 393 0 +341 393 392 0 +342 343 394 0 +342 394 393 0 +343 344 395 0 +343 395 394 0 +344 345 396 0 +344 396 395 0 +345 346 397 0 +345 397 396 0 +346 347 398 0 +346 398 397 0 +347 348 399 0 +347 399 398 0 +348 349 400 0 +348 400 399 0 +349 350 401 0 +349 401 400 0 +350 351 402 0 +350 402 401 0 +351 352 403 0 +351 403 402 0 +352 353 404 0 +352 404 403 0 +353 354 405 0 +353 405 404 0 +354 355 406 0 +354 406 405 0 +355 356 407 0 +355 407 406 0 +356 357 408 0 +356 408 407 0 +358 359 410 0 +358 410 409 0 +359 360 411 0 +359 411 410 0 +360 361 412 0 +360 412 411 0 +361 362 413 0 +361 413 412 0 +362 363 414 0 +362 414 413 0 +363 364 415 0 +363 415 414 0 +364 365 416 0 +364 416 415 0 +365 366 417 0 +365 417 416 0 +366 367 418 0 +366 418 417 0 +367 368 419 0 +367 419 418 0 +368 369 420 0 +368 420 419 0 +369 370 421 0 +369 421 420 0 +370 371 422 0 +370 422 421 0 +371 372 423 0 +371 423 422 0 +372 373 424 0 +372 424 423 0 +373 374 425 0 +373 425 424 0 +374 375 426 0 +374 426 425 0 +375 376 427 0 +375 427 426 0 +376 377 428 0 +376 428 427 0 +377 378 429 0 +377 429 428 0 +378 379 430 0 +378 430 429 0 +379 380 431 0 +379 431 430 0 +380 381 432 0 +380 432 431 0 +381 382 433 0 +381 433 432 0 +382 383 434 0 +382 434 433 0 +383 384 435 0 +383 435 434 0 +384 385 436 0 +384 436 435 0 +385 386 437 0 +385 437 436 0 +386 387 438 0 +386 438 437 0 +387 388 439 0 +387 439 438 0 +388 389 440 0 +388 440 439 0 +389 390 441 0 +389 441 440 0 +390 391 442 0 +390 442 441 0 +391 392 443 0 +391 443 442 0 +392 393 444 0 +392 444 443 0 +393 394 445 0 +393 445 444 0 +394 395 446 0 +394 446 445 0 +395 396 447 0 +395 447 446 0 +396 397 448 0 +396 448 447 0 +397 398 449 0 +397 449 448 0 +398 399 450 0 +398 450 449 0 +399 400 451 0 +399 451 450 0 +400 401 452 0 +400 452 451 0 +401 402 453 0 +401 453 452 0 +402 403 454 0 +402 454 453 0 +403 404 455 0 +403 455 454 0 +404 405 456 0 +404 456 455 0 +405 406 457 0 +405 457 456 0 +406 407 458 0 +406 458 457 0 +407 408 459 0 +407 459 458 0 +409 410 461 0 +409 461 460 0 +410 411 462 0 +410 462 461 0 +411 412 463 0 +411 463 462 0 +412 413 464 0 +412 464 463 0 +413 414 465 0 +413 465 464 0 +414 415 466 0 +414 466 465 0 +415 416 467 0 +415 467 466 0 +416 417 468 0 +416 468 467 0 +417 418 469 0 +417 469 468 0 +418 419 470 0 +418 470 469 0 +419 420 471 0 +419 471 470 0 +420 421 472 0 +420 472 471 0 +421 422 473 0 +421 473 472 0 +422 423 474 0 +422 474 473 0 +423 424 475 0 +423 475 474 0 +424 425 476 0 +424 476 475 0 +425 426 477 0 +425 477 476 0 +426 427 478 0 +426 478 477 0 +427 428 479 0 +427 479 478 0 +428 429 480 0 +428 480 479 0 +429 430 481 0 +429 481 480 0 +430 431 482 0 +430 482 481 0 +431 432 483 0 +431 483 482 0 +432 433 484 0 +432 484 483 0 +433 434 485 0 +433 485 484 0 +434 435 486 0 +434 486 485 0 +435 436 487 0 +435 487 486 0 +436 437 488 0 +436 488 487 0 +437 438 489 0 +437 489 488 0 +438 439 490 0 +438 490 489 0 +439 440 491 0 +439 491 490 0 +440 441 492 0 +440 492 491 0 +441 442 493 0 +441 493 492 0 +442 443 494 0 +442 494 493 0 +443 444 495 0 +443 495 494 0 +444 445 496 0 +444 496 495 0 +445 446 497 0 +445 497 496 0 +446 447 498 0 +446 498 497 0 +447 448 499 0 +447 499 498 0 +448 449 500 0 +448 500 499 0 +449 450 501 0 +449 501 500 0 +450 451 502 0 +450 502 501 0 +451 452 503 0 +451 503 502 0 +452 453 504 0 +452 504 503 0 +453 454 505 0 +453 505 504 0 +454 455 506 0 +454 506 505 0 +455 456 507 0 +455 507 506 0 +456 457 508 0 +456 508 507 0 +457 458 509 0 +457 509 508 0 +458 459 510 0 +458 510 509 0 +460 461 512 0 +460 512 511 0 +461 462 513 0 +461 513 512 0 +462 463 514 0 +462 514 513 0 +463 464 515 0 +463 515 514 0 +464 465 516 0 +464 516 515 0 +465 466 517 0 +465 517 516 0 +466 467 518 0 +466 518 517 0 +467 468 519 0 +467 519 518 0 +468 469 520 0 +468 520 519 0 +469 470 521 0 +469 521 520 0 +470 471 522 0 +470 522 521 0 +471 472 523 0 +471 523 522 0 +472 473 524 0 +472 524 523 0 +473 474 525 0 +473 525 524 0 +474 475 526 0 +474 526 525 0 +475 476 527 0 +475 527 526 0 +476 477 528 0 +476 528 527 0 +477 478 529 0 +477 529 528 0 +478 479 530 0 +478 530 529 0 +479 480 531 0 +479 531 530 0 +480 481 532 0 +480 532 531 0 +481 482 533 0 +481 533 532 0 +482 483 534 0 +482 534 533 0 +483 484 535 0 +483 535 534 0 +484 485 536 0 +484 536 535 0 +485 486 537 0 +485 537 536 0 +486 487 538 0 +486 538 537 0 +487 488 539 0 +487 539 538 0 +488 489 540 0 +488 540 539 0 +489 490 541 0 +489 541 540 0 +490 491 542 0 +490 542 541 0 +491 492 543 0 +491 543 542 0 +492 493 544 0 +492 544 543 0 +493 494 545 0 +493 545 544 0 +494 495 546 0 +494 546 545 0 +495 496 547 0 +495 547 546 0 +496 497 548 0 +496 548 547 0 +497 498 549 0 +497 549 548 0 +498 499 550 0 +498 550 549 0 +499 500 551 0 +499 551 550 0 +500 501 552 0 +500 552 551 0 +501 502 553 0 +501 553 552 0 +502 503 554 0 +502 554 553 0 +503 504 555 0 +503 555 554 0 +504 505 556 0 +504 556 555 0 +505 506 557 0 +505 557 556 0 +506 507 558 0 +506 558 557 0 +507 508 559 0 +507 559 558 0 +508 509 560 0 +508 560 559 0 +509 510 561 0 +509 561 560 0 +511 512 563 0 +511 563 562 0 +512 513 564 0 +512 564 563 0 +513 514 565 0 +513 565 564 0 +514 515 566 0 +514 566 565 0 +515 516 567 0 +515 567 566 0 +516 517 568 0 +516 568 567 0 +517 518 569 0 +517 569 568 0 +518 519 570 0 +518 570 569 0 +519 520 571 0 +519 571 570 0 +520 521 572 0 +520 572 571 0 +521 522 573 0 +521 573 572 0 +522 523 574 0 +522 574 573 0 +523 524 575 0 +523 575 574 0 +524 525 576 0 +524 576 575 0 +525 526 577 0 +525 577 576 0 +526 527 578 0 +526 578 577 0 +527 528 579 0 +527 579 578 0 +528 529 580 0 +528 580 579 0 +529 530 581 0 +529 581 580 0 +530 531 582 0 +530 582 581 0 +531 532 583 0 +531 583 582 0 +532 533 584 0 +532 584 583 0 +533 534 585 0 +533 585 584 0 +534 535 586 0 +534 586 585 0 +535 536 587 0 +535 587 586 0 +536 537 588 0 +536 588 587 0 +537 538 589 0 +537 589 588 0 +538 539 590 0 +538 590 589 0 +539 540 591 0 +539 591 590 0 +540 541 592 0 +540 592 591 0 +541 542 593 0 +541 593 592 0 +542 543 594 0 +542 594 593 0 +543 544 595 0 +543 595 594 0 +544 545 596 0 +544 596 595 0 +545 546 597 0 +545 597 596 0 +546 547 598 0 +546 598 597 0 +547 548 599 0 +547 599 598 0 +548 549 600 0 +548 600 599 0 +549 550 601 0 +549 601 600 0 +550 551 602 0 +550 602 601 0 +551 552 603 0 +551 603 602 0 +552 553 604 0 +552 604 603 0 +553 554 605 0 +553 605 604 0 +554 555 606 0 +554 606 605 0 +555 556 607 0 +555 607 606 0 +556 557 608 0 +556 608 607 0 +557 558 609 0 +557 609 608 0 +558 559 610 0 +558 610 609 0 +559 560 611 0 +559 611 610 0 +560 561 612 0 +560 612 611 0 +562 563 614 0 +562 614 613 0 +563 564 615 0 +563 615 614 0 +564 565 616 0 +564 616 615 0 +565 566 617 0 +565 617 616 0 +566 567 618 0 +566 618 617 0 +567 568 619 0 +567 619 618 0 +568 569 620 0 +568 620 619 0 +569 570 621 0 +569 621 620 0 +570 571 622 0 +570 622 621 0 +571 572 623 0 +571 623 622 0 +572 573 624 0 +572 624 623 0 +573 574 625 0 +573 625 624 0 +574 575 626 0 +574 626 625 0 +575 576 627 0 +575 627 626 0 +576 577 628 0 +576 628 627 0 +577 578 629 0 +577 629 628 0 +578 579 630 0 +578 630 629 0 +579 580 631 0 +579 631 630 0 +580 581 632 0 +580 632 631 0 +581 582 633 0 +581 633 632 0 +582 583 634 0 +582 634 633 0 +583 584 635 0 +583 635 634 0 +584 585 636 0 +584 636 635 0 +585 586 637 0 +585 637 636 0 +586 587 638 0 +586 638 637 0 +587 588 639 0 +587 639 638 0 +588 589 640 0 +588 640 639 0 +589 590 641 0 +589 641 640 0 +590 591 642 0 +590 642 641 0 +591 592 643 0 +591 643 642 0 +592 593 644 0 +592 644 643 0 +593 594 645 0 +593 645 644 0 +594 595 646 0 +594 646 645 0 +595 596 647 0 +595 647 646 0 +596 597 648 0 +596 648 647 0 +597 598 649 0 +597 649 648 0 +598 599 650 0 +598 650 649 0 +599 600 651 0 +599 651 650 0 +600 601 652 0 +600 652 651 0 +601 602 653 0 +601 653 652 0 +602 603 654 0 +602 654 653 0 +603 604 655 0 +603 655 654 0 +604 605 656 0 +604 656 655 0 +605 606 657 0 +605 657 656 0 +606 607 658 0 +606 658 657 0 +607 608 659 0 +607 659 658 0 +608 609 660 0 +608 660 659 0 +609 610 661 0 +609 661 660 0 +610 611 662 0 +610 662 661 0 +611 612 663 0 +611 663 662 0 +613 614 665 0 +613 665 664 0 +614 615 666 0 +614 666 665 0 +615 616 667 0 +615 667 666 0 +616 617 668 0 +616 668 667 0 +617 618 669 0 +617 669 668 0 +618 619 670 0 +618 670 669 0 +619 620 671 0 +619 671 670 0 +620 621 672 0 +620 672 671 0 +621 622 673 0 +621 673 672 0 +622 623 674 0 +622 674 673 0 +623 624 675 0 +623 675 674 0 +624 625 676 0 +624 676 675 0 +625 626 677 0 +625 677 676 0 +626 627 678 0 +626 678 677 0 +627 628 679 0 +627 679 678 0 +628 629 680 0 +628 680 679 0 +629 630 681 0 +629 681 680 0 +630 631 682 0 +630 682 681 0 +631 632 683 0 +631 683 682 0 +632 633 684 0 +632 684 683 0 +633 634 685 0 +633 685 684 0 +634 635 686 0 +634 686 685 0 +635 636 687 0 +635 687 686 0 +636 637 688 0 +636 688 687 0 +637 638 689 0 +637 689 688 0 +638 639 690 0 +638 690 689 0 +639 640 691 0 +639 691 690 0 +640 641 692 0 +640 692 691 0 +641 642 693 0 +641 693 692 0 +642 643 694 0 +642 694 693 0 +643 644 695 0 +643 695 694 0 +644 645 696 0 +644 696 695 0 +645 646 697 0 +645 697 696 0 +646 647 698 0 +646 698 697 0 +647 648 699 0 +647 699 698 0 +648 649 700 0 +648 700 699 0 +649 650 701 0 +649 701 700 0 +650 651 702 0 +650 702 701 0 +651 652 703 0 +651 703 702 0 +652 653 704 0 +652 704 703 0 +653 654 705 0 +653 705 704 0 +654 655 706 0 +654 706 705 0 +655 656 707 0 +655 707 706 0 +656 657 708 0 +656 708 707 0 +657 658 709 0 +657 709 708 0 +658 659 710 0 +658 710 709 0 +659 660 711 0 +659 711 710 0 +660 661 712 0 +660 712 711 0 +661 662 713 0 +661 713 712 0 +662 663 714 0 +662 714 713 0 +664 665 716 0 +664 716 715 0 +665 666 717 0 +665 717 716 0 +666 667 718 0 +666 718 717 0 +667 668 719 0 +667 719 718 0 +668 669 720 0 +668 720 719 0 +669 670 721 0 +669 721 720 0 +670 671 722 0 +670 722 721 0 +671 672 723 0 +671 723 722 0 +672 673 724 0 +672 724 723 0 +673 674 725 0 +673 725 724 0 +674 675 726 0 +674 726 725 0 +675 676 727 0 +675 727 726 0 +676 677 728 0 +676 728 727 0 +677 678 729 0 +677 729 728 0 +678 679 730 0 +678 730 729 0 +679 680 731 0 +679 731 730 0 +680 681 732 0 +680 732 731 0 +681 682 733 0 +681 733 732 0 +682 683 734 0 +682 734 733 0 +683 684 735 0 +683 735 734 0 +684 685 736 0 +684 736 735 0 +685 686 737 0 +685 737 736 0 +686 687 738 0 +686 738 737 0 +687 688 739 0 +687 739 738 0 +688 689 740 0 +688 740 739 0 +689 690 741 0 +689 741 740 0 +690 691 742 0 +690 742 741 0 +691 692 743 0 +691 743 742 0 +692 693 744 0 +692 744 743 0 +693 694 745 0 +693 745 744 0 +694 695 746 0 +694 746 745 0 +695 696 747 0 +695 747 746 0 +696 697 748 0 +696 748 747 0 +697 698 749 0 +697 749 748 0 +698 699 750 0 +698 750 749 0 +699 700 751 0 +699 751 750 0 +700 701 752 0 +700 752 751 0 +701 702 753 0 +701 753 752 0 +702 703 754 0 +702 754 753 0 +703 704 755 0 +703 755 754 0 +704 705 756 0 +704 756 755 0 +705 706 757 0 +705 757 756 0 +706 707 758 0 +706 758 757 0 +707 708 759 0 +707 759 758 0 +708 709 760 0 +708 760 759 0 +709 710 761 0 +709 761 760 0 +710 711 762 0 +710 762 761 0 +711 712 763 0 +711 763 762 0 +712 713 764 0 +712 764 763 0 +713 714 765 0 +713 765 764 0 +715 716 767 0 +715 767 766 0 +716 717 768 0 +716 768 767 0 +717 718 769 0 +717 769 768 0 +718 719 770 0 +718 770 769 0 +719 720 771 0 +719 771 770 0 +720 721 772 0 +720 772 771 0 +721 722 773 0 +721 773 772 0 +722 723 774 0 +722 774 773 0 +723 724 775 0 +723 775 774 0 +724 725 776 0 +724 776 775 0 +725 726 777 0 +725 777 776 0 +726 727 778 0 +726 778 777 0 +727 728 779 0 +727 779 778 0 +728 729 780 0 +728 780 779 0 +729 730 781 0 +729 781 780 0 +730 731 782 0 +730 782 781 0 +731 732 783 0 +731 783 782 0 +732 733 784 0 +732 784 783 0 +733 734 785 0 +733 785 784 0 +734 735 786 0 +734 786 785 0 +735 736 787 0 +735 787 786 0 +736 737 788 0 +736 788 787 0 +737 738 789 0 +737 789 788 0 +738 739 790 0 +738 790 789 0 +739 740 791 0 +739 791 790 0 +740 741 792 0 +740 792 791 0 +741 742 793 0 +741 793 792 0 +742 743 794 0 +742 794 793 0 +743 744 795 0 +743 795 794 0 +744 745 796 0 +744 796 795 0 +745 746 797 0 +745 797 796 0 +746 747 798 0 +746 798 797 0 +747 748 799 0 +747 799 798 0 +748 749 800 0 +748 800 799 0 +749 750 801 0 +749 801 800 0 +750 751 802 0 +750 802 801 0 +751 752 803 0 +751 803 802 0 +752 753 804 0 +752 804 803 0 +753 754 805 0 +753 805 804 0 +754 755 806 0 +754 806 805 0 +755 756 807 0 +755 807 806 0 +756 757 808 0 +756 808 807 0 +757 758 809 0 +757 809 808 0 +758 759 810 0 +758 810 809 0 +759 760 811 0 +759 811 810 0 +760 761 812 0 +760 812 811 0 +761 762 813 0 +761 813 812 0 +762 763 814 0 +762 814 813 0 +763 764 815 0 +763 815 814 0 +764 765 816 0 +764 816 815 0 +766 767 818 0 +766 818 817 0 +767 768 819 0 +767 819 818 0 +768 769 820 0 +768 820 819 0 +769 770 821 0 +769 821 820 0 +770 771 822 0 +770 822 821 0 +771 772 823 0 +771 823 822 0 +772 773 824 0 +772 824 823 0 +773 774 825 0 +773 825 824 0 +774 775 826 0 +774 826 825 0 +775 776 827 0 +775 827 826 0 +776 777 828 0 +776 828 827 0 +777 778 829 0 +777 829 828 0 +778 779 830 0 +778 830 829 0 +779 780 831 0 +779 831 830 0 +780 781 832 0 +780 832 831 0 +781 782 833 0 +781 833 832 0 +782 783 834 0 +782 834 833 0 +783 784 835 0 +783 835 834 0 +784 785 836 0 +784 836 835 0 +785 786 837 0 +785 837 836 0 +786 787 838 0 +786 838 837 0 +787 788 839 0 +787 839 838 0 +788 789 840 0 +788 840 839 0 +789 790 841 0 +789 841 840 0 +790 791 842 0 +790 842 841 0 +791 792 843 0 +791 843 842 0 +792 793 844 0 +792 844 843 0 +793 794 845 0 +793 845 844 0 +794 795 846 0 +794 846 845 0 +795 796 847 0 +795 847 846 0 +796 797 848 0 +796 848 847 0 +797 798 849 0 +797 849 848 0 +798 799 850 0 +798 850 849 0 +799 800 851 0 +799 851 850 0 +800 801 852 0 +800 852 851 0 +801 802 853 0 +801 853 852 0 +802 803 854 0 +802 854 853 0 +803 804 855 0 +803 855 854 0 +804 805 856 0 +804 856 855 0 +805 806 857 0 +805 857 856 0 +806 807 858 0 +806 858 857 0 +807 808 859 0 +807 859 858 0 +808 809 860 0 +808 860 859 0 +809 810 861 0 +809 861 860 0 +810 811 862 0 +810 862 861 0 +811 812 863 0 +811 863 862 0 +812 813 864 0 +812 864 863 0 +813 814 865 0 +813 865 864 0 +814 815 866 0 +814 866 865 0 +815 816 867 0 +815 867 866 0 +817 818 869 0 +817 869 868 0 +818 819 870 0 +818 870 869 0 +819 820 871 0 +819 871 870 0 +820 821 872 0 +820 872 871 0 +821 822 873 0 +821 873 872 0 +822 823 874 0 +822 874 873 0 +823 824 875 0 +823 875 874 0 +824 825 876 0 +824 876 875 0 +825 826 877 0 +825 877 876 0 +826 827 878 0 +826 878 877 0 +827 828 879 0 +827 879 878 0 +828 829 880 0 +828 880 879 0 +829 830 881 0 +829 881 880 0 +830 831 882 0 +830 882 881 0 +831 832 883 0 +831 883 882 0 +832 833 884 0 +832 884 883 0 +833 834 885 0 +833 885 884 0 +834 835 886 0 +834 886 885 0 +835 836 887 0 +835 887 886 0 +836 837 888 0 +836 888 887 0 +837 838 889 0 +837 889 888 0 +838 839 890 0 +838 890 889 0 +839 840 891 0 +839 891 890 0 +840 841 892 0 +840 892 891 0 +841 842 893 0 +841 893 892 0 +842 843 894 0 +842 894 893 0 +843 844 895 0 +843 895 894 0 +844 845 896 0 +844 896 895 0 +845 846 897 0 +845 897 896 0 +846 847 898 0 +846 898 897 0 +847 848 899 0 +847 899 898 0 +848 849 900 0 +848 900 899 0 +849 850 901 0 +849 901 900 0 +850 851 902 0 +850 902 901 0 +851 852 903 0 +851 903 902 0 +852 853 904 0 +852 904 903 0 +853 854 905 0 +853 905 904 0 +854 855 906 0 +854 906 905 0 +855 856 907 0 +855 907 906 0 +856 857 908 0 +856 908 907 0 +857 858 909 0 +857 909 908 0 +858 859 910 0 +858 910 909 0 +859 860 911 0 +859 911 910 0 +860 861 912 0 +860 912 911 0 +861 862 913 0 +861 913 912 0 +862 863 914 0 +862 914 913 0 +863 864 915 0 +863 915 914 0 +864 865 916 0 +864 916 915 0 +865 866 917 0 +865 917 916 0 +866 867 918 0 +866 918 917 0 +868 869 920 0 +868 920 919 0 +869 870 921 0 +869 921 920 0 +870 871 922 0 +870 922 921 0 +871 872 923 0 +871 923 922 0 +872 873 924 0 +872 924 923 0 +873 874 925 0 +873 925 924 0 +874 875 926 0 +874 926 925 0 +875 876 927 0 +875 927 926 0 +876 877 928 0 +876 928 927 0 +877 878 929 0 +877 929 928 0 +878 879 930 0 +878 930 929 0 +879 880 931 0 +879 931 930 0 +880 881 932 0 +880 932 931 0 +881 882 933 0 +881 933 932 0 +882 883 934 0 +882 934 933 0 +883 884 935 0 +883 935 934 0 +884 885 936 0 +884 936 935 0 +885 886 937 0 +885 937 936 0 +886 887 938 0 +886 938 937 0 +887 888 939 0 +887 939 938 0 +888 889 940 0 +888 940 939 0 +889 890 941 0 +889 941 940 0 +890 891 942 0 +890 942 941 0 +891 892 943 0 +891 943 942 0 +892 893 944 0 +892 944 943 0 +893 894 945 0 +893 945 944 0 +894 895 946 0 +894 946 945 0 +895 896 947 0 +895 947 946 0 +896 897 948 0 +896 948 947 0 +897 898 949 0 +897 949 948 0 +898 899 950 0 +898 950 949 0 +899 900 951 0 +899 951 950 0 +900 901 952 0 +900 952 951 0 +901 902 953 0 +901 953 952 0 +902 903 954 0 +902 954 953 0 +903 904 955 0 +903 955 954 0 +904 905 956 0 +904 956 955 0 +905 906 957 0 +905 957 956 0 +906 907 958 0 +906 958 957 0 +907 908 959 0 +907 959 958 0 +908 909 960 0 +908 960 959 0 +909 910 961 0 +909 961 960 0 +910 911 962 0 +910 962 961 0 +911 912 963 0 +911 963 962 0 +912 913 964 0 +912 964 963 0 +913 914 965 0 +913 965 964 0 +914 915 966 0 +914 966 965 0 +915 916 967 0 +915 967 966 0 +916 917 968 0 +916 968 967 0 +917 918 969 0 +917 969 968 0 +919 920 971 0 +919 971 970 0 +920 921 972 0 +920 972 971 0 +921 922 973 0 +921 973 972 0 +922 923 974 0 +922 974 973 0 +923 924 975 0 +923 975 974 0 +924 925 976 0 +924 976 975 0 +925 926 977 0 +925 977 976 0 +926 927 978 0 +926 978 977 0 +927 928 979 0 +927 979 978 0 +928 929 980 0 +928 980 979 0 +929 930 981 0 +929 981 980 0 +930 931 982 0 +930 982 981 0 +931 932 983 0 +931 983 982 0 +932 933 984 0 +932 984 983 0 +933 934 985 0 +933 985 984 0 +934 935 986 0 +934 986 985 0 +935 936 987 0 +935 987 986 0 +936 937 988 0 +936 988 987 0 +937 938 989 0 +937 989 988 0 +938 939 990 0 +938 990 989 0 +939 940 991 0 +939 991 990 0 +940 941 992 0 +940 992 991 0 +941 942 993 0 +941 993 992 0 +942 943 994 0 +942 994 993 0 +943 944 995 0 +943 995 994 0 +944 945 996 0 +944 996 995 0 +945 946 997 0 +945 997 996 0 +946 947 998 0 +946 998 997 0 +947 948 999 0 +947 999 998 0 +948 949 1000 0 +948 1000 999 0 +949 950 1001 0 +949 1001 1000 0 +950 951 1002 0 +950 1002 1001 0 +951 952 1003 0 +951 1003 1002 0 +952 953 1004 0 +952 1004 1003 0 +953 954 1005 0 +953 1005 1004 0 +954 955 1006 0 +954 1006 1005 0 +955 956 1007 0 +955 1007 1006 0 +956 957 1008 0 +956 1008 1007 0 +957 958 1009 0 +957 1009 1008 0 +958 959 1010 0 +958 1010 1009 0 +959 960 1011 0 +959 1011 1010 0 +960 961 1012 0 +960 1012 1011 0 +961 962 1013 0 +961 1013 1012 0 +962 963 1014 0 +962 1014 1013 0 +963 964 1015 0 +963 1015 1014 0 +964 965 1016 0 +964 1016 1015 0 +965 966 1017 0 +965 1017 1016 0 +966 967 1018 0 +966 1018 1017 0 +967 968 1019 0 +967 1019 1018 0 +968 969 1020 0 +968 1020 1019 0 +970 971 1022 0 +970 1022 1021 0 +971 972 1023 0 +971 1023 1022 0 +972 973 1024 0 +972 1024 1023 0 +973 974 1025 0 +973 1025 1024 0 +974 975 1026 0 +974 1026 1025 0 +975 976 1027 0 +975 1027 1026 0 +976 977 1028 0 +976 1028 1027 0 +977 978 1029 0 +977 1029 1028 0 +978 979 1030 0 +978 1030 1029 0 +979 980 1031 0 +979 1031 1030 0 +980 981 1032 0 +980 1032 1031 0 +981 982 1033 0 +981 1033 1032 0 +982 983 1034 0 +982 1034 1033 0 +983 984 1035 0 +983 1035 1034 0 +984 985 1036 0 +984 1036 1035 0 +985 986 1037 0 +985 1037 1036 0 +986 987 1038 0 +986 1038 1037 0 +987 988 1039 0 +987 1039 1038 0 +988 989 1040 0 +988 1040 1039 0 +989 990 1041 0 +989 1041 1040 0 +990 991 1042 0 +990 1042 1041 0 +991 992 1043 0 +991 1043 1042 0 +992 993 1044 0 +992 1044 1043 0 +993 994 1045 0 +993 1045 1044 0 +994 995 1046 0 +994 1046 1045 0 +995 996 1047 0 +995 1047 1046 0 +996 997 1048 0 +996 1048 1047 0 +997 998 1049 0 +997 1049 1048 0 +998 999 1050 0 +998 1050 1049 0 +999 1000 1051 0 +999 1051 1050 0 +1000 1001 1052 0 +1000 1052 1051 0 +1001 1002 1053 0 +1001 1053 1052 0 +1002 1003 1054 0 +1002 1054 1053 0 +1003 1004 1055 0 +1003 1055 1054 0 +1004 1005 1056 0 +1004 1056 1055 0 +1005 1006 1057 0 +1005 1057 1056 0 +1006 1007 1058 0 +1006 1058 1057 0 +1007 1008 1059 0 +1007 1059 1058 0 +1008 1009 1060 0 +1008 1060 1059 0 +1009 1010 1061 0 +1009 1061 1060 0 +1010 1011 1062 0 +1010 1062 1061 0 +1011 1012 1063 0 +1011 1063 1062 0 +1012 1013 1064 0 +1012 1064 1063 0 +1013 1014 1065 0 +1013 1065 1064 0 +1014 1015 1066 0 +1014 1066 1065 0 +1015 1016 1067 0 +1015 1067 1066 0 +1016 1017 1068 0 +1016 1068 1067 0 +1017 1018 1069 0 +1017 1069 1068 0 +1018 1019 1070 0 +1018 1070 1069 0 +1019 1020 1071 0 +1019 1071 1070 0 + +End \ No newline at end of file diff --git a/Data/Mesh/nonlinearmembrane_rectangle.mesh b/Data/Mesh/nonlinearmembrane_rectangle.mesh new file mode 100644 index 0000000000000000000000000000000000000000..0e6876dad8407ad12d1d976a59384829c5595231 --- /dev/null +++ b/Data/Mesh/nonlinearmembrane_rectangle.mesh @@ -0,0 +1,566 @@ +MeshVersionFormatted 2 +Dimension +3 +Vertices +186 +0 0 0 1 +0.2 0 0 1 +0.4 0 0 1 +0.6 0 0 1 +0.8 0 0 1 +1 0 0 1 +1.2 0 0 1 +1.4 0 0 1 +1.6 0 0 1 +1.8 0 0 1 +2 0 0 1 +2.2 0 0 1 +2.4 0 0 1 +2.6 0 0 1 +2.8 0 0 1 +3 0 0 1 +3.2 0 0 1 +3.4 0 0 1 +3.6 0 0 1 +3.8 0 0 1 +4 0 0 1 +4.2 0 0 1 +4.4 0 0 1 +4.6 0 0 1 +4.8 0 0 1 +5 0 0 1 +5.2 0 0 1 +5.4 0 0 1 +5.6 0 0 1 +5.8 0 0 1 +6 0 0 1 +0 0.2 0 1 +0.2 0.2 0 1 +0.4 0.2 0 1 +0.6 0.2 0 1 +0.8 0.2 0 1 +1 0.2 0 1 +1.2 0.2 0 1 +1.4 0.2 0 1 +1.6 0.2 0 1 +1.8 0.2 0 1 +2 0.2 0 1 +2.2 0.2 0 1 +2.4 0.2 0 1 +2.6 0.2 0 1 +2.8 0.2 0 1 +3 0.2 0 1 +3.2 0.2 0 1 +3.4 0.2 0 1 +3.6 0.2 0 1 +3.8 0.2 0 1 +4 0.2 0 1 +4.2 0.2 0 1 +4.4 0.2 0 1 +4.6 0.2 0 1 +4.8 0.2 0 1 +5 0.2 0 1 +5.2 0.2 0 1 +5.4 0.2 0 1 +5.6 0.2 0 1 +5.8 0.2 0 1 +6 0.2 0 1 +0 0.4 0 1 +0.2 0.4 0 1 +0.4 0.4 0 1 +0.6 0.4 0 1 +0.8 0.4 0 1 +1 0.4 0 1 +1.2 0.4 0 1 +1.4 0.4 0 1 +1.6 0.4 0 1 +1.8 0.4 0 1 +2 0.4 0 1 +2.2 0.4 0 1 +2.4 0.4 0 1 +2.6 0.4 0 1 +2.8 0.4 0 1 +3 0.4 0 1 +3.2 0.4 0 1 +3.4 0.4 0 1 +3.6 0.4 0 1 +3.8 0.4 0 1 +4 0.4 0 1 +4.2 0.4 0 1 +4.4 0.4 0 1 +4.6 0.4 0 1 +4.8 0.4 0 1 +5 0.4 0 1 +5.2 0.4 0 1 +5.4 0.4 0 1 +5.6 0.4 0 1 +5.8 0.4 0 1 +6 0.4 0 1 +0 0.6 0 1 +0.2 0.6 0 1 +0.4 0.6 0 1 +0.6 0.6 0 1 +0.8 0.6 0 1 +1 0.6 0 1 +1.2 0.6 0 1 +1.4 0.6 0 1 +1.6 0.6 0 1 +1.8 0.6 0 1 +2 0.6 0 1 +2.2 0.6 0 1 +2.4 0.6 0 1 +2.6 0.6 0 1 +2.8 0.6 0 1 +3 0.6 0 1 +3.2 0.6 0 1 +3.4 0.6 0 1 +3.6 0.6 0 1 +3.8 0.6 0 1 +4 0.6 0 1 +4.2 0.6 0 1 +4.4 0.6 0 1 +4.6 0.6 0 1 +4.8 0.6 0 1 +5 0.6 0 1 +5.2 0.6 0 1 +5.4 0.6 0 1 +5.6 0.6 0 1 +5.8 0.6 0 1 +6 0.6 0 1 +0 0.8 0 1 +0.2 0.8 0 1 +0.4 0.8 0 1 +0.6 0.8 0 1 +0.8 0.8 0 1 +1 0.8 0 1 +1.2 0.8 0 1 +1.4 0.8 0 1 +1.6 0.8 0 1 +1.8 0.8 0 1 +2 0.8 0 1 +2.2 0.8 0 1 +2.4 0.8 0 1 +2.6 0.8 0 1 +2.8 0.8 0 1 +3 0.8 0 1 +3.2 0.8 0 1 +3.4 0.8 0 1 +3.6 0.8 0 1 +3.8 0.8 0 1 +4 0.8 0 1 +4.2 0.8 0 1 +4.4 0.8 0 1 +4.6 0.8 0 1 +4.8 0.8 0 1 +5 0.8 0 1 +5.2 0.8 0 1 +5.4 0.8 0 1 +5.6 0.8 0 1 +5.8 0.8 0 1 +6 0.8 0 1 +0 1 0 1 +0.2 1 0 1 +0.4 1 0 1 +0.6 1 0 1 +0.8 1 0 1 +1 1 0 1 +1.2 1 0 1 +1.4 1 0 1 +1.6 1 0 1 +1.8 1 0 1 +2 1 0 1 +2.2 1 0 1 +2.4 1 0 1 +2.6 1 0 1 +2.8 1 0 1 +3 1 0 1 +3.2 1 0 1 +3.4 1 0 1 +3.6 1 0 1 +3.8 1 0 1 +4 1 0 1 +4.2 1 0 1 +4.4 1 0 1 +4.6 1 0 1 +4.8 1 0 1 +5 1 0 1 +5.2 1 0 1 +5.4 1 0 1 +5.6 1 0 1 +5.8 1 0 1 +6 1 0 1 +Edges +70 +1 2 1 +2 3 1 +3 4 1 +4 5 1 +5 6 1 +6 7 1 +7 8 1 +8 9 1 +9 10 1 +10 11 1 +11 12 1 +12 13 1 +13 14 1 +14 15 1 +15 16 1 +16 17 1 +17 18 1 +18 19 1 +19 20 1 +20 21 1 +21 22 1 +22 23 1 +23 24 1 +24 25 1 +25 26 1 +26 27 1 +27 28 1 +28 29 1 +29 30 1 +30 31 1 +31 62 2 +62 93 2 +93 124 2 +124 155 2 +155 186 2 +157 156 3 +158 157 3 +159 158 3 +160 159 3 +161 160 3 +162 161 3 +163 162 3 +164 163 3 +165 164 3 +166 165 3 +167 166 3 +168 167 3 +169 168 3 +170 169 3 +171 170 3 +172 171 3 +173 172 3 +174 173 3 +175 174 3 +176 175 3 +177 176 3 +178 177 3 +179 178 3 +180 179 3 +181 180 3 +182 181 3 +183 182 3 +184 183 3 +185 184 3 +186 185 3 +32 1 4 +63 32 4 +94 63 4 +125 94 4 +156 125 4 +Triangles +300 +1 2 33 0 +1 33 32 0 +2 3 34 0 +2 34 33 0 +3 4 35 0 +3 35 34 0 +4 5 36 0 +4 36 35 0 +5 6 37 0 +5 37 36 0 +6 7 38 0 +6 38 37 0 +7 8 39 0 +7 39 38 0 +8 9 40 0 +8 40 39 0 +9 10 41 0 +9 41 40 0 +10 11 42 0 +10 42 41 0 +11 12 43 0 +11 43 42 0 +12 13 44 0 +12 44 43 0 +13 14 45 0 +13 45 44 0 +14 15 46 0 +14 46 45 0 +15 16 47 0 +15 47 46 0 +16 17 48 0 +16 48 47 0 +17 18 49 0 +17 49 48 0 +18 19 50 0 +18 50 49 0 +19 20 51 0 +19 51 50 0 +20 21 52 0 +20 52 51 0 +21 22 53 0 +21 53 52 0 +22 23 54 0 +22 54 53 0 +23 24 55 0 +23 55 54 0 +24 25 56 0 +24 56 55 0 +25 26 57 0 +25 57 56 0 +26 27 58 0 +26 58 57 0 +27 28 59 0 +27 59 58 0 +28 29 60 0 +28 60 59 0 +29 30 61 0 +29 61 60 0 +30 31 62 0 +30 62 61 0 +32 33 64 0 +32 64 63 0 +33 34 65 0 +33 65 64 0 +34 35 66 0 +34 66 65 0 +35 36 67 0 +35 67 66 0 +36 37 68 0 +36 68 67 0 +37 38 69 0 +37 69 68 0 +38 39 70 0 +38 70 69 0 +39 40 71 0 +39 71 70 0 +40 41 72 0 +40 72 71 0 +41 42 73 0 +41 73 72 0 +42 43 74 0 +42 74 73 0 +43 44 75 0 +43 75 74 0 +44 45 76 0 +44 76 75 0 +45 46 77 0 +45 77 76 0 +46 47 78 0 +46 78 77 0 +47 48 79 0 +47 79 78 0 +48 49 80 0 +48 80 79 0 +49 50 81 0 +49 81 80 0 +50 51 82 0 +50 82 81 0 +51 52 83 0 +51 83 82 0 +52 53 84 0 +52 84 83 0 +53 54 85 0 +53 85 84 0 +54 55 86 0 +54 86 85 0 +55 56 87 0 +55 87 86 0 +56 57 88 0 +56 88 87 0 +57 58 89 0 +57 89 88 0 +58 59 90 0 +58 90 89 0 +59 60 91 0 +59 91 90 0 +60 61 92 0 +60 92 91 0 +61 62 93 0 +61 93 92 0 +63 64 95 0 +63 95 94 0 +64 65 96 0 +64 96 95 0 +65 66 97 0 +65 97 96 0 +66 67 98 0 +66 98 97 0 +67 68 99 0 +67 99 98 0 +68 69 100 0 +68 100 99 0 +69 70 101 0 +69 101 100 0 +70 71 102 0 +70 102 101 0 +71 72 103 0 +71 103 102 0 +72 73 104 0 +72 104 103 0 +73 74 105 0 +73 105 104 0 +74 75 106 0 +74 106 105 0 +75 76 107 0 +75 107 106 0 +76 77 108 0 +76 108 107 0 +77 78 109 0 +77 109 108 0 +78 79 110 0 +78 110 109 0 +79 80 111 0 +79 111 110 0 +80 81 112 0 +80 112 111 0 +81 82 113 0 +81 113 112 0 +82 83 114 0 +82 114 113 0 +83 84 115 0 +83 115 114 0 +84 85 116 0 +84 116 115 0 +85 86 117 0 +85 117 116 0 +86 87 118 0 +86 118 117 0 +87 88 119 0 +87 119 118 0 +88 89 120 0 +88 120 119 0 +89 90 121 0 +89 121 120 0 +90 91 122 0 +90 122 121 0 +91 92 123 0 +91 123 122 0 +92 93 124 0 +92 124 123 0 +94 95 126 0 +94 126 125 0 +95 96 127 0 +95 127 126 0 +96 97 128 0 +96 128 127 0 +97 98 129 0 +97 129 128 0 +98 99 130 0 +98 130 129 0 +99 100 131 0 +99 131 130 0 +100 101 132 0 +100 132 131 0 +101 102 133 0 +101 133 132 0 +102 103 134 0 +102 134 133 0 +103 104 135 0 +103 135 134 0 +104 105 136 0 +104 136 135 0 +105 106 137 0 +105 137 136 0 +106 107 138 0 +106 138 137 0 +107 108 139 0 +107 139 138 0 +108 109 140 0 +108 140 139 0 +109 110 141 0 +109 141 140 0 +110 111 142 0 +110 142 141 0 +111 112 143 0 +111 143 142 0 +112 113 144 0 +112 144 143 0 +113 114 145 0 +113 145 144 0 +114 115 146 0 +114 146 145 0 +115 116 147 0 +115 147 146 0 +116 117 148 0 +116 148 147 0 +117 118 149 0 +117 149 148 0 +118 119 150 0 +118 150 149 0 +119 120 151 0 +119 151 150 0 +120 121 152 0 +120 152 151 0 +121 122 153 0 +121 153 152 0 +122 123 154 0 +122 154 153 0 +123 124 155 0 +123 155 154 0 +125 126 157 0 +125 157 156 0 +126 127 158 0 +126 158 157 0 +127 128 159 0 +127 159 158 0 +128 129 160 0 +128 160 159 0 +129 130 161 0 +129 161 160 0 +130 131 162 0 +130 162 161 0 +131 132 163 0 +131 163 162 0 +132 133 164 0 +132 164 163 0 +133 134 165 0 +133 165 164 0 +134 135 166 0 +134 166 165 0 +135 136 167 0 +135 167 166 0 +136 137 168 0 +136 168 167 0 +137 138 169 0 +137 169 168 0 +138 139 170 0 +138 170 169 0 +139 140 171 0 +139 171 170 0 +140 141 172 0 +140 172 171 0 +141 142 173 0 +141 173 172 0 +142 143 174 0 +142 174 173 0 +143 144 175 0 +143 175 174 0 +144 145 176 0 +144 176 175 0 +145 146 177 0 +145 177 176 0 +146 147 178 0 +146 178 177 0 +147 148 179 0 +147 179 178 0 +148 149 180 0 +148 180 179 0 +149 150 181 0 +149 181 180 0 +150 151 182 0 +150 182 181 0 +151 152 183 0 +151 183 182 0 +152 153 184 0 +152 184 183 0 +153 154 185 0 +153 185 184 0 +154 155 186 0 +154 186 185 0 +End diff --git a/Data/Mesh/nonlinearmembrane_square.mesh b/Data/Mesh/nonlinearmembrane_square.mesh new file mode 100644 index 0000000000000000000000000000000000000000..74d343fe459d2bc4a1555b0c536ac476a32b8a89 --- /dev/null +++ b/Data/Mesh/nonlinearmembrane_square.mesh @@ -0,0 +1,308 @@ +MeshVersionFormatted 2 +Dimension +3 +Vertices +100 +0 0 0 2 +1 0 0 2 +1 1 0 2 +0 1 0 2 +0.11111111111091 0 0 2 +0.22222222222171 0 0 2 +0.3333333333325 0 0 2 +0.44444444444329 0 0 2 +0.55555555555438 0 0 2 +0.66666666666579 0 0 2 +0.77777777777719 0 0 2 +0.8888888888886 0 0 2 +1 0.11111111111091 0 2 +1 0.22222222222171 0 2 +1 0.3333333333325 0 2 +1 0.44444444444329 0 2 +1 0.55555555555438 0 2 +1 0.66666666666579 0 2 +1 0.77777777777719 0 2 +1 0.8888888888886 0 2 +0.88888888888843 1 0 2 +0.77777777777793 1 0 2 +0.66666666666759 1 0 2 +0.55555555555725 1 0 2 +0.44444444444629 1 0 2 +0.33333333333472 1 0 2 +0.22222222222315 1 0 2 +0.11111111111157 1 0 2 +0 0.88888888888843 0 2 +0 0.77777777777793 0 2 +0 0.66666666666759 0 2 +0 0.55555555555725 0 2 +0 0.44444444444629 0 2 +0 0.33333333333472 0 2 +0 0.22222222222315 0 2 +0 0.11111111111157 0 2 +0.11111111111098 0.1111111111115 0 1 +0.11111111111106 0.22222222222299 0 1 +0.11111111111113 0.33333333333447 0 1 +0.1111111111112 0.44444444444596 0 1 +0.11111111111128 0.55555555555693 0 1 +0.11111111111135 0.66666666666739 0 1 +0.11111111111143 0.77777777777785 0 1 +0.1111111111115 0.88888888888845 0 1 +0.22222222222187 0.11111111111143 0 1 +0.22222222222203 0.22222222222283 0 1 +0.22222222222219 0.33333333333423 0 1 +0.22222222222235 0.44444444444563 0 1 +0.22222222222251 0.55555555555661 0 1 +0.22222222222267 0.66666666666719 0 1 +0.22222222222283 0.77777777777777 0 1 +0.22222222222299 0.88888888888846 0 1 +0.33333333333275 0.11111111111135 0 1 +0.33333333333299 0.22222222222267 0 1 +0.33333333333324 0.33333333333398 0 1 +0.33333333333349 0.44444444444529 0 1 +0.33333333333373 0.5555555555563 0 1 +0.33333333333398 0.66666666666699 0 1 +0.33333333333423 0.77777777777769 0 1 +0.33333333333447 0.88888888888848 0 1 +0.44444444444362 0.11111111111128 0 1 +0.44444444444396 0.22222222222251 0 1 +0.44444444444429 0.33333333333373 0 1 +0.44444444444462 0.44444444444496 0 1 +0.44444444444496 0.55555555555598 0 1 +0.44444444444529 0.66666666666679 0 1 +0.44444444444563 0.7777777777776 0 1 +0.44444444444596 0.8888888888885 0 1 +0.5555555555547 0.1111111111112 0 1 +0.55555555555502 0.22222222222235 0 1 +0.55555555555534 0.33333333333349 0 1 +0.55555555555566 0.44444444444462 0 1 +0.55555555555598 0.55555555555566 0 1 +0.5555555555563 0.66666666666659 0 1 +0.55555555555661 0.77777777777752 0 1 +0.55555555555693 0.88888888888852 0 1 +0.66666666666599 0.11111111111113 0 1 +0.66666666666619 0.22222222222219 0 1 +0.66666666666639 0.33333333333324 0 1 +0.66666666666659 0.44444444444429 0 1 +0.66666666666679 0.55555555555534 0 1 +0.66666666666699 0.66666666666639 0 1 +0.66666666666719 0.77777777777744 0 1 +0.66666666666739 0.88888888888854 0 1 +0.77777777777727 0.11111111111106 0 1 +0.77777777777736 0.22222222222203 0 1 +0.77777777777744 0.33333333333299 0 1 +0.77777777777752 0.44444444444396 0 1 +0.7777777777776 0.55555555555502 0 1 +0.77777777777769 0.66666666666619 0 1 +0.77777777777777 0.77777777777736 0 1 +0.77777777777785 0.88888888888856 0 1 +0.88888888888858 0.11111111111098 0 1 +0.88888888888856 0.22222222222187 0 1 +0.88888888888854 0.33333333333275 0 1 +0.88888888888852 0.44444444444362 0 1 +0.8888888888885 0.5555555555547 0 1 +0.88888888888848 0.66666666666599 0 1 +0.88888888888846 0.77777777777727 0 1 +0.88888888888845 0.88888888888858 0 1 +Edges +36 +1 5 2 +5 6 2 +6 7 2 +7 8 2 +8 9 2 +9 10 2 +10 11 2 +11 12 2 +12 2 2 +2 13 2 +13 14 2 +14 15 2 +15 16 2 +16 17 2 +17 18 2 +18 19 2 +19 20 2 +20 3 2 +3 21 2 +21 22 2 +22 23 2 +23 24 2 +24 25 2 +25 26 2 +26 27 2 +27 28 2 +28 4 2 +4 29 2 +29 30 2 +30 31 2 +31 32 2 +32 33 2 +33 34 2 +34 35 2 +35 36 2 +36 1 2 +Triangles +162 +1 5 36 1 +36 5 37 1 +36 37 35 1 +35 37 38 1 +35 38 34 1 +34 38 39 1 +34 39 33 1 +33 39 40 1 +33 40 32 1 +32 40 41 1 +32 41 31 1 +31 41 42 1 +31 42 30 1 +30 42 43 1 +30 43 29 1 +29 43 44 1 +29 44 4 1 +4 44 28 1 +5 6 37 1 +37 6 45 1 +37 45 38 1 +38 45 46 1 +38 46 39 1 +39 46 47 1 +39 47 40 1 +40 47 48 1 +40 48 41 1 +41 48 49 1 +41 49 42 1 +42 49 50 1 +42 50 43 1 +43 50 51 1 +43 51 44 1 +44 51 52 1 +44 52 28 1 +28 52 27 1 +6 7 45 1 +45 7 53 1 +45 53 46 1 +46 53 54 1 +46 54 47 1 +47 54 55 1 +47 55 48 1 +48 55 56 1 +48 56 49 1 +49 56 57 1 +49 57 50 1 +50 57 58 1 +50 58 51 1 +51 58 59 1 +51 59 52 1 +52 59 60 1 +52 60 27 1 +27 60 26 1 +7 8 53 1 +53 8 61 1 +53 61 54 1 +54 61 62 1 +54 62 55 1 +55 62 63 1 +55 63 56 1 +56 63 64 1 +56 64 57 1 +57 64 65 1 +57 65 58 1 +58 65 66 1 +58 66 59 1 +59 66 67 1 +59 67 60 1 +60 67 68 1 +60 68 26 1 +26 68 25 1 +8 9 61 1 +61 9 69 1 +61 69 62 1 +62 69 70 1 +62 70 63 1 +63 70 71 1 +63 71 64 1 +64 71 72 1 +64 72 65 1 +65 72 73 1 +65 73 66 1 +66 73 74 1 +66 74 67 1 +67 74 75 1 +67 75 68 1 +68 75 76 1 +68 76 25 1 +25 76 24 1 +9 10 69 1 +69 10 77 1 +69 77 70 1 +70 77 78 1 +70 78 71 1 +71 78 79 1 +71 79 72 1 +72 79 80 1 +72 80 73 1 +73 80 81 1 +73 81 74 1 +74 81 82 1 +74 82 75 1 +75 82 83 1 +75 83 76 1 +76 83 84 1 +76 84 24 1 +24 84 23 1 +10 11 77 1 +77 11 85 1 +77 85 78 1 +78 85 86 1 +78 86 79 1 +79 86 87 1 +79 87 80 1 +80 87 88 1 +80 88 81 1 +81 88 89 1 +81 89 82 1 +82 89 90 1 +82 90 83 1 +83 90 91 1 +83 91 84 1 +84 91 92 1 +84 92 23 1 +23 92 22 1 +11 12 85 1 +85 12 93 1 +85 93 86 1 +86 93 94 1 +86 94 87 1 +87 94 95 1 +87 95 88 1 +88 95 96 1 +88 96 89 1 +89 96 97 1 +89 97 90 1 +90 97 98 1 +90 98 91 1 +91 98 99 1 +91 99 92 1 +92 99 100 1 +92 100 22 1 +22 100 21 1 +12 2 93 1 +93 2 13 1 +93 13 94 1 +94 13 14 1 +94 14 95 1 +95 14 15 1 +95 15 96 1 +96 15 16 1 +96 16 97 1 +97 16 17 1 +97 17 98 1 +98 17 18 1 +98 18 99 1 +99 18 19 1 +99 19 100 1 +100 19 20 1 +100 20 21 1 +21 20 3 1 +End diff --git a/Sources/ModelInstances/Elasticity/demo_input_elasticity_copy.lua b/Sources/ModelInstances/Elasticity/demo_input_elasticity_copy.lua new file mode 100644 index 0000000000000000000000000000000000000000..8970926d8352f672ec8783d9d8d039677352fadd --- /dev/null +++ b/Sources/ModelInstances/Elasticity/demo_input_elasticity_copy.lua @@ -0,0 +1,575 @@ +-- Comment lines are introduced by "--". +-- In a section (i.e. within braces), all entries must be separated by a comma. + +-- transient +transient = { + + -- Tells which policy is used to describe time evolution. + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'constant_time_step', 'variable_time_step'}}) + time_evolution_policy = "constant_time_step", + + + -- Time at the beginning of the code (in seconds). + -- Expected format: VALUE + -- Constraint: v >= 0. + init_time = 0., + + -- Time step between two iterations, in seconds. + -- Expected format: VALUE + -- Constraint: v > 0. + timeStep = 0.1, + + -- Minimum time step between two iterations, in seconds. + -- Expected format: VALUE + -- Constraint: v > 0. + minimum_time_step = 0.1, + + -- Maximum time, if set to zero run a static case. + -- Expected format: VALUE + -- Constraint: v >= 0. + timeMax = 0. +} + + +-- NumberingSubset1 +NumberingSubset1 = { + -- Name of the numbering subset (not really used; at the moment I just need one input parameter to ground + -- the possible values to choose elsewhere). + -- Expected format: "VALUE" + name = "monolithic", + + -- Whether a vector defined on this numbering subset might be used to compute a movemesh. If true, a + -- FEltSpace featuring this numbering subset will compute additional quantities to enable fast computation. + -- This should be false for most numbering subsets, and when it's true the sole unknown involved should be a + -- displacement. + -- Expected format: 'true' or 'false' (without the quote) + do_move_mesh = false + +} + + + + +-- Unknown1: displacement. +Unknown1 = { + -- Name of the unknown (used for displays in output). + -- Expected format: "VALUE" + name = "solid_displacement", + + -- Index of the god of dof into which the finite element space is defined. + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'scalar', 'vectorial'}) + nature = "vectorial" +} + + + + + +-- Mesh1 +Mesh1 = { + -- Path of the mesh file to use. + -- Expected format: "VALUE" + mesh = "${MOREFEM_ROOT}/Data/Mesh/elasticity_Nx50_Ny20_force_label.mesh", + + -- Format of the input mesh. + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'Ensight', 'Medit'}) + format = "Medit", + + -- Highest dimension of the input mesh. This dimension might be lower than the one effectively read in the mesh file; + -- in which case Coords will be reduced provided all the dropped values are 0. If not, an exception is thrown. + -- Expected format: VALUE + -- Constraint: v <= 3 and v > 0 + dimension = 2, + + -- Space unit of the mesh. + -- Expected format: VALUE + space_unit = 1. +} + + +-- Domain1 - The 2D elements of the bar. +Domain1 = { + -- Index of the geometric mesh upon which the domain is defined (as defined in the present file). Might be + -- left empty if domain not limited to one mesh; at most one value is expected here. + -- Expected format: {VALUE1} + mesh_index = { 1 }, + + -- List of dimensions encompassed by the domain. Might be left empty if no restriction at all upon + -- dimensions. + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: value_in(v, {0, 1, 2, 3}) + dimension_list = { 2 }, + + -- List of mesh labels encompassed by the domain. Might be left empty if no restriction at all upon mesh + -- labels. This parameter does not make sense if no mesh is defined for the domain. + -- Expected format: {VALUE1, VALUE2, ...} + mesh_label_list = { }, + + -- List of geometric element types considered in the domain. Might be left empty if no restriction upon + -- these. + -- Expected format: {"VALUE1", "VALUE2", ...} + geometric_element_type_list = { } +} + + +-- Domain2 - The border upon which surfacic forces are applied. +Domain2 = { + -- Index of the geometric mesh upon which the domain is defined (as defined in the present file). Might be + -- left empty if domain not limited to one mesh; at most one value is expected here. + -- Expected format: {VALUE1} + -- Expected format: VALUE + mesh_index = { 1 }, + + -- List of dimensions encompassed by the domain. Might be left empty if no restriction at all upon + -- dimensions. + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: value_in(v, {0, 1, 2, 3}) + dimension_list = { 1 }, + + -- List of mesh labels encompassed by the domain. Might be left empty if no restriction at all upon mesh + -- labels. This parameter does not make sense if no mesh is defined for the domain. + -- Expected format: {VALUE1, VALUE2, ...} + mesh_label_list = { 2 }, + + -- List of geometric element types considered in the domain. Might be left empty if no restriction upon + -- these. + -- Expected format: {"VALUE1", "VALUE2", ...} + geometric_element_type_list = { } +} + + +-- Domain3 - upon which Dirichlet boundary condition is applied. +Domain3 = { + -- Index of the geometric mesh upon which the domain is defined (as defined in the present file). Might be + -- left empty if domain not limited to one mesh; at most one value is expected here. + -- Expected format: {VALUE1} + -- Expected format: VALUE + mesh_index = { 1 }, + + -- List of dimensions encompassed by the domain. Might be left empty if no restriction at all upon + -- dimensions. + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: value_in(v, {0, 1, 2, 3}) + dimension_list = { }, + + -- List of mesh labels encompassed by the domain. Might be left empty if no restriction at all upon mesh + -- labels. This parameter does not make sense if no mesh is defined for the domain. + -- Expected format: {VALUE1, VALUE2, ...} + mesh_label_list = { 1 }, + + -- List of geometric element types considered in the domain. Might be left empty if no restriction upon + -- these. + -- Expected format: {"VALUE1", "VALUE2", ...} + geometric_element_type_list = { } +} + + +Domain4 = { + -- Index of the geometric mesh upon which the domain is defined (as defined in the present file). Might be + -- left empty if domain not limited to one mesh; at most one value is expected here. + -- Expected format: {VALUE1} + -- Expected format: VALUE + mesh_index = { 1 }, + + -- List of dimensions encompassed by the domain. Might be left empty if no restriction at all upon + -- dimensions. + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: value_in(v, {0, 1, 2, 3}) + dimension_list = { }, + + -- List of mesh labels encompassed by the domain. Might be left empty if no restriction at all upon mesh + -- labels. This parameter does not make sense if no mesh is defined for the domain. + -- Expected format: {VALUE1, VALUE2, ...} + mesh_label_list = { }, + + -- List of geometric element types considered in the domain. Might be left empty if no restriction upon + -- these. + -- Expected format: {"VALUE1", "VALUE2", ...} + geometric_element_type_list = { } +} + + + +EssentialBoundaryCondition1 = { + + -- Name of the boundary condition (must be unique). + -- Expected format: "VALUE" + name = "sole", + + -- Comp1, Comp2 or Comp3 + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'Comp1', 'Comp2', 'Comp3', 'Comp12', 'Comp23', 'Comp13', 'Comp123'}) + component = 'Comp12', + + -- Name of the unknown addressed by the boundary condition. + -- Expected format: "VALUE" + unknown = 'solid_displacement', + + -- Values at each of the relevant component. + -- Expected format: {VALUE1, VALUE2, ...} + value = { 0., 0. }, + + -- Index of the domain onto which essential boundary condition is defined. + -- Expected format: VALUE + domain_index = 3, + + -- Whether the values of the boundary condition may vary over time. + -- Expected format: 'true' or 'false' (without the quote) + is_mutable = false, + + -- Whether a dof of this boundary condition may also belong to another one. This highlights an ill-defined + -- model in most cases, but I nonetheless need it for FSI/ALE. + -- Expected format: 'true' or 'false' (without the quote) + may_overlap = false + +} -- EssentialBoundaryCondition1 + + + +-- FiniteElementSpace1 +FiniteElementSpace1 = { + -- Index of the god of dof into which the finite element space is defined. + -- Expected format: VALUE + god_of_dof_index = 1, + + -- Index of the domain onto which the finite element space is defined. This domain must be unidimensional. + -- Expected format: VALUE + domain_index = 1, + + -- List of all unknowns defined in the finite element space. Unknowns here must be defined in this file as + -- an 'Unknown' block; expected name/identifier is the name given there. + -- Expected format: {"VALUE1", "VALUE2", ...} + unknown_list = {"solid_displacement"}, + + + -- List of the shape function to use for each unknown; + -- Expected format: {"VALUE1", "VALUE2", ...} + shape_function_list = {'P1b'}, + + -- List of the numbering subset to use for each unknown; + -- Expected format: {VALUE1, VALUE2, ...} + numbering_subset_list = { 1 } +} + +-- FiniteElementSpace2 +FiniteElementSpace2 = { + -- Index of the god of dof into which the finite element space is defined. + -- Expected format: VALUE + god_of_dof_index = 1, + + -- Index of the domain onto which the finite element space is defined. This domain must be unidimensional. + -- Expected format: VALUE + domain_index = 2, + + -- List of all unknowns defined in the finite element space. Unknowns here must be defined in this file as + -- an 'Unknown' block; expected name/identifier is the name given there. + -- Expected format: {"VALUE1", "VALUE2", ...} + unknown_list = {'solid_displacement'}, + + -- List of the shape function to use for each unknown; + -- Expected format: {"VALUE1", "VALUE2", ...} + shape_function_list = {'P1'}, + + -- List of the numbering subset to use for each unknown; + -- Expected format: {VALUE1, VALUE2, ...} + numbering_subset_list = { 1 } +} + + + +-- Petsc +Petsc1 = { + -- Absolute tolerance + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: v > 0. + absoluteTolerance = 1e-10, + + -- gmresStart + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: v >= 0 + gmresRestart = 200, + + -- Maximum iteration + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: v > 0 + maxIteration = 1000, + + -- List of preconditioner: { none jacobi sor lu bjacobi ilu asm cholesky }. + -- To use mumps: preconditioner = lu + -- Expected format: {"VALUE1", "VALUE2", ...} + -- Constraint: value_in(v, 'lu') + preconditioner = 'lu', + + -- Relative tolerance + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: v > 0. + relativeTolerance = 1e-6, + + -- Step size tolerance + -- Expected format: {VALUE1, VALUE2, ...} + -- Constraint: v > 0. + stepSizeTolerance = 1e-8, + + -- List of solver: { chebychev cg gmres preonly bicg python }; + -- To use Mumps choose preonly. + -- Expected format: {"VALUE1", "VALUE2", ...} + -- Constraint: value_in(v, {'Mumps', 'Umfpack'}) + solver = 'Mumps' +} + +Solid = { + + -- For 2D operators, which approximation to use. + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'irrelevant', 'plane_strain', 'plane_stress'}) + PlaneStressStrain = "plane_strain", + + VolumicMass = { + + -- How is given the parameter (as a constant, as a Lua function, per quadrature point, etc...) + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'ignore', 'constant', 'lua_function', 'at_quadrature_point', 'piecewise_constant_by_domain'}) + nature = "constant", + + -- Value of the volumic mass in the case nature is 'constant'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: VALUE + -- Constraint: v > 0. + scalar_value = 1., + + -- Value of the volumic mass in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function = none, + + -- Domain indices of the parameter in the case nature is 'piecewise_constant_by_domain'. The various + -- domains given here must not intersect. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_id = { }, + + -- Value of the parameter in the case nature is 'piecewise_constant_by_domain'. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_value = { } + + }, -- VolumicMass + + + + + PoissonRatio = { + + -- How is given the parameter (as a constant, as a Lua function, per quadrature point, etc...) + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'ignore', 'constant', 'lua_function', 'at_quadrature_point', 'piecewise_constant_by_domain'}) + nature = "constant", + + -- Value of the volumic mass in the case nature is 'constant'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: VALUE + -- Constraint: v > 0. + scalar_value = 0.3, + + -- Value of the volumic mass in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function = none, + + -- Domain indices of the parameter in the case nature is 'piecewise_constant_by_domain'. The various + -- domains given here must not intersect. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_id = {}, + + -- Value of the parameter in the case nature is 'piecewise_constant_by_domain'. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_value = {}, + + }, -- PoissonRatio + + YoungModulus = { + + -- How is given the parameter (as a constant, as a Lua function, per quadrature point, etc...) + -- Expected format: "VALUE" + -- Constraint: value_in(v, {'ignore', 'constant', 'lua_function', 'at_quadrature_point', 'piecewise_constant_by_domain'}) + nature = 'constant', + + -- Value of the volumic mass in the case nature is 'constant'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: VALUE + -- Constraint: v > 0. + scalar_value = 3.e7, + + -- Value of the volumic mass in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function = none, + + -- Domain indices of the parameter in the case nature is 'piecewise_constant_by_domain'. The various + -- domains given here must not intersect. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_id = {}, + + -- Value of the parameter in the case nature is 'piecewise_constant_by_domain'. + -- Expected format: {VALUE1, VALUE2, ...} + piecewise_constant_domain_value = {}, + + }, -- YoungModulus + +} -- Solid + +TransientSource1 = { + + -- How is given the transient source value (as a constant, as a Lua function, per quadrature point, etc...) + -- Expected format: {"VALUE1", "VALUE2", ...} + -- Constraint: value_in(v, {'ignore', 'constant', 'lua_function'}) + nature = {"ignore", "ignore", "ignore"}, + + -- Value of the transient source in the case nature is 'constant'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: {VALUE1, VALUE2, ...} + scalar_value = { 0. , -0.0001, 0.}, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_x = none, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_y = none, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_z = none + +} -- TransientSource1 + +-- Surfacic source +TransientSource2 = { + + -- How is given the transient source value (as a constant, as a Lua function, per quadrature point, etc...) + -- Expected format: {"VALUE1", "VALUE2", ...} + -- Constraint: value_in(v, {'ignore', 'constant', 'lua_function'}) + nature = {"constant", "constant", "constant"}, + + -- Value of the transient source in the case nature is 'constant'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Expected format: {VALUE1, VALUE2, ...} + scalar_value = { 1.e6, 0., 0.}, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_x = none, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_y = function (x, y, z) + return -0.5 * math.cos(x) * math.sin(y + .1); + end, + + -- Value of the transient source in the case nature is 'lua_function'(and also initial value if nature is + -- 'at_quadrature_point'; irrelevant otherwise). + -- Function expects as arguments global coordinates (coordinates in the mesh). + -- Expected format: Function in Lua language, for instance: + -- -- function(arg1, arg2, arg3) + -- return arg1 + arg2 - + -- arg3 + -- end + -- sin, cos and tan require a 'math.' + -- preffix. + -- If you do not wish to provide one, put anything you want (e.g. 'none'): the + -- content is not interpreted by LuaOptionFile until an actual use of the underlying function. + lua_function_z = none, + +} -- TransientSource2 + + +-- Result +Result = { + -- Directory in which all the results will be written. This path may use the environment variable + -- MOREFEM_RESULT_DIR, which is either provided in user's environment or automatically set to + -- '/Volumes/Data/${USER}/MoReFEM/Results' in MoReFEM initialization step. + -- Expected format: "VALUE" + -- output_directory = "/home/sebastien/Data/MoReFEM/Results/Hyperelasticity", -- LINUX + output_directory = "/Users/antoineo/Codes/MoReFEM/Results/Elasticity", + + -- Enables to skip some printing in the console. Can be used to WriteSolution every n time. + -- Expected format: VALUE + display_value = 1 +} + diff --git a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.cpp b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.cpp index 9a59daac022fde605418e05595f104fbdb17c374..7ff2b1c307d989cf408949af344922d55b8067cb 100644 --- a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.cpp +++ b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.cpp @@ -36,13 +36,15 @@ namespace MoReFEM elementary_data_type&& a_elementary_data, const scalar_parameter& youngs_modulus, const scalar_parameter& poisson_ratio, - const scalar_parameter& thickness) + const scalar_parameter& thickness, + const scalar_parameter& pretension) : NonlinearLocalVariationalOperator(a_unknown_storage, a_test_unknown_storage, std::move(a_elementary_data)), matrix_parent(), vector_parent(), youngs_modulus_(youngs_modulus), poisson_ratio_(poisson_ratio), - thickness_(thickness) + thickness_(thickness), + pretension_(pretension) { const auto& elementary_data = GetElementaryData(); @@ -75,7 +77,7 @@ namespace MoReFEM { 6, 6 }, // tangent_matrix { 3, 2 }, // displacement_gradient { 3, 2 }, // covariant_basis, - // { 3, 2 }, // contravariant_basis, + { 3, 2 }, // contravariant_basis, { 2, 3 }, // transposed_covariant_basis, { 2, 2 }, // covariant_metric_tensor, { 2, 2 }, // contravariant_metric_tensor, @@ -83,6 +85,9 @@ namespace MoReFEM { 2, Nnode }, //transposed_dphi { Nnode_test, 2 }, // dphi_test_mult_gradient_based_block { Nnode_test, Nnode }, // block_contribution + { 3, 3 }, // invert_generalized_covariant_basis + { 3, 3 }, // test_pk_in_ref_basis + { 3, 3 }, // test_gl_in_ref_basis }}); vector_parent::InitLocalVectorStorage @@ -129,7 +134,8 @@ namespace MoReFEM decltype(auto) youngs_modulus = GetYoungsModulus(); decltype(auto) poisson_ratio = GetPoissonRatio(); - decltype(auto) thickness = GetThickness(); + decltype(auto) thickness = GetThickness(); + decltype(auto) pretension = GetPretension(); const auto& geom_elt = elementary_data.GetCurrentGeomElt(); @@ -154,10 +160,11 @@ namespace MoReFEM const double youngs_modulus_value = youngs_modulus.GetValue(quad_pt, geom_elt); const double poisson_ratio_value = poisson_ratio.GetValue(quad_pt, geom_elt); const double thickness_value = thickness.GetValue(quad_pt, geom_elt); + const auto pretension_value = pretension.GetValue(quad_pt, geom_elt); ComputeContravariantBasis(infos_at_quad_pt, determinant); ComputeDisplacementGradient(infos_at_quad_pt); - ComputeGreenLagrange(); + ComputeGreenLagrange(pretension_value); ComputeTangentTensor(youngs_modulus_value, poisson_ratio_value); ComputeSecondPiolaKirchhoff(); ComputeDe(); @@ -169,7 +176,7 @@ namespace MoReFEM const auto& dphi_test = grad_felt_phi_test; - assert(dphi.GetM() == Nnode); + assert(dphi.GetM() == Nnode); assert(dphi_test.GetM() == Nnode_test); if (parent::DoAssembleIntoMatrix()) @@ -227,6 +234,8 @@ namespace MoReFEM } } } + +// std::cout << "matrix_result" << std::endl << matrix_result << std::endl; } const auto int_Nsurface_comp = static_cast(Nsurface_comp); @@ -264,7 +273,8 @@ namespace MoReFEM } ++quad_pt_index; - } // loop over quadrature points + } // loop over quadrature points + } @@ -304,7 +314,7 @@ namespace MoReFEM } } } - +// std::cout << "covariant basis" << std::endl << covariant_basis << std::endl; auto& covariant_metric_tensor = GetLocalMatrix(); @@ -319,8 +329,8 @@ namespace MoReFEM determinant); // #1259: Actually never used!!! Only determinant seems relevant. - // auto& contravariant_basis = GetLocalMatrix(); -// Seldon::Mlt(1., covariant_basis, contravariant_metric_tensor, contravariant_basis); + auto& contravariant_basis = GetLocalMatrix(); + Seldon::Mlt(1., covariant_basis, contravariant_metric_tensor, contravariant_basis); } @@ -331,13 +341,13 @@ namespace MoReFEM displacement_gradient.Zero(); auto& local_displacement = GetFormerLocalDisplacement(); - - const auto& dphi_geo = infos_at_quad_pt.GetGradientRefGeometricPhi(); + + const auto& dphi = infos_at_quad_pt.GetGradientRefFEltPhi(); constexpr auto Nsurface_comp = 2ul; constexpr auto euclidean_dimension = 3ul; - const unsigned int Nshape_function = static_cast(dphi_geo.GetM()); + const unsigned int Nshape_function = static_cast(dphi.GetM()); for (unsigned int component_shape_function = 0u; component_shape_function < Nsurface_comp; ++component_shape_function) { @@ -351,20 +361,20 @@ namespace MoReFEM { const int int_coord_index = static_cast(coord_index); - const auto local_displacement_index = - euclidean_dimension * shape_fct_index + coord_index; +// const auto local_displacement_index = +// euclidean_dimension * shape_fct_index + coord_index; + const auto local_displacement_index = Nshape_function * coord_index + shape_fct_index; assert(local_displacement_index < local_displacement.size()); displacement_gradient(int_coord_index, int_component_shape_function) - += local_displacement[local_displacement_index] - * dphi_geo(int_shape_fct_index, int_component_shape_function); + += local_displacement[local_displacement_index] * dphi(int_shape_fct_index, int_component_shape_function); } } } } - void NonlinearMembrane::ComputeGreenLagrange() + void NonlinearMembrane::ComputeGreenLagrange(const double pretension) { auto& green_lagrange = GetLocalVector(); const auto& displacement_gradient = @@ -390,6 +400,18 @@ namespace MoReFEM disp_grad_1 * cov_basis_0 + disp_grad_0 * disp_grad_1; } + + //PRETENSION + + if (pretension != 0.) + { + const auto& covariant_metric_tensor = GetLocalMatrix(); + green_lagrange(0) += pretension * covariant_metric_tensor(0, 0); + green_lagrange(1) += pretension * covariant_metric_tensor(1, 1); + green_lagrange(2) += 2 * pretension * covariant_metric_tensor(0, 1); + + } + } @@ -426,6 +448,8 @@ namespace MoReFEM tangent_tensor(1, 0) = tangent_tensor(0, 1); tangent_tensor(2, 0) = tangent_tensor(0, 2); tangent_tensor(2, 1) = tangent_tensor(1, 2); + +// std::cout << "tt" << std::endl << tangent_tensor << std::endl; } @@ -436,6 +460,8 @@ namespace MoReFEM const auto& tangent_tensor = GetLocalMatrix(); Seldon::Mlt(tangent_tensor, green_lagrange, second_PK); + +// std::cout << "pk" << std::endl << second_PK << std::endl; } @@ -466,6 +492,8 @@ namespace MoReFEM // Transposed Wrappers::Seldon::Transpose(De_membrane, transposed_De_membrane); + +// std::cout << "de" << std::endl << De_membrane << std::endl; } @@ -501,12 +529,14 @@ namespace MoReFEM tangent_matrix(2 * i, 2 * i + 1) += second_pk_2; tangent_matrix(2 * i + 1, 2 * i) += second_pk_2; } +// std::cout << "tm" << std::endl << tangent_matrix << std::endl; } if (parent::DoAssembleIntoVector()) { rhs_part.Zero(); Seldon::Mlt(transposed_De_membrane, second_PK, rhs_part); +// std::cout << "rhs_part" << std::endl << rhs_part << std::endl; } } diff --git a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hpp b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hpp index 7c9ef4e3923a45a6b3af823e147c52b7e5315a49..0ab382781ed8c6bcc76a7fa39c09f795cf7fba2a 100644 --- a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hpp +++ b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hpp @@ -41,7 +41,7 @@ namespace MoReFEM */ class NonlinearMembrane final : public NonlinearLocalVariationalOperator, - public Crtp::LocalMatrixStorage, + public Crtp::LocalMatrixStorage, public Crtp::LocalVectorStorage { @@ -60,7 +60,7 @@ namespace MoReFEM using scalar_parameter = ScalarParameter; //! Alias to the parent that provides LocalMatrixStorage. - using matrix_parent = Crtp::LocalMatrixStorage; + using matrix_parent = Crtp::LocalMatrixStorage; //! Alias to the parent that provides LocalVectorStorage. using vector_parent = Crtp::LocalVectorStorage; @@ -90,7 +90,8 @@ namespace MoReFEM elementary_data_type&& elementary_data, const scalar_parameter& youngs_modulus, const scalar_parameter& poisson_ratio, - const scalar_parameter& thickness); + const scalar_parameter& thickness, + const scalar_parameter& pretension); //! Destructor. ~NonlinearMembrane(); @@ -157,7 +158,7 @@ namespace MoReFEM /*! * \brief Compute Green-Lagrange. */ - void ComputeGreenLagrange(); + void ComputeGreenLagrange(const double pretension); /*! * \brief Compute locally the tangent tensor. @@ -199,6 +200,8 @@ namespace MoReFEM * \return Thickness of the solid. */ const scalar_parameter& GetThickness() const noexcept; + + const scalar_parameter& GetPretension() const noexcept; private: @@ -213,6 +216,8 @@ namespace MoReFEM //! Thickness. const scalar_parameter& thickness_; + + const scalar_parameter& pretension_; ///@} @@ -234,7 +239,7 @@ namespace MoReFEM tangent_matrix, displacement_gradient, covariant_basis, - // contravariant_basis, + contravariant_basis, transposed_covariant_basis, covariant_metric_tensor, contravariant_metric_tensor, @@ -242,6 +247,9 @@ namespace MoReFEM transposed_dphi, dphi_test_mult_gradient_based_block, block_contribution, + invert_generalized_covariant_basis, + test_pk_in_ref_basis, + test_gl_in_ref_basis, // block_matrix1 = 0, // block_matrix2, // block_matrix3, diff --git a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hxx b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hxx index c55ae77335e4def82e073f6e93217e480a7eec96..70e53296e81dce9cfad5038215302dee296e3fda 100644 --- a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hxx +++ b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/Local/NonlinearMembrane.hxx @@ -43,6 +43,12 @@ namespace MoReFEM { return thickness_; } + + inline const NonlinearMembrane::scalar_parameter& NonlinearMembrane + ::GetPretension() const noexcept + { + return pretension_; + } inline const std::vector& diff --git a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.cpp b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.cpp index 4d6c820fd6cc4b181d48c8f3281fadb5c288c1b8..ac590e5ea6e01f35a6c034b02abd9bf59a377cfe 100644 --- a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.cpp +++ b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.cpp @@ -26,6 +26,7 @@ namespace MoReFEM const scalar_parameter& youngs_modulus, const scalar_parameter& poisson_ratio, const scalar_parameter& thickness, + const scalar_parameter& pretension, const QuadratureRulePerTopology* const quadrature_rule_per_topology) : parent(felt_space, @@ -36,7 +37,8 @@ namespace MoReFEM DoComputeProcessorWiseLocal2Global::yes, youngs_modulus, poisson_ratio, - thickness) + thickness, + pretension) { } diff --git a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.hpp b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.hpp index 9155293b0f573284a7380cb64690b7b17ddaef3b..ce56f1c86aefcf4a768e8d376b313a12fafdfd4d 100644 --- a/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.hpp +++ b/Sources/OperatorInstances/VariationalOperator/NonlinearForm/NonlinearMembrane.hpp @@ -97,6 +97,7 @@ namespace MoReFEM const scalar_parameter& youngs_modulus, const scalar_parameter& poisson_ratio, const scalar_parameter& thickness, + const scalar_parameter& pretension, const QuadratureRulePerTopology* const quadrature_rule_per_topology = nullptr); //! Destructor.