From c82f6efd88aa8dfafae5012b1a016bedd22b8466 Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 17:36:34 -0600 Subject: [PATCH 1/6] Changing material property type in Linear spring to be compatible with Isolators. This also helps use the MaterialRealCMMAux Auxkernel to get forces and deformations of the spring. Refs #376. --- include/kernels/StressDivergenceSpring.h | 10 +++--- include/materials/LinearSpring.h | 43 +++++++++++++----------- src/kernels/StressDivergenceSpring.C | 10 +++--- src/materials/LinearSpring.C | 32 ++++++++++++++---- 4 files changed, 58 insertions(+), 37 deletions(-) diff --git a/include/kernels/StressDivergenceSpring.h b/include/kernels/StressDivergenceSpring.h index ddd61f8cd0..9a749f782e 100644 --- a/include/kernels/StressDivergenceSpring.h +++ b/include/kernels/StressDivergenceSpring.h @@ -50,19 +50,19 @@ class StressDivergenceSpring : public Kernel std::vector _rot_var; /// Spring forces - const MaterialProperty & _spring_forces_global; + const MaterialProperty & _spring_forces_global; /// Spring moments - const MaterialProperty & _spring_moments_global; + const MaterialProperty & _spring_moments_global; /// Displacement stiffness matrix - const MaterialProperty & _kdd; + const MaterialProperty & _kdd; /// Rotation stiffness matrix - const MaterialProperty & _krr; + const MaterialProperty & _krr; /// Rotation stiffness matrix - const MaterialProperty & _total_global_to_local_rotation; + const MaterialProperty & _total_global_to_local_rotation; }; #endif // STRESSDIVERGENCESPRING_H diff --git a/include/materials/LinearSpring.h b/include/materials/LinearSpring.h index 5479b4c379..add626353a 100644 --- a/include/materials/LinearSpring.h +++ b/include/materials/LinearSpring.h @@ -29,6 +29,9 @@ class LinearSpring : public Material LinearSpring(const InputParameters & parameters); protected: + /// Initialize stateful properties + virtual void initQpStatefulProperties() override; + /// Compute properties for each qp virtual void computeQpProperties() override; @@ -61,10 +64,10 @@ class LinearSpring : public Material std::vector _disp_num; /// Deformations in the spring calculated at a time t in the spring local coordinates - MaterialProperty & _deformations; + MaterialProperty & _deformations; /// Rotations in the spring calculated at a time t in the spring local coordinates - MaterialProperty & _rotations; + MaterialProperty & _rotations; /// Axial stiffness of the spring (local x direction) const VariableValue & _kx; @@ -85,58 +88,58 @@ class LinearSpring : public Material const VariableValue & _krz; /// Global displacement at node 0 of the spring element - RealVectorValue _global_disp0; // node 0 + ColumnMajorMatrix _global_disp0; // node 0 /// Global displacement at node 1 of the spring element - RealVectorValue _global_disp1; // node 1 + ColumnMajorMatrix _global_disp1; // node 1 /// Global rotation at node 0 of the spring element - RealVectorValue _global_rot0; // node 0 + ColumnMajorMatrix _global_rot0; // node 0 /// Global rotation at node 1 of the spring element - RealVectorValue _global_rot1; // node 1 + ColumnMajorMatrix _global_rot1; // node 1 /// Local displacement at node 0 of the spring element - RealVectorValue _local_disp0; // node 0 + ColumnMajorMatrix _local_disp0; // node 0 /// Local displacement at node 1 of the spring element - RealVectorValue _local_disp1; // node 1 + ColumnMajorMatrix _local_disp1; // node 1 /// Local rotation at node 0 of the spring element - RealVectorValue _local_rot0; // node 0 + ColumnMajorMatrix _local_rot0; // node 0 /// Local rotation at node 1 of the spring element - RealVectorValue _local_rot1; // node 1 + ColumnMajorMatrix _local_rot1; // node 1 /// Spring forces in the local coordinate system - RealVectorValue _spring_forces_local; + ColumnMajorMatrix _spring_forces_local; // Spring moments in the local coordinate system - RealVectorValue _spring_moments_local; + ColumnMajorMatrix _spring_moments_local; /// Spring displacement stiffness matrix in the local coordinate system - RankTwoTensor _kdd_local; + ColumnMajorMatrix _kdd_local; /// Spring rotational stiffness matrix in the local coordinate system - RankTwoTensor _krr_local; + ColumnMajorMatrix _krr_local; /// Spring forces in the global coordinate system - MaterialProperty & _spring_forces_global; + MaterialProperty & _spring_forces_global; /// Spring moments in the global coordinate system - MaterialProperty & _spring_moments_global; + MaterialProperty & _spring_moments_global; /// Spring displacement stiffness matrix in the global coordinate system - MaterialProperty & _kdd; + MaterialProperty & _kdd; /// Spring rotational stiffness matrix in the global coordinate system - MaterialProperty & _krr; + MaterialProperty & _krr; /// Rotational transformation from global coordinate system to spring local configuration at t = 0 - RankTwoTensor _original_global_to_local_rotation; + ColumnMajorMatrix _original_global_to_local_rotation; /// Rotational transformation from global coordinate system to spring local configuration at any time - MaterialProperty & _total_global_to_local_rotation; + MaterialProperty & _total_global_to_local_rotation; }; #endif // LINEARSPRING_H diff --git a/src/kernels/StressDivergenceSpring.C b/src/kernels/StressDivergenceSpring.C index 18b7f554e7..42a2500639 100644 --- a/src/kernels/StressDivergenceSpring.C +++ b/src/kernels/StressDivergenceSpring.C @@ -53,12 +53,12 @@ StressDivergenceSpring::StressDivergenceSpring(const InputParameters & parameter _disp_var(_ndisp), _nrot(coupledComponents("rotations")), _rot_var(_nrot), - _spring_forces_global(getMaterialPropertyByName("global_forces")), - _spring_moments_global(getMaterialPropertyByName("global_moments")), - _kdd(getMaterialPropertyByName("displacement_stiffness_matrix")), - _krr(getMaterialPropertyByName("rotation_stiffness_matrix")), + _spring_forces_global(getMaterialPropertyByName("global_forces")), + _spring_moments_global(getMaterialPropertyByName("global_moments")), + _kdd(getMaterialPropertyByName("displacement_stiffness_matrix")), + _krr(getMaterialPropertyByName("rotation_stiffness_matrix")), _total_global_to_local_rotation( - getMaterialPropertyByName("total_global_to_local_rotation")) + getMaterialPropertyByName("total_global_to_local_rotation")) { if (_component > 5) mooseError("Error in StressDivergenceSpring block ", diff --git a/src/materials/LinearSpring.C b/src/materials/LinearSpring.C index d1b06a4baa..c6e77e5abb 100644 --- a/src/materials/LinearSpring.C +++ b/src/materials/LinearSpring.C @@ -46,20 +46,20 @@ LinearSpring::LinearSpring(const InputParameters & parameters) _ndisp(coupledComponents("displacements")), _rot_num(3), _disp_num(3), - _deformations(declareProperty("deformations")), - _rotations(declareProperty("rotations")), + _deformations(declareProperty("deformations")), + _rotations(declareProperty("rotations")), _kx(coupledValue("kx")), _ky(coupledValue("ky")), _kz(coupledValue("kz")), _krx(coupledValue("krx")), _kry(coupledValue("kry")), _krz(coupledValue("krz")), - _spring_forces_global(declareProperty("global_forces")), - _spring_moments_global(declareProperty("global_moments")), - _kdd(declareProperty("displacement_stiffness_matrix")), - _krr(declareProperty("rotation_stiffness_matrix")), + _spring_forces_global(declareProperty("global_forces")), + _spring_moments_global(declareProperty("global_moments")), + _kdd(declareProperty("displacement_stiffness_matrix")), + _krr(declareProperty("rotation_stiffness_matrix")), _total_global_to_local_rotation( - declareProperty("total_global_to_local_rotation")) + declareProperty("total_global_to_local_rotation")) { // Checking for consistency between length of the provided displacements and rotations vector if (_ndisp != _nrot) @@ -77,6 +77,24 @@ LinearSpring::LinearSpring(const InputParameters & parameters) } } +void +LinearSpring::initQpStatefulProperties() +{ + _deformations[_qp].reshape(3, 1); + _rotations[_qp].reshape(3, 1); + _spring_forces_global[_qp].reshape(3, 1); + _spring_forces_global[_qp].zero(); + _spring_moments_global[_qp].reshape(3, 1); + _kdd[_qp].reshape(3, 3); + _krr[_qp].reshape(3, 3); + _total_global_to_local_rotation[_qp].reshape(3, 3); + _original_global_to_local_rotation.reshape(3, 3); + _global_disp0.reshape(3, 1); + _global_disp1.reshape(3, 1); + _global_rot0.reshape(3, 1); + _global_rot1.reshape(3, 1); +} + void LinearSpring::computeQpProperties() { From 45dd85c67237487853a095b783df414e45d5d204 Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 17:38:20 -0600 Subject: [PATCH 2/6] Added MaterialRealCMMAux auxkernel. Refs #376. --- include/auxkernels/MaterialRealCMMAux.h | 49 +++++++++++++++++++++++++ src/auxkernels/MaterialRealCMMAux.C | 45 +++++++++++++++++++++++ 2 files changed, 94 insertions(+) create mode 100644 include/auxkernels/MaterialRealCMMAux.h create mode 100644 src/auxkernels/MaterialRealCMMAux.C diff --git a/include/auxkernels/MaterialRealCMMAux.h b/include/auxkernels/MaterialRealCMMAux.h new file mode 100644 index 0000000000..cd579212d1 --- /dev/null +++ b/include/auxkernels/MaterialRealCMMAux.h @@ -0,0 +1,49 @@ +/*************************************************/ +/* DO NOT MODIFY THIS HEADER */ +/* */ +/* MASTODON */ +/* */ +/* (c) 2015 Battelle Energy Alliance, LLC */ +/* ALL RIGHTS RESERVED */ +/* */ +/* Prepared by Battelle Energy Alliance, LLC */ +/* With the U. S. Department of Energy */ +/* */ +/* See COPYRIGHT for full restrictions */ +/*************************************************/ + +#pragma once + +// MOOSE includes +#include "MaterialAuxBase.h" + +// Forward declarations +class MaterialRealCMMAux; + +template <> +InputParameters validParams(); + +/** + * AuxKernel for outputting a DenseMatrix material property component to an AuxVariable + */ +class MaterialRealCMMAux : public MaterialAuxBase +{ +public: + static InputParameters validParams(); + + /** + * Class constructor + * @param parameters The input parameters for this AuxKernel + */ + MaterialRealCMMAux(const InputParameters & parameters); + +protected: + /// Returns the component of the tensor for output + virtual Real getRealValue() override; + + /// The row index to output + unsigned int _row; + + /// The column index to output + unsigned int _col; +}; diff --git a/src/auxkernels/MaterialRealCMMAux.C b/src/auxkernels/MaterialRealCMMAux.C new file mode 100644 index 0000000000..ee1e655677 --- /dev/null +++ b/src/auxkernels/MaterialRealCMMAux.C @@ -0,0 +1,45 @@ +/*************************************************/ +/* DO NOT MODIFY THIS HEADER */ +/* */ +/* MASTODON */ +/* */ +/* (c) 2015 Battelle Energy Alliance, LLC */ +/* ALL RIGHTS RESERVED */ +/* */ +/* Prepared by Battelle Energy Alliance, LLC */ +/* With the U. S. Department of Energy */ +/* */ +/* See COPYRIGHT for full restrictions */ +/*************************************************/ + +#include "MaterialRealCMMAux.h" + +registerMooseObject("MastodonApp", MaterialRealCMMAux); + +defineLegacyParams(MaterialRealCMMAux); + +InputParameters +MaterialRealCMMAux::validParams() +{ + InputParameters params = MaterialAuxBase<>::validParams(); + params.addClassDescription( + "Populate an auxiliary variable with an entry from a ColumnMajorMatrix material property."); + params.addParam("row", 0, "The row component to consider for this kernel"); + params.addParam("column", 0, "The column component to consider for this kernel"); + return params; +} + +MaterialRealCMMAux::MaterialRealCMMAux(const InputParameters & parameters) + : MaterialAuxBase(parameters), + _row(getParam("row")), + _col(getParam("column")) +{ +} + +Real +MaterialRealCMMAux::getRealValue() +{ + if (_row >= _prop[_qp].n() || _col >= _prop[_qp].m()) + mooseError("In ", name(), " the row or column parameter is out of bounds."); + return _prop[_qp](_row, _col); +} From 7748e8cf3e6426db7d3c5c883cad676956ab5240 Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 17:40:02 -0600 Subject: [PATCH 3/6] Adding tests for MateralRealCMMAux. Refs #376. --- .../gold/spring_static_out.csv | 802 ++++++++++++++++++ .../materialrealcmm/spring_static.i | 303 +++++++ test/tests/auxkernels/materialrealcmm/tests | 11 + 3 files changed, 1116 insertions(+) create mode 100644 test/tests/auxkernels/materialrealcmm/gold/spring_static_out.csv create mode 100644 test/tests/auxkernels/materialrealcmm/spring_static.i create mode 100644 test/tests/auxkernels/materialrealcmm/tests diff --git a/test/tests/auxkernels/materialrealcmm/gold/spring_static_out.csv b/test/tests/auxkernels/materialrealcmm/gold/spring_static_out.csv new file mode 100644 index 0000000000..af95be5d8d --- /dev/null +++ b/test/tests/auxkernels/materialrealcmm/gold/spring_static_out.csv @@ -0,0 +1,802 @@ +time,disp_x,disp_y,disp_z,gf_x,rot_x,rot_y,rot_z +0,0,0,0,0,0,0,0 +0.005,0.005,0.005,0.005,0.005,0.005,0.005,0.005 +0.01,0.01,0.01,0.01,0.01,0.01,0.01,0.01 +0.015,0.015,0.015,0.015,0.015,0.015,0.015,0.015 +0.02,0.02,0.02,0.02,0.02,0.02,0.02,0.02 +0.025,0.025,0.025,0.025,0.025,0.025,0.025,0.025 +0.03,0.03,0.03,0.03,0.03,0.03,0.03,0.03 +0.035,0.035,0.035,0.035,0.035,0.035,0.035,0.035 +0.04,0.04,0.04,0.04,0.04,0.04,0.04,0.04 +0.045,0.045,0.045,0.045,0.045,0.045,0.045,0.045 +0.05,0.05,0.05,0.05,0.05,0.05,0.05,0.05 +0.055,0.055,0.055,0.055,0.055,0.055,0.055,0.055 +0.06,0.06,0.06,0.06,0.06,0.06,0.06,0.06 +0.065,0.065,0.065,0.065,0.065,0.065,0.065,0.065 +0.07,0.07,0.07,0.07,0.07,0.07,0.07,0.07 +0.075,0.075,0.075,0.075,0.075,0.075,0.075,0.075 +0.08,0.08,0.08,0.08,0.08,0.08,0.08,0.08 +0.085,0.085,0.085,0.085,0.085,0.085,0.085,0.085 +0.09,0.09,0.09,0.09,0.09,0.09,0.09,0.09 +0.095,0.095,0.095,0.095,0.095,0.095,0.095,0.095 +0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1 +0.105,0.105,0.105,0.105,0.105,0.105,0.105,0.105 +0.11,0.11,0.11,0.11,0.11,0.11,0.11,0.11 +0.115,0.115,0.115,0.115,0.115,0.115,0.115,0.115 +0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12 +0.125,0.125,0.125,0.125,0.125,0.125,0.125,0.125 +0.13,0.13,0.13,0.13,0.13,0.13,0.13,0.13 +0.135,0.135,0.135,0.135,0.135,0.135,0.135,0.135 +0.14,0.14,0.14,0.14,0.14,0.14,0.14,0.14 +0.145,0.145,0.145,0.145,0.145,0.145,0.145,0.145 +0.15,0.15,0.15,0.15,0.15,0.15,0.15,0.15 +0.155,0.155,0.155,0.155,0.155,0.155,0.155,0.155 +0.16,0.16,0.16,0.16,0.16,0.16,0.16,0.16 +0.165,0.165,0.165,0.165,0.165,0.165,0.165,0.165 +0.17,0.17,0.17,0.17,0.17,0.17,0.17,0.17 +0.175,0.175,0.175,0.175,0.175,0.175,0.175,0.175 +0.18,0.18,0.18,0.18,0.18,0.18,0.18,0.18 +0.185,0.185,0.185,0.185,0.185,0.185,0.185,0.185 +0.19,0.19,0.19,0.19,0.19,0.19,0.19,0.19 +0.195,0.195,0.195,0.195,0.195,0.195,0.195,0.195 +0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2 +0.205,0.205,0.205,0.205,0.205,0.205,0.205,0.205 +0.21,0.21,0.21,0.21,0.21,0.21,0.21,0.21 +0.215,0.215,0.215,0.215,0.215,0.215,0.215,0.215 +0.22,0.22,0.22,0.22,0.22,0.22,0.22,0.22 +0.225,0.225,0.225,0.225,0.225,0.225,0.225,0.225 +0.23,0.23,0.23,0.23,0.23,0.23,0.23,0.23 +0.235,0.235,0.235,0.235,0.235,0.235,0.235,0.235 +0.24,0.24,0.24,0.24,0.24,0.24,0.24,0.24 +0.245,0.245,0.245,0.245,0.245,0.245,0.245,0.245 +0.25,0.25,0.25,0.25,0.25,0.25,0.25,0.25 +0.255,0.255,0.255,0.255,0.255,0.255,0.255,0.255 +0.26,0.26,0.26,0.26,0.26,0.26,0.26,0.26 +0.265,0.265,0.265,0.265,0.265,0.265,0.265,0.265 +0.27,0.27,0.27,0.27,0.27,0.27,0.27,0.27 +0.275,0.275,0.275,0.275,0.275,0.275,0.275,0.275 +0.28,0.28,0.28,0.28,0.28,0.28,0.28,0.28 +0.285,0.285,0.285,0.285,0.285,0.285,0.285,0.285 +0.29,0.29,0.29,0.29,0.29,0.29,0.29,0.29 +0.295,0.295,0.295,0.295,0.295,0.295,0.295,0.295 +0.3,0.3,0.3,0.3,0.3,0.3,0.3,0.3 +0.305,0.305,0.305,0.305,0.305,0.305,0.305,0.305 +0.31,0.31,0.31,0.31,0.31,0.31,0.31,0.31 +0.315,0.315,0.315,0.315,0.315,0.315,0.315,0.315 +0.32,0.32,0.32,0.32,0.32,0.32,0.32,0.32 +0.325,0.325,0.325,0.325,0.325,0.325,0.325,0.325 +0.33,0.33,0.33,0.33,0.33,0.33,0.33,0.33 +0.335,0.335,0.335,0.335,0.335,0.335,0.335,0.335 +0.34,0.34,0.34,0.34,0.34,0.34,0.34,0.34 +0.345,0.345,0.345,0.345,0.345,0.345,0.345,0.345 +0.35,0.35,0.35,0.35,0.35,0.35,0.35,0.35 +0.355,0.355,0.355,0.355,0.355,0.355,0.355,0.355 +0.36,0.36,0.36,0.36,0.36,0.36,0.36,0.36 +0.365,0.365,0.365,0.365,0.365,0.365,0.365,0.365 +0.37,0.37,0.37,0.37,0.37,0.37,0.37,0.37 +0.375,0.375,0.375,0.375,0.375,0.375,0.375,0.375 +0.38,0.38,0.38,0.38,0.38,0.38,0.38,0.38 +0.385,0.385,0.385,0.385,0.385,0.385,0.385,0.385 +0.39,0.39,0.39,0.39,0.39,0.39,0.39,0.39 +0.395,0.395,0.395,0.395,0.395,0.395,0.395,0.395 +0.4,0.4,0.4,0.4,0.4,0.4,0.4,0.4 +0.405,0.405,0.405,0.405,0.405,0.405,0.405,0.405 +0.41,0.41,0.41,0.41,0.41,0.41,0.41,0.41 +0.415,0.415,0.415,0.415,0.415,0.415,0.415,0.415 +0.42,0.42,0.42,0.42,0.42,0.42,0.42,0.42 +0.425,0.425,0.425,0.425,0.425,0.425,0.425,0.425 +0.43,0.43,0.43,0.43,0.43,0.43,0.43,0.43 +0.435,0.435,0.435,0.435,0.435,0.435,0.435,0.435 +0.44,0.44,0.44,0.44,0.44,0.44,0.44,0.44 +0.445,0.445,0.445,0.445,0.445,0.445,0.445,0.445 +0.45,0.45,0.45,0.45,0.45,0.45,0.45,0.45 +0.455,0.455,0.455,0.455,0.455,0.455,0.455,0.455 +0.46,0.46,0.46,0.46,0.46,0.46,0.46,0.46 +0.465,0.465,0.465,0.465,0.465,0.465,0.465,0.465 +0.47,0.47,0.47,0.47,0.47,0.47,0.47,0.47 +0.475,0.475,0.475,0.475,0.475,0.475,0.475,0.475 +0.48,0.48,0.48,0.48,0.48,0.48,0.48,0.48 +0.485,0.485,0.485,0.485,0.485,0.485,0.485,0.485 +0.49,0.49,0.49,0.49,0.49,0.49,0.49,0.49 +0.495,0.495,0.495,0.495,0.495,0.495,0.495,0.495 +0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5 +0.505,0.505,0.505,0.505,0.505,0.505,0.505,0.505 +0.51,0.51,0.51,0.51,0.51,0.51,0.51,0.51 +0.515,0.515,0.515,0.515,0.515,0.515,0.515,0.515 +0.52,0.52,0.52,0.52,0.52,0.52,0.52,0.52 +0.525,0.525,0.525,0.525,0.525,0.525,0.525,0.525 +0.53,0.53,0.53,0.53,0.53,0.53,0.53,0.53 +0.535,0.535,0.535,0.535,0.535,0.535,0.535,0.535 +0.54,0.54,0.54,0.54,0.54,0.54,0.54,0.54 +0.545,0.545,0.545,0.545,0.545,0.545,0.545,0.545 +0.55,0.55,0.55,0.55,0.55,0.55,0.55,0.55 +0.555,0.555,0.555,0.555,0.555,0.555,0.555,0.555 +0.56,0.56,0.56,0.56,0.56,0.56,0.56,0.56 +0.565,0.565,0.565,0.565,0.565,0.565,0.565,0.565 +0.57,0.57,0.57,0.57,0.57,0.57,0.57,0.57 +0.575,0.575,0.575,0.575,0.575,0.575,0.575,0.575 +0.58,0.58,0.58,0.58,0.58,0.58,0.58,0.58 +0.585,0.585,0.585,0.585,0.585,0.585,0.585,0.585 +0.59,0.59,0.59,0.59,0.59,0.59,0.59,0.59 +0.595,0.595,0.595,0.595,0.595,0.595,0.595,0.595 +0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6 +0.605,0.605,0.605,0.605,0.605,0.605,0.605,0.605 +0.61,0.61,0.61,0.61,0.61,0.61,0.61,0.61 +0.615,0.615,0.615,0.615,0.615,0.615,0.615,0.615 +0.62,0.62,0.62,0.62,0.62,0.62,0.62,0.62 +0.625,0.625,0.625,0.625,0.625,0.625,0.625,0.625 +0.63,0.63,0.63,0.63,0.63,0.63,0.63,0.63 +0.635,0.635,0.635,0.635,0.635,0.635,0.635,0.635 +0.64,0.64,0.64,0.64,0.64,0.64,0.64,0.64 +0.645,0.645,0.645,0.645,0.645,0.645,0.645,0.645 +0.65,0.65,0.65,0.65,0.65,0.65,0.65,0.65 +0.655,0.655,0.655,0.655,0.655,0.655,0.655,0.655 +0.66,0.66,0.66,0.66,0.66,0.66,0.66,0.66 +0.665,0.665,0.665,0.665,0.665,0.665,0.665,0.665 +0.67,0.67,0.67,0.67,0.67,0.67,0.67,0.67 +0.675,0.675,0.675,0.675,0.675,0.675,0.675,0.675 +0.68,0.68,0.68,0.68,0.68,0.68,0.68,0.68 +0.685,0.685,0.685,0.685,0.685,0.685,0.685,0.685 +0.69,0.69,0.69,0.69,0.69,0.69,0.69,0.69 +0.695,0.695,0.695,0.695,0.695,0.695,0.695,0.695 +0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7 +0.705,0.705,0.705,0.705,0.705,0.705,0.705,0.705 +0.71,0.71,0.71,0.71,0.71,0.71,0.71,0.71 +0.715,0.715,0.715,0.715,0.715,0.715,0.715,0.715 +0.72,0.72,0.72,0.72,0.72,0.72,0.72,0.72 +0.725,0.725,0.725,0.725,0.725,0.725,0.725,0.725 +0.73,0.73,0.73,0.73,0.73,0.73,0.73,0.73 +0.735,0.735,0.735,0.735,0.735,0.735,0.735,0.735 +0.74,0.74,0.74,0.74,0.74,0.74,0.74,0.74 +0.745,0.745,0.745,0.745,0.745,0.745,0.745,0.745 +0.75,0.75,0.75,0.75,0.75,0.75,0.75,0.75 +0.755,0.755,0.755,0.755,0.755,0.755,0.755,0.755 +0.76,0.76,0.76,0.76,0.76,0.76,0.76,0.76 +0.765,0.765,0.765,0.765,0.765,0.765,0.765,0.765 +0.77,0.77,0.77,0.77,0.77,0.77,0.77,0.77 +0.775,0.775,0.775,0.775,0.775,0.775,0.775,0.775 +0.78,0.78,0.78,0.78,0.78,0.78,0.78,0.78 +0.785,0.785,0.785,0.785,0.785,0.785,0.785,0.785 +0.79,0.79,0.79,0.79,0.79,0.79,0.79,0.79 +0.795,0.795,0.795,0.795,0.795,0.795,0.795,0.795 +0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8 +0.805,0.805,0.805,0.805,0.805,0.805,0.805,0.805 +0.81,0.81,0.81,0.81,0.81,0.81,0.81,0.81 +0.815,0.815,0.815,0.815,0.815,0.815,0.815,0.815 +0.82,0.82,0.82,0.82,0.82,0.82,0.82,0.82 +0.825,0.825,0.825,0.825,0.825,0.825,0.825,0.825 +0.83,0.83,0.83,0.83,0.83,0.83,0.83,0.83 +0.835,0.835,0.835,0.835,0.835,0.835,0.835,0.835 +0.84,0.84,0.84,0.84,0.84,0.84,0.84,0.84 +0.845,0.845,0.845,0.845,0.845,0.845,0.845,0.845 +0.85,0.85,0.85,0.85,0.85,0.85,0.85,0.85 +0.855,0.855,0.855,0.855,0.855,0.855,0.855,0.855 +0.86,0.86,0.86,0.86,0.86,0.86,0.86,0.86 +0.865,0.865,0.865,0.865,0.865,0.865,0.865,0.865 +0.87,0.87,0.87,0.87,0.87,0.87,0.87,0.87 +0.875,0.875,0.875,0.875,0.875,0.875,0.875,0.875 +0.88,0.88,0.88,0.88,0.88,0.88,0.88,0.88 +0.885,0.885,0.885,0.885,0.885,0.885,0.885,0.885 +0.89,0.89,0.89,0.89,0.89,0.89,0.89,0.89 +0.895,0.895,0.895,0.895,0.895,0.895,0.895,0.895 +0.9,0.9,0.9,0.9,0.9,0.9,0.9,0.9 +0.905,0.905,0.905,0.905,0.905,0.905,0.905,0.905 +0.91,0.91,0.91,0.91,0.91,0.91,0.91,0.91 +0.915,0.915,0.915,0.915,0.915,0.915,0.915,0.915 +0.92,0.92,0.92,0.92,0.92,0.92,0.92,0.92 +0.925,0.925,0.925,0.925,0.925,0.925,0.925,0.925 +0.93,0.93,0.93,0.93,0.93,0.93,0.93,0.93 +0.935,0.935,0.935,0.935,0.935,0.935,0.935,0.935 +0.94,0.94,0.94,0.94,0.94,0.94,0.94,0.94 +0.945,0.945,0.945,0.945,0.945,0.945,0.945,0.945 +0.95,0.95,0.95,0.95,0.95,0.95,0.95,0.95 +0.955,0.955,0.955,0.955,0.955,0.955,0.955,0.955 +0.96,0.96,0.96,0.96,0.96,0.96,0.96,0.96 +0.965,0.965,0.965,0.965,0.965,0.965,0.965,0.965 +0.97,0.97,0.97,0.97,0.97,0.97,0.97,0.97 +0.975,0.975,0.975,0.975,0.975,0.975,0.975,0.975 +0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98 +0.985,0.985,0.985,0.985,0.985,0.985,0.985,0.985 +0.99,0.99,0.99,0.99,0.99,0.99,0.99,0.99 +0.995,0.995,0.995,0.995,0.995,0.995,0.995,0.995 +1,1,1,1,1,1,1,1 +1.005,0.995,0.995,0.995,0.995,0.995,0.995,0.995 +1.01,0.99,0.99,0.99,0.99,0.99,0.99,0.99 +1.015,0.985,0.985,0.985,0.985,0.985,0.985,0.985 +1.02,0.98,0.98,0.98,0.98,0.98,0.98,0.98 +1.025,0.975,0.975,0.975,0.975,0.975,0.975,0.975 +1.03,0.97,0.97,0.97,0.97,0.97,0.97,0.97 +1.035,0.965,0.965,0.965,0.965,0.965,0.965,0.965 +1.04,0.96,0.96,0.96,0.96,0.96,0.96,0.96 +1.045,0.955,0.955,0.955,0.955,0.955,0.955,0.955 +1.05,0.95,0.95,0.95,0.95,0.95,0.95,0.95 +1.055,0.945,0.945,0.945,0.945,0.945,0.945,0.945 +1.06,0.94,0.94,0.94,0.94,0.94,0.94,0.94 +1.065,0.935,0.935,0.935,0.935,0.935,0.935,0.935 +1.07,0.93,0.93,0.93,0.93,0.93,0.93,0.93 +1.075,0.925,0.925,0.925,0.925,0.925,0.925,0.925 +1.08,0.92,0.92,0.92,0.92,0.92,0.92,0.92 +1.085,0.915,0.915,0.915,0.915,0.915,0.915,0.915 +1.09,0.91,0.91,0.91,0.91,0.91,0.91,0.91 +1.095,0.905,0.905,0.905,0.905,0.905,0.905,0.905 +1.1,0.9,0.9,0.9,0.9,0.9,0.9,0.9 +1.105,0.895,0.895,0.895,0.895,0.895,0.895,0.895 +1.11,0.89,0.89,0.89,0.89,0.89,0.89,0.89 +1.115,0.885,0.885,0.885,0.885,0.885,0.885,0.885 +1.12,0.88,0.88,0.88,0.88,0.88,0.88,0.88 +1.125,0.875,0.875,0.875,0.875,0.875,0.875,0.875 +1.13,0.87,0.87,0.87,0.87,0.87,0.87,0.87 +1.135,0.865,0.865,0.865,0.865,0.865,0.865,0.865 +1.14,0.86,0.86,0.86,0.86,0.86,0.86,0.86 +1.145,0.855,0.855,0.855,0.855,0.855,0.855,0.855 +1.15,0.85,0.85,0.85,0.85,0.85,0.85,0.85 +1.155,0.845,0.845,0.845,0.845,0.845,0.845,0.845 +1.16,0.84,0.84,0.84,0.84,0.84,0.84,0.84 +1.165,0.835,0.835,0.835,0.835,0.835,0.835,0.835 +1.17,0.83,0.83,0.83,0.83,0.83,0.83,0.83 +1.175,0.825,0.825,0.825,0.825,0.825,0.825,0.825 +1.18,0.82,0.82,0.82,0.82,0.82,0.82,0.82 +1.185,0.815,0.815,0.815,0.815,0.815,0.815,0.815 +1.19,0.81,0.81,0.81,0.81,0.81,0.81,0.81 +1.195,0.805,0.805,0.805,0.805,0.805,0.805,0.805 +1.2,0.8,0.8,0.8,0.8,0.8,0.8,0.8 +1.205,0.795,0.795,0.795,0.795,0.795,0.795,0.795 +1.21,0.79,0.79,0.79,0.79,0.79,0.79,0.79 +1.215,0.785,0.785,0.785,0.785,0.785,0.785,0.785 +1.22,0.78,0.78,0.78,0.78,0.78,0.78,0.78 +1.225,0.775,0.775,0.775,0.775,0.775,0.775,0.775 +1.23,0.77,0.77,0.77,0.77,0.77,0.77,0.77 +1.235,0.765,0.765,0.765,0.765,0.765,0.765,0.765 +1.24,0.76,0.76,0.76,0.76,0.76,0.76,0.76 +1.245,0.755,0.755,0.755,0.755,0.755,0.755,0.755 +1.25,0.75,0.75,0.75,0.75,0.75,0.75,0.75 +1.255,0.745,0.745,0.745,0.745,0.745,0.745,0.745 +1.26,0.74,0.74,0.74,0.74,0.74,0.74,0.74 +1.265,0.735,0.735,0.73500000000001,0.735,0.735,0.735,0.73500000000001 +1.27,0.73000000000001,0.73000000000001,0.73000000000001,0.73000000000001,0.73000000000001,0.73000000000001,0.73000000000001 +1.275,0.72500000000001,0.72500000000001,0.72500000000001,0.72500000000001,0.72500000000001,0.72500000000001,0.72500000000001 +1.28,0.72000000000001,0.72000000000001,0.72000000000001,0.72000000000001,0.72000000000001,0.72000000000001,0.72000000000001 +1.285,0.71500000000001,0.71500000000001,0.71500000000001,0.71500000000001,0.71500000000001,0.71500000000001,0.71500000000001 +1.29,0.71000000000001,0.71000000000001,0.71000000000001,0.71000000000001,0.71000000000001,0.71000000000001,0.71000000000001 +1.295,0.70500000000001,0.70500000000001,0.70500000000001,0.70500000000001,0.70500000000001,0.70500000000001,0.70500000000001 +1.3,0.70000000000001,0.70000000000001,0.70000000000001,0.70000000000001,0.70000000000001,0.70000000000001,0.70000000000001 +1.305,0.69500000000001,0.69500000000001,0.69500000000001,0.69500000000001,0.69500000000001,0.69500000000001,0.69500000000001 +1.31,0.69000000000001,0.69000000000001,0.69000000000001,0.69000000000001,0.69000000000001,0.69000000000001,0.69000000000001 +1.315,0.68500000000001,0.68500000000001,0.68500000000001,0.68500000000001,0.68500000000001,0.68500000000001,0.68500000000001 +1.32,0.68000000000001,0.68000000000001,0.68000000000001,0.68000000000001,0.68000000000001,0.68000000000001,0.68000000000001 +1.325,0.67500000000001,0.67500000000001,0.67500000000001,0.67500000000001,0.67500000000001,0.67500000000001,0.67500000000001 +1.33,0.67000000000001,0.67000000000001,0.67000000000001,0.67000000000001,0.67000000000001,0.67000000000001,0.67000000000001 +1.335,0.66500000000001,0.66500000000001,0.66500000000001,0.66500000000001,0.66500000000001,0.66500000000001,0.66500000000001 +1.34,0.66000000000001,0.66000000000001,0.66000000000001,0.66000000000001,0.66000000000001,0.66000000000001,0.66000000000001 +1.345,0.65500000000001,0.65500000000001,0.65500000000001,0.65500000000001,0.65500000000001,0.65500000000001,0.65500000000001 +1.35,0.65000000000001,0.65000000000001,0.65000000000001,0.65000000000001,0.65000000000001,0.65000000000001,0.65000000000001 +1.355,0.64500000000001,0.64500000000001,0.64500000000001,0.64500000000001,0.64500000000001,0.64500000000001,0.64500000000001 +1.36,0.64000000000001,0.64000000000001,0.64000000000001,0.64000000000001,0.64000000000001,0.64000000000001,0.64000000000001 +1.365,0.63500000000001,0.63500000000001,0.63500000000001,0.63500000000001,0.63500000000001,0.63500000000001,0.63500000000001 +1.37,0.63000000000001,0.63000000000001,0.63000000000001,0.63000000000001,0.63000000000001,0.63000000000001,0.63000000000001 +1.375,0.62500000000001,0.62500000000001,0.62500000000001,0.62500000000001,0.62500000000001,0.62500000000001,0.62500000000001 +1.38,0.62000000000001,0.62000000000001,0.62000000000001,0.62000000000001,0.62000000000001,0.62000000000001,0.62000000000001 +1.385,0.61500000000001,0.61500000000001,0.61500000000001,0.61500000000001,0.61500000000001,0.61500000000001,0.61500000000001 +1.39,0.61000000000001,0.61000000000001,0.61000000000001,0.61000000000001,0.61000000000001,0.61000000000001,0.61000000000001 +1.395,0.60500000000001,0.60500000000001,0.60500000000001,0.60500000000001,0.60500000000001,0.60500000000001,0.60500000000001 +1.4,0.60000000000001,0.60000000000001,0.60000000000001,0.60000000000001,0.60000000000001,0.60000000000001,0.60000000000001 +1.405,0.59500000000001,0.59500000000001,0.59500000000001,0.59500000000001,0.59500000000001,0.59500000000001,0.59500000000001 +1.41,0.59000000000001,0.59000000000001,0.59000000000001,0.59000000000001,0.59000000000001,0.59000000000001,0.59000000000001 +1.415,0.58500000000001,0.58500000000001,0.58500000000001,0.58500000000001,0.58500000000001,0.58500000000001,0.58500000000001 +1.42,0.58000000000001,0.58000000000001,0.58000000000001,0.58000000000001,0.58000000000001,0.58000000000001,0.58000000000001 +1.425,0.57500000000001,0.57500000000001,0.57500000000001,0.57500000000001,0.57500000000001,0.57500000000001,0.57500000000001 +1.43,0.57000000000001,0.57000000000001,0.57000000000001,0.57000000000001,0.57000000000001,0.57000000000001,0.57000000000001 +1.435,0.56500000000001,0.56500000000001,0.56500000000001,0.56500000000001,0.56500000000001,0.56500000000001,0.56500000000001 +1.44,0.56000000000001,0.56000000000001,0.56000000000001,0.56000000000001,0.56000000000001,0.56000000000001,0.56000000000001 +1.445,0.55500000000001,0.55500000000001,0.55500000000001,0.55500000000001,0.55500000000001,0.55500000000001,0.55500000000001 +1.45,0.55000000000001,0.55000000000001,0.55000000000001,0.55000000000001,0.55000000000001,0.55000000000001,0.55000000000001 +1.455,0.54500000000001,0.54500000000001,0.54500000000001,0.54500000000001,0.54500000000001,0.54500000000001,0.54500000000001 +1.46,0.54000000000001,0.54000000000001,0.54000000000001,0.54000000000001,0.54000000000001,0.54000000000001,0.54000000000001 +1.465,0.53500000000001,0.53500000000001,0.53500000000001,0.53500000000001,0.53500000000001,0.53500000000001,0.53500000000001 +1.47,0.53000000000001,0.53000000000001,0.53000000000001,0.53000000000001,0.53000000000001,0.53000000000001,0.53000000000001 +1.475,0.52500000000001,0.52500000000001,0.52500000000001,0.52500000000001,0.52500000000001,0.52500000000001,0.52500000000001 +1.48,0.52000000000001,0.52000000000001,0.52000000000001,0.52000000000001,0.52000000000001,0.52000000000001,0.52000000000001 +1.485,0.51500000000001,0.51500000000001,0.51500000000001,0.51500000000001,0.51500000000001,0.51500000000001,0.51500000000001 +1.49,0.51000000000001,0.51000000000001,0.51000000000001,0.51000000000001,0.51000000000001,0.51000000000001,0.51000000000001 +1.495,0.50500000000001,0.50500000000001,0.50500000000001,0.50500000000001,0.50500000000001,0.50500000000001,0.50500000000001 +1.5,0.50000000000001,0.50000000000001,0.50000000000001,0.50000000000001,0.50000000000001,0.50000000000001,0.50000000000001 +1.505,0.49500000000001,0.49500000000001,0.49500000000001,0.49500000000001,0.49500000000001,0.49500000000001,0.49500000000001 +1.51,0.49000000000001,0.49000000000001,0.49000000000001,0.49000000000001,0.49000000000001,0.49000000000001,0.49000000000001 +1.515,0.48500000000001,0.48500000000001,0.48500000000001,0.48500000000001,0.48500000000001,0.48500000000001,0.48500000000001 +1.52,0.48000000000001,0.48000000000001,0.48000000000001,0.48000000000001,0.48000000000001,0.48000000000001,0.48000000000001 +1.525,0.47500000000001,0.47500000000001,0.47500000000001,0.47500000000001,0.47500000000001,0.47500000000001,0.47500000000001 +1.53,0.47000000000001,0.47000000000001,0.47000000000001,0.47000000000001,0.47000000000001,0.47000000000001,0.47000000000001 +1.535,0.46500000000001,0.46500000000001,0.46500000000001,0.46500000000001,0.46500000000001,0.46500000000001,0.46500000000001 +1.54,0.46000000000001,0.46000000000001,0.46000000000001,0.46000000000001,0.46000000000001,0.46000000000001,0.46000000000001 +1.545,0.45500000000001,0.45500000000001,0.45500000000001,0.45500000000001,0.45500000000001,0.45500000000001,0.45500000000001 +1.55,0.45000000000001,0.45000000000001,0.45000000000001,0.45000000000001,0.45000000000001,0.45000000000001,0.45000000000001 +1.555,0.44500000000001,0.44500000000001,0.44500000000001,0.44500000000001,0.44500000000001,0.44500000000001,0.44500000000001 +1.56,0.44000000000001,0.44000000000001,0.44000000000001,0.44000000000001,0.44000000000001,0.44000000000001,0.44000000000001 +1.565,0.43500000000001,0.43500000000001,0.43500000000001,0.43500000000001,0.43500000000001,0.43500000000001,0.43500000000001 +1.57,0.43000000000001,0.43000000000001,0.43000000000001,0.43000000000001,0.43000000000001,0.43000000000001,0.43000000000001 +1.575,0.42500000000001,0.42500000000001,0.42500000000001,0.42500000000001,0.42500000000001,0.42500000000001,0.42500000000001 +1.58,0.42000000000001,0.42000000000001,0.42000000000001,0.42000000000001,0.42000000000001,0.42000000000001,0.42000000000001 +1.585,0.41500000000001,0.41500000000001,0.41500000000001,0.41500000000001,0.41500000000001,0.41500000000001,0.41500000000001 +1.59,0.41000000000001,0.41000000000001,0.41000000000001,0.41000000000001,0.41000000000001,0.41000000000001,0.41000000000001 +1.595,0.40500000000001,0.40500000000001,0.40500000000001,0.40500000000001,0.40500000000001,0.40500000000001,0.40500000000001 +1.6,0.40000000000001,0.40000000000001,0.40000000000001,0.40000000000001,0.40000000000001,0.40000000000001,0.40000000000001 +1.605,0.39500000000001,0.39500000000001,0.39500000000001,0.39500000000001,0.39500000000001,0.39500000000001,0.39500000000001 +1.61,0.39000000000001,0.39000000000001,0.39000000000001,0.39000000000001,0.39000000000001,0.39000000000001,0.39000000000001 +1.615,0.38500000000001,0.38500000000001,0.38500000000001,0.38500000000001,0.38500000000001,0.38500000000001,0.38500000000001 +1.62,0.38000000000001,0.38000000000001,0.38000000000001,0.38000000000001,0.38000000000001,0.38000000000001,0.38000000000001 +1.625,0.37500000000001,0.37500000000001,0.37500000000001,0.37500000000001,0.37500000000001,0.37500000000001,0.37500000000001 +1.63,0.37000000000001,0.37000000000001,0.37000000000001,0.37000000000001,0.37000000000001,0.37000000000001,0.37000000000001 +1.635,0.36500000000001,0.36500000000001,0.36500000000001,0.36500000000001,0.36500000000001,0.36500000000001,0.36500000000001 +1.64,0.36000000000001,0.36000000000001,0.36000000000001,0.36000000000001,0.36000000000001,0.36000000000001,0.36000000000001 +1.645,0.35500000000001,0.35500000000001,0.35500000000001,0.35500000000001,0.35500000000001,0.35500000000001,0.35500000000001 +1.65,0.35000000000001,0.35000000000001,0.35000000000001,0.35000000000001,0.35000000000001,0.35000000000001,0.35000000000001 +1.655,0.34500000000001,0.34500000000001,0.34500000000001,0.34500000000001,0.34500000000001,0.34500000000001,0.34500000000001 +1.66,0.34000000000001,0.34000000000001,0.34000000000001,0.34000000000001,0.34000000000001,0.34000000000001,0.34000000000001 +1.665,0.33500000000001,0.33500000000001,0.33500000000001,0.33500000000001,0.33500000000001,0.33500000000001,0.33500000000001 +1.67,0.33000000000001,0.33000000000001,0.33000000000001,0.33000000000001,0.33000000000001,0.33000000000001,0.33000000000001 +1.675,0.32500000000001,0.32500000000001,0.32500000000001,0.32500000000001,0.32500000000001,0.32500000000001,0.32500000000001 +1.68,0.32000000000001,0.32000000000001,0.32000000000001,0.32000000000001,0.32000000000001,0.32000000000001,0.32000000000001 +1.685,0.31500000000001,0.31500000000001,0.31500000000001,0.31500000000001,0.31500000000001,0.31500000000001,0.31500000000001 +1.69,0.31000000000001,0.31000000000001,0.31000000000001,0.31000000000001,0.31000000000001,0.31000000000001,0.31000000000001 +1.695,0.30500000000001,0.30500000000001,0.30500000000001,0.30500000000001,0.30500000000001,0.30500000000001,0.30500000000001 +1.7,0.30000000000001,0.30000000000001,0.30000000000001,0.30000000000001,0.30000000000001,0.30000000000001,0.30000000000001 +1.705,0.29500000000001,0.29500000000001,0.29500000000001,0.29500000000001,0.29500000000001,0.29500000000001,0.29500000000001 +1.71,0.29000000000001,0.29000000000001,0.29000000000001,0.29000000000001,0.29000000000001,0.29000000000001,0.29000000000001 +1.715,0.28500000000001,0.28500000000001,0.28500000000001,0.28500000000001,0.28500000000001,0.28500000000001,0.28500000000001 +1.72,0.28000000000001,0.28000000000001,0.28000000000001,0.28000000000001,0.28000000000001,0.28000000000001,0.28000000000001 +1.725,0.27500000000001,0.27500000000001,0.27500000000001,0.27500000000001,0.27500000000001,0.27500000000001,0.27500000000001 +1.73,0.27000000000001,0.27000000000001,0.27000000000001,0.27000000000001,0.27000000000001,0.27000000000001,0.27000000000001 +1.735,0.26500000000002,0.26500000000002,0.26500000000002,0.26500000000002,0.26500000000001,0.26500000000001,0.26500000000002 +1.74,0.26000000000002,0.26000000000002,0.26000000000002,0.26000000000002,0.26000000000002,0.26000000000002,0.26000000000002 +1.745,0.25500000000002,0.25500000000002,0.25500000000002,0.25500000000002,0.25500000000002,0.25500000000002,0.25500000000002 +1.75,0.25000000000002,0.25000000000002,0.25000000000002,0.25000000000002,0.25000000000002,0.25000000000002,0.25000000000002 +1.755,0.24500000000002,0.24500000000002,0.24500000000002,0.24500000000002,0.24500000000002,0.24500000000002,0.24500000000002 +1.76,0.24000000000002,0.24000000000002,0.24000000000002,0.24000000000002,0.24000000000002,0.24000000000002,0.24000000000002 +1.765,0.23500000000002,0.23500000000002,0.23500000000002,0.23500000000002,0.23500000000002,0.23500000000002,0.23500000000002 +1.77,0.23000000000002,0.23000000000002,0.23000000000002,0.23000000000002,0.23000000000002,0.23000000000002,0.23000000000002 +1.775,0.22500000000002,0.22500000000002,0.22500000000002,0.22500000000002,0.22500000000002,0.22500000000002,0.22500000000002 +1.78,0.22000000000002,0.22000000000002,0.22000000000002,0.22000000000002,0.22000000000002,0.22000000000002,0.22000000000002 +1.785,0.21500000000002,0.21500000000002,0.21500000000002,0.21500000000002,0.21500000000002,0.21500000000002,0.21500000000002 +1.79,0.21000000000002,0.21000000000002,0.21000000000002,0.21000000000002,0.21000000000002,0.21000000000002,0.21000000000002 +1.795,0.20500000000002,0.20500000000002,0.20500000000002,0.20500000000002,0.20500000000002,0.20500000000002,0.20500000000002 +1.8,0.20000000000002,0.20000000000002,0.20000000000002,0.20000000000002,0.20000000000002,0.20000000000002,0.20000000000002 +1.805,0.19500000000002,0.19500000000002,0.19500000000002,0.19500000000002,0.19500000000002,0.19500000000002,0.19500000000002 +1.81,0.19000000000002,0.19000000000002,0.19000000000002,0.19000000000002,0.19000000000002,0.19000000000002,0.19000000000002 +1.815,0.18500000000002,0.18500000000002,0.18500000000002,0.18500000000002,0.18500000000002,0.18500000000002,0.18500000000002 +1.82,0.18000000000002,0.18000000000002,0.18000000000002,0.18000000000002,0.18000000000002,0.18000000000002,0.18000000000002 +1.825,0.17500000000002,0.17500000000002,0.17500000000002,0.17500000000002,0.17500000000002,0.17500000000002,0.17500000000002 +1.83,0.17000000000002,0.17000000000002,0.17000000000002,0.17000000000002,0.17000000000002,0.17000000000002,0.17000000000002 +1.835,0.16500000000002,0.16500000000002,0.16500000000002,0.16500000000002,0.16500000000002,0.16500000000002,0.16500000000002 +1.84,0.16000000000002,0.16000000000002,0.16000000000002,0.16000000000002,0.16000000000002,0.16000000000002,0.16000000000002 +1.845,0.15500000000002,0.15500000000002,0.15500000000002,0.15500000000002,0.15500000000002,0.15500000000002,0.15500000000002 +1.85,0.15000000000002,0.15000000000002,0.15000000000002,0.15000000000002,0.15000000000002,0.15000000000002,0.15000000000002 +1.855,0.14500000000002,0.14500000000002,0.14500000000002,0.14500000000002,0.14500000000002,0.14500000000002,0.14500000000002 +1.86,0.14000000000002,0.14000000000002,0.14000000000002,0.14000000000002,0.14000000000002,0.14000000000002,0.14000000000002 +1.865,0.13500000000002,0.13500000000002,0.13500000000002,0.13500000000002,0.13500000000002,0.13500000000002,0.13500000000002 +1.87,0.13000000000002,0.13000000000002,0.13000000000002,0.13000000000002,0.13000000000002,0.13000000000002,0.13000000000002 +1.875,0.12500000000002,0.12500000000002,0.12500000000002,0.12500000000002,0.12500000000002,0.12500000000002,0.12500000000002 +1.88,0.12000000000002,0.12000000000002,0.12000000000002,0.12000000000002,0.12000000000002,0.12000000000002,0.12000000000002 +1.885,0.11500000000002,0.11500000000002,0.11500000000002,0.11500000000002,0.11500000000002,0.11500000000002,0.11500000000002 +1.89,0.11000000000002,0.11000000000002,0.11000000000002,0.11000000000002,0.11000000000002,0.11000000000002,0.11000000000002 +1.895,0.10500000000002,0.10500000000002,0.10500000000002,0.10500000000002,0.10500000000002,0.10500000000002,0.10500000000002 +1.9,0.10000000000002,0.10000000000002,0.10000000000002,0.10000000000002,0.10000000000002,0.10000000000002,0.10000000000002 +1.905,0.095000000000019,0.095000000000019,0.095000000000019,0.095000000000019,0.095000000000019,0.095000000000019,0.095000000000019 +1.91,0.090000000000019,0.090000000000019,0.090000000000019,0.090000000000019,0.090000000000019,0.090000000000019,0.090000000000019 +1.915,0.085000000000019,0.085000000000019,0.085000000000019,0.085000000000019,0.085000000000019,0.085000000000019,0.085000000000019 +1.92,0.080000000000019,0.080000000000019,0.080000000000019,0.080000000000019,0.080000000000019,0.080000000000019,0.080000000000019 +1.925,0.075000000000019,0.075000000000019,0.075000000000019,0.075000000000019,0.075000000000019,0.075000000000019,0.075000000000019 +1.93,0.070000000000019,0.070000000000019,0.070000000000019,0.070000000000019,0.070000000000019,0.070000000000019,0.070000000000019 +1.935,0.065000000000019,0.065000000000019,0.065000000000019,0.065000000000019,0.065000000000019,0.065000000000019,0.065000000000019 +1.94,0.060000000000019,0.060000000000019,0.060000000000019,0.060000000000019,0.060000000000019,0.060000000000019,0.060000000000019 +1.945,0.055000000000019,0.055000000000019,0.05500000000002,0.055000000000019,0.055000000000019,0.055000000000019,0.055000000000019 +1.95,0.05000000000002,0.05000000000002,0.05000000000002,0.05000000000002,0.05000000000002,0.05000000000002,0.05000000000002 +1.955,0.04500000000002,0.04500000000002,0.04500000000002,0.04500000000002,0.04500000000002,0.04500000000002,0.04500000000002 +1.96,0.04000000000002,0.04000000000002,0.04000000000002,0.04000000000002,0.04000000000002,0.04000000000002,0.04000000000002 +1.965,0.03500000000002,0.03500000000002,0.03500000000002,0.03500000000002,0.03500000000002,0.03500000000002,0.03500000000002 +1.97,0.03000000000002,0.03000000000002,0.03000000000002,0.03000000000002,0.03000000000002,0.03000000000002,0.03000000000002 +1.975,0.02500000000002,0.02500000000002,0.02500000000002,0.02500000000002,0.02500000000002,0.02500000000002,0.02500000000002 +1.98,0.02000000000002,0.02000000000002,0.02000000000002,0.02000000000002,0.02000000000002,0.02000000000002,0.02000000000002 +1.985,0.01500000000002,0.01500000000002,0.01500000000002,0.01500000000002,0.01500000000002,0.01500000000002,0.01500000000002 +1.99,0.01000000000002,0.01000000000002,0.010000000000021,0.01000000000002,0.01000000000002,0.01000000000002,0.01000000000002 +1.995,0.0050000000000205,0.0050000000000205,0.0050000000000206,0.0050000000000205,0.0050000000000205,0.0050000000000205,0.0050000000000206 +2,2.0650148258028e-14,2.0650148258028e-14,2.0723874005757e-14,2.0650148258028e-14,2.0604178085915e-14,2.0604178085915e-14,2.0605912809391e-14 +2.005,-0.0049999999999795,-0.0049999999999795,-0.0049999999999795,-0.0049999999999795,-0.0049999999999795,-0.0049999999999795,-0.0049999999999795 +2.01,-0.0099999999999794,-0.0099999999999794,-0.0099999999999794,-0.0099999999999794,-0.0099999999999794,-0.0099999999999794,-0.0099999999999794 +2.015,-0.014999999999979,-0.014999999999979,-0.014999999999979,-0.014999999999979,-0.014999999999979,-0.014999999999979,-0.014999999999979 +2.02,-0.019999999999979,-0.019999999999979,-0.019999999999979,-0.019999999999979,-0.019999999999979,-0.019999999999979,-0.019999999999979 +2.025,-0.024999999999979,-0.024999999999979,-0.024999999999979,-0.024999999999979,-0.024999999999979,-0.024999999999979,-0.024999999999979 +2.03,-0.029999999999979,-0.029999999999979,-0.029999999999979,-0.029999999999979,-0.029999999999979,-0.029999999999979,-0.029999999999979 +2.035,-0.034999999999979,-0.034999999999979,-0.034999999999979,-0.034999999999979,-0.034999999999979,-0.034999999999979,-0.034999999999979 +2.04,-0.039999999999979,-0.039999999999979,-0.039999999999979,-0.039999999999979,-0.039999999999979,-0.039999999999979,-0.039999999999979 +2.045,-0.044999999999979,-0.044999999999979,-0.044999999999979,-0.044999999999979,-0.044999999999979,-0.044999999999979,-0.044999999999979 +2.05,-0.049999999999979,-0.049999999999979,-0.049999999999979,-0.049999999999979,-0.049999999999979,-0.049999999999979,-0.049999999999979 +2.055,-0.054999999999978,-0.054999999999978,-0.054999999999978,-0.054999999999978,-0.054999999999978,-0.054999999999978,-0.054999999999978 +2.06,-0.059999999999978,-0.059999999999978,-0.059999999999978,-0.059999999999978,-0.059999999999978,-0.059999999999978,-0.059999999999978 +2.065,-0.064999999999978,-0.064999999999978,-0.064999999999978,-0.064999999999978,-0.064999999999978,-0.064999999999978,-0.064999999999978 +2.07,-0.069999999999978,-0.069999999999978,-0.069999999999978,-0.069999999999978,-0.069999999999978,-0.069999999999978,-0.069999999999978 +2.075,-0.074999999999978,-0.074999999999978,-0.074999999999978,-0.074999999999978,-0.074999999999978,-0.074999999999978,-0.074999999999978 +2.08,-0.079999999999978,-0.079999999999978,-0.079999999999978,-0.079999999999978,-0.079999999999978,-0.079999999999978,-0.079999999999978 +2.085,-0.084999999999978,-0.084999999999978,-0.084999999999978,-0.084999999999978,-0.084999999999978,-0.084999999999978,-0.084999999999978 +2.09,-0.089999999999978,-0.089999999999978,-0.089999999999978,-0.089999999999978,-0.089999999999978,-0.089999999999978,-0.089999999999978 +2.095,-0.094999999999978,-0.094999999999978,-0.094999999999978,-0.094999999999978,-0.094999999999978,-0.094999999999978,-0.094999999999978 +2.1,-0.099999999999977,-0.099999999999977,-0.099999999999977,-0.099999999999977,-0.099999999999977,-0.099999999999977,-0.099999999999977 +2.105,-0.10499999999998,-0.10499999999998,-0.10499999999998,-0.10499999999998,-0.10499999999998,-0.10499999999998,-0.10499999999998 +2.11,-0.10999999999998,-0.10999999999998,-0.10999999999998,-0.10999999999998,-0.10999999999998,-0.10999999999998,-0.10999999999998 +2.115,-0.11499999999998,-0.11499999999998,-0.11499999999998,-0.11499999999998,-0.11499999999998,-0.11499999999998,-0.11499999999998 +2.12,-0.11999999999998,-0.11999999999998,-0.11999999999998,-0.11999999999998,-0.11999999999998,-0.11999999999998,-0.11999999999998 +2.125,-0.12499999999998,-0.12499999999998,-0.12499999999998,-0.12499999999998,-0.12499999999998,-0.12499999999998,-0.12499999999998 +2.13,-0.12999999999998,-0.12999999999998,-0.12999999999998,-0.12999999999998,-0.12999999999998,-0.12999999999998,-0.12999999999998 +2.135,-0.13499999999998,-0.13499999999998,-0.13499999999998,-0.13499999999998,-0.13499999999998,-0.13499999999998,-0.13499999999998 +2.14,-0.13999999999998,-0.13999999999998,-0.13999999999998,-0.13999999999998,-0.13999999999998,-0.13999999999998,-0.13999999999998 +2.145,-0.14499999999998,-0.14499999999998,-0.14499999999998,-0.14499999999998,-0.14499999999998,-0.14499999999998,-0.14499999999998 +2.15,-0.14999999999998,-0.14999999999998,-0.14999999999998,-0.14999999999998,-0.14999999999998,-0.14999999999998,-0.14999999999998 +2.155,-0.15499999999998,-0.15499999999998,-0.15499999999998,-0.15499999999998,-0.15499999999998,-0.15499999999998,-0.15499999999998 +2.16,-0.15999999999998,-0.15999999999998,-0.15999999999998,-0.15999999999998,-0.15999999999998,-0.15999999999998,-0.15999999999998 +2.165,-0.16499999999998,-0.16499999999998,-0.16499999999998,-0.16499999999998,-0.16499999999998,-0.16499999999998,-0.16499999999998 +2.17,-0.16999999999998,-0.16999999999998,-0.16999999999998,-0.16999999999998,-0.16999999999998,-0.16999999999998,-0.16999999999998 +2.175,-0.17499999999998,-0.17499999999998,-0.17499999999998,-0.17499999999998,-0.17499999999998,-0.17499999999998,-0.17499999999998 +2.18,-0.17999999999998,-0.17999999999998,-0.17999999999998,-0.17999999999998,-0.17999999999998,-0.17999999999998,-0.17999999999998 +2.185,-0.18499999999998,-0.18499999999998,-0.18499999999998,-0.18499999999998,-0.18499999999998,-0.18499999999998,-0.18499999999998 +2.19,-0.18999999999998,-0.18999999999998,-0.18999999999998,-0.18999999999998,-0.18999999999998,-0.18999999999998,-0.18999999999998 +2.195,-0.19499999999998,-0.19499999999998,-0.19499999999998,-0.19499999999998,-0.19499999999998,-0.19499999999998,-0.19499999999998 +2.2,-0.19999999999998,-0.19999999999998,-0.19999999999998,-0.19999999999998,-0.19999999999998,-0.19999999999998,-0.19999999999998 +2.205,-0.20499999999998,-0.20499999999998,-0.20499999999998,-0.20499999999998,-0.20499999999998,-0.20499999999998,-0.20499999999998 +2.21,-0.20999999999998,-0.20999999999998,-0.20999999999998,-0.20999999999998,-0.20999999999998,-0.20999999999998,-0.20999999999998 +2.215,-0.21499999999997,-0.21499999999997,-0.21499999999997,-0.21499999999997,-0.21499999999997,-0.21499999999997,-0.21499999999997 +2.22,-0.21999999999997,-0.21999999999997,-0.21999999999997,-0.21999999999997,-0.21999999999997,-0.21999999999997,-0.21999999999997 +2.225,-0.22499999999997,-0.22499999999997,-0.22499999999997,-0.22499999999997,-0.22499999999997,-0.22499999999997,-0.22499999999997 +2.23,-0.22999999999997,-0.22999999999997,-0.22999999999997,-0.22999999999997,-0.22999999999997,-0.22999999999997,-0.22999999999997 +2.235,-0.23499999999997,-0.23499999999997,-0.23499999999997,-0.23499999999997,-0.23499999999997,-0.23499999999997,-0.23499999999997 +2.24,-0.23999999999997,-0.23999999999997,-0.23999999999997,-0.23999999999997,-0.23999999999997,-0.23999999999997,-0.23999999999997 +2.245,-0.24499999999997,-0.24499999999997,-0.24499999999997,-0.24499999999997,-0.24499999999997,-0.24499999999997,-0.24499999999997 +2.25,-0.24999999999997,-0.24999999999997,-0.24999999999997,-0.24999999999997,-0.24999999999997,-0.24999999999997,-0.24999999999997 +2.255,-0.25499999999997,-0.25499999999997,-0.25499999999997,-0.25499999999997,-0.25499999999997,-0.25499999999997,-0.25499999999997 +2.26,-0.25999999999997,-0.25999999999997,-0.25999999999997,-0.25999999999997,-0.25999999999997,-0.25999999999997,-0.25999999999997 +2.265,-0.26499999999997,-0.26499999999997,-0.26499999999997,-0.26499999999997,-0.26499999999997,-0.26499999999997,-0.26499999999997 +2.27,-0.26999999999997,-0.26999999999997,-0.26999999999997,-0.26999999999997,-0.26999999999997,-0.26999999999997,-0.26999999999997 +2.275,-0.27499999999997,-0.27499999999997,-0.27499999999997,-0.27499999999997,-0.27499999999997,-0.27499999999997,-0.27499999999997 +2.28,-0.27999999999997,-0.27999999999997,-0.27999999999997,-0.27999999999997,-0.27999999999997,-0.27999999999997,-0.27999999999997 +2.285,-0.28499999999997,-0.28499999999997,-0.28499999999997,-0.28499999999997,-0.28499999999997,-0.28499999999997,-0.28499999999997 +2.29,-0.28999999999997,-0.28999999999997,-0.28999999999997,-0.28999999999997,-0.28999999999997,-0.28999999999997,-0.28999999999997 +2.295,-0.29499999999997,-0.29499999999997,-0.29499999999997,-0.29499999999997,-0.29499999999997,-0.29499999999997,-0.29499999999997 +2.3,-0.29999999999997,-0.29999999999997,-0.29999999999997,-0.29999999999997,-0.29999999999997,-0.29999999999997,-0.29999999999997 +2.305,-0.30499999999997,-0.30499999999997,-0.30499999999997,-0.30499999999997,-0.30499999999997,-0.30499999999997,-0.30499999999997 +2.31,-0.30999999999997,-0.30999999999997,-0.30999999999997,-0.30999999999997,-0.30999999999997,-0.30999999999997,-0.30999999999997 +2.315,-0.31499999999997,-0.31499999999997,-0.31499999999997,-0.31499999999997,-0.31499999999997,-0.31499999999997,-0.31499999999997 +2.32,-0.31999999999997,-0.31999999999997,-0.31999999999997,-0.31999999999997,-0.31999999999997,-0.31999999999997,-0.31999999999997 +2.325,-0.32499999999997,-0.32499999999997,-0.32499999999997,-0.32499999999997,-0.32499999999997,-0.32499999999997,-0.32499999999997 +2.33,-0.32999999999997,-0.32999999999997,-0.32999999999997,-0.32999999999997,-0.32999999999997,-0.32999999999997,-0.32999999999997 +2.335,-0.33499999999997,-0.33499999999997,-0.33499999999997,-0.33499999999997,-0.33499999999997,-0.33499999999997,-0.33499999999997 +2.34,-0.33999999999997,-0.33999999999997,-0.33999999999997,-0.33999999999997,-0.33999999999997,-0.33999999999997,-0.33999999999997 +2.345,-0.34499999999997,-0.34499999999997,-0.34499999999997,-0.34499999999997,-0.34499999999997,-0.34499999999997,-0.34499999999997 +2.35,-0.34999999999997,-0.34999999999997,-0.34999999999997,-0.34999999999997,-0.34999999999997,-0.34999999999997,-0.34999999999997 +2.355,-0.35499999999997,-0.35499999999997,-0.35499999999997,-0.35499999999997,-0.35499999999997,-0.35499999999997,-0.35499999999997 +2.36,-0.35999999999997,-0.35999999999997,-0.35999999999997,-0.35999999999997,-0.35999999999997,-0.35999999999997,-0.35999999999997 +2.365,-0.36499999999997,-0.36499999999997,-0.36499999999997,-0.36499999999997,-0.36499999999997,-0.36499999999997,-0.36499999999997 +2.37,-0.36999999999997,-0.36999999999997,-0.36999999999997,-0.36999999999997,-0.36999999999997,-0.36999999999997,-0.36999999999997 +2.375,-0.37499999999997,-0.37499999999997,-0.37499999999997,-0.37499999999997,-0.37499999999997,-0.37499999999997,-0.37499999999997 +2.38,-0.37999999999997,-0.37999999999997,-0.37999999999997,-0.37999999999997,-0.37999999999997,-0.37999999999997,-0.37999999999997 +2.385,-0.38499999999997,-0.38499999999997,-0.38499999999997,-0.38499999999997,-0.38499999999997,-0.38499999999997,-0.38499999999997 +2.39,-0.38999999999997,-0.38999999999997,-0.38999999999997,-0.38999999999997,-0.38999999999997,-0.38999999999997,-0.38999999999997 +2.395,-0.39499999999997,-0.39499999999997,-0.39499999999997,-0.39499999999997,-0.39499999999997,-0.39499999999997,-0.39499999999997 +2.4,-0.39999999999997,-0.39999999999997,-0.39999999999997,-0.39999999999997,-0.39999999999997,-0.39999999999997,-0.39999999999997 +2.405,-0.40499999999997,-0.40499999999997,-0.40499999999997,-0.40499999999997,-0.40499999999997,-0.40499999999997,-0.40499999999997 +2.41,-0.40999999999997,-0.40999999999997,-0.40999999999997,-0.40999999999997,-0.40999999999997,-0.40999999999997,-0.40999999999997 +2.415,-0.41499999999997,-0.41499999999997,-0.41499999999997,-0.41499999999997,-0.41499999999997,-0.41499999999997,-0.41499999999997 +2.42,-0.41999999999997,-0.41999999999997,-0.41999999999997,-0.41999999999997,-0.41999999999997,-0.41999999999997,-0.41999999999997 +2.425,-0.42499999999997,-0.42499999999997,-0.42499999999997,-0.42499999999997,-0.42499999999997,-0.42499999999997,-0.42499999999997 +2.43,-0.42999999999997,-0.42999999999997,-0.42999999999997,-0.42999999999997,-0.42999999999997,-0.42999999999997,-0.42999999999997 +2.435,-0.43499999999997,-0.43499999999997,-0.43499999999997,-0.43499999999997,-0.43499999999997,-0.43499999999997,-0.43499999999997 +2.44,-0.43999999999997,-0.43999999999997,-0.43999999999997,-0.43999999999997,-0.43999999999997,-0.43999999999997,-0.43999999999997 +2.445,-0.44499999999997,-0.44499999999997,-0.44499999999997,-0.44499999999997,-0.44499999999997,-0.44499999999997,-0.44499999999997 +2.45,-0.44999999999997,-0.44999999999997,-0.44999999999997,-0.44999999999997,-0.44999999999997,-0.44999999999997,-0.44999999999997 +2.455,-0.45499999999997,-0.45499999999997,-0.45499999999997,-0.45499999999997,-0.45499999999997,-0.45499999999997,-0.45499999999997 +2.46,-0.45999999999997,-0.45999999999997,-0.45999999999997,-0.45999999999997,-0.45999999999997,-0.45999999999997,-0.45999999999997 +2.465,-0.46499999999997,-0.46499999999997,-0.46499999999997,-0.46499999999997,-0.46499999999997,-0.46499999999997,-0.46499999999997 +2.47,-0.46999999999997,-0.46999999999997,-0.46999999999997,-0.46999999999997,-0.46999999999997,-0.46999999999997,-0.46999999999997 +2.475,-0.47499999999997,-0.47499999999997,-0.47499999999997,-0.47499999999997,-0.47499999999997,-0.47499999999997,-0.47499999999997 +2.48,-0.47999999999997,-0.47999999999997,-0.47999999999997,-0.47999999999997,-0.47999999999997,-0.47999999999997,-0.47999999999997 +2.485,-0.48499999999997,-0.48499999999997,-0.48499999999997,-0.48499999999997,-0.48499999999997,-0.48499999999997,-0.48499999999997 +2.49,-0.48999999999997,-0.48999999999997,-0.48999999999997,-0.48999999999997,-0.48999999999997,-0.48999999999997,-0.48999999999997 +2.495,-0.49499999999997,-0.49499999999997,-0.49499999999997,-0.49499999999997,-0.49499999999997,-0.49499999999997,-0.49499999999997 +2.5,-0.49999999999997,-0.49999999999997,-0.49999999999997,-0.49999999999997,-0.49999999999997,-0.49999999999997,-0.49999999999997 +2.505,-0.50499999999997,-0.50499999999997,-0.50499999999997,-0.50499999999997,-0.50499999999997,-0.50499999999997,-0.50499999999997 +2.51,-0.50999999999997,-0.50999999999997,-0.50999999999997,-0.50999999999997,-0.50999999999997,-0.50999999999997,-0.50999999999997 +2.515,-0.51499999999997,-0.51499999999997,-0.51499999999997,-0.51499999999997,-0.51499999999997,-0.51499999999997,-0.51499999999997 +2.52,-0.51999999999997,-0.51999999999997,-0.51999999999997,-0.51999999999997,-0.51999999999997,-0.51999999999997,-0.51999999999997 +2.525,-0.52499999999997,-0.52499999999997,-0.52499999999997,-0.52499999999997,-0.52499999999997,-0.52499999999997,-0.52499999999997 +2.53,-0.52999999999997,-0.52999999999997,-0.52999999999997,-0.52999999999997,-0.52999999999997,-0.52999999999997,-0.52999999999997 +2.535,-0.53499999999997,-0.53499999999997,-0.53499999999997,-0.53499999999997,-0.53499999999997,-0.53499999999997,-0.53499999999997 +2.54,-0.53999999999997,-0.53999999999997,-0.53999999999997,-0.53999999999997,-0.53999999999997,-0.53999999999997,-0.53999999999997 +2.545,-0.54499999999997,-0.54499999999997,-0.54499999999997,-0.54499999999997,-0.54499999999997,-0.54499999999997,-0.54499999999997 +2.55,-0.54999999999997,-0.54999999999997,-0.54999999999997,-0.54999999999997,-0.54999999999997,-0.54999999999997,-0.54999999999997 +2.555,-0.55499999999997,-0.55499999999997,-0.55499999999997,-0.55499999999997,-0.55499999999997,-0.55499999999997,-0.55499999999997 +2.56,-0.55999999999997,-0.55999999999997,-0.55999999999997,-0.55999999999997,-0.55999999999997,-0.55999999999997,-0.55999999999997 +2.565,-0.56499999999997,-0.56499999999997,-0.56499999999997,-0.56499999999997,-0.56499999999997,-0.56499999999997,-0.56499999999997 +2.57,-0.56999999999997,-0.56999999999997,-0.56999999999997,-0.56999999999997,-0.56999999999997,-0.56999999999997,-0.56999999999997 +2.575,-0.57499999999997,-0.57499999999997,-0.57499999999997,-0.57499999999997,-0.57499999999997,-0.57499999999997,-0.57499999999997 +2.58,-0.57999999999997,-0.57999999999997,-0.57999999999997,-0.57999999999997,-0.57999999999997,-0.57999999999997,-0.57999999999997 +2.585,-0.58499999999997,-0.58499999999997,-0.58499999999997,-0.58499999999997,-0.58499999999997,-0.58499999999997,-0.58499999999997 +2.59,-0.58999999999997,-0.58999999999997,-0.58999999999997,-0.58999999999997,-0.58999999999997,-0.58999999999997,-0.58999999999997 +2.595,-0.59499999999997,-0.59499999999997,-0.59499999999997,-0.59499999999997,-0.59499999999997,-0.59499999999997,-0.59499999999997 +2.6,-0.59999999999997,-0.59999999999997,-0.59999999999997,-0.59999999999997,-0.59999999999997,-0.59999999999997,-0.59999999999997 +2.605,-0.60499999999997,-0.60499999999997,-0.60499999999997,-0.60499999999997,-0.60499999999997,-0.60499999999997,-0.60499999999997 +2.61,-0.60999999999997,-0.60999999999997,-0.60999999999997,-0.60999999999997,-0.60999999999997,-0.60999999999997,-0.60999999999997 +2.615,-0.61499999999997,-0.61499999999997,-0.61499999999997,-0.61499999999997,-0.61499999999997,-0.61499999999997,-0.61499999999997 +2.62,-0.61999999999997,-0.61999999999997,-0.61999999999997,-0.61999999999997,-0.61999999999997,-0.61999999999997,-0.61999999999997 +2.625,-0.62499999999997,-0.62499999999997,-0.62499999999997,-0.62499999999997,-0.62499999999997,-0.62499999999997,-0.62499999999997 +2.63,-0.62999999999997,-0.62999999999997,-0.62999999999997,-0.62999999999997,-0.62999999999997,-0.62999999999997,-0.62999999999997 +2.635,-0.63499999999997,-0.63499999999997,-0.63499999999997,-0.63499999999997,-0.63499999999997,-0.63499999999997,-0.63499999999997 +2.64,-0.63999999999997,-0.63999999999997,-0.63999999999997,-0.63999999999997,-0.63999999999997,-0.63999999999997,-0.63999999999997 +2.645,-0.64499999999997,-0.64499999999997,-0.64499999999997,-0.64499999999997,-0.64499999999997,-0.64499999999997,-0.64499999999997 +2.65,-0.64999999999997,-0.64999999999997,-0.64999999999997,-0.64999999999997,-0.64999999999997,-0.64999999999997,-0.64999999999997 +2.655,-0.65499999999997,-0.65499999999997,-0.65499999999997,-0.65499999999997,-0.65499999999997,-0.65499999999997,-0.65499999999997 +2.66,-0.65999999999997,-0.65999999999997,-0.65999999999997,-0.65999999999997,-0.65999999999997,-0.65999999999997,-0.65999999999997 +2.665,-0.66499999999997,-0.66499999999997,-0.66499999999997,-0.66499999999997,-0.66499999999997,-0.66499999999997,-0.66499999999997 +2.67,-0.66999999999997,-0.66999999999997,-0.66999999999997,-0.66999999999997,-0.66999999999997,-0.66999999999997,-0.66999999999997 +2.675,-0.67499999999997,-0.67499999999997,-0.67499999999997,-0.67499999999997,-0.67499999999997,-0.67499999999997,-0.67499999999997 +2.68,-0.67999999999997,-0.67999999999997,-0.67999999999997,-0.67999999999997,-0.67999999999997,-0.67999999999997,-0.67999999999997 +2.685,-0.68499999999996,-0.68499999999996,-0.68499999999996,-0.68499999999996,-0.68499999999996,-0.68499999999996,-0.68499999999996 +2.69,-0.68999999999996,-0.68999999999996,-0.68999999999996,-0.68999999999996,-0.68999999999996,-0.68999999999996,-0.68999999999996 +2.695,-0.69499999999996,-0.69499999999996,-0.69499999999996,-0.69499999999996,-0.69499999999996,-0.69499999999996,-0.69499999999996 +2.7,-0.69999999999996,-0.69999999999996,-0.69999999999996,-0.69999999999996,-0.69999999999996,-0.69999999999996,-0.69999999999996 +2.705,-0.70499999999996,-0.70499999999996,-0.70499999999996,-0.70499999999996,-0.70499999999996,-0.70499999999996,-0.70499999999996 +2.71,-0.70999999999996,-0.70999999999996,-0.70999999999996,-0.70999999999996,-0.70999999999996,-0.70999999999996,-0.70999999999996 +2.715,-0.71499999999996,-0.71499999999996,-0.71499999999996,-0.71499999999996,-0.71499999999996,-0.71499999999996,-0.71499999999996 +2.72,-0.71999999999996,-0.71999999999996,-0.71999999999996,-0.71999999999996,-0.71999999999996,-0.71999999999996,-0.71999999999996 +2.725,-0.72499999999996,-0.72499999999996,-0.72499999999996,-0.72499999999996,-0.72499999999996,-0.72499999999996,-0.72499999999996 +2.73,-0.72999999999996,-0.72999999999996,-0.72999999999996,-0.72999999999996,-0.72999999999996,-0.72999999999996,-0.72999999999996 +2.735,-0.73499999999996,-0.73499999999996,-0.73499999999996,-0.73499999999996,-0.73499999999996,-0.73499999999996,-0.73499999999996 +2.74,-0.73999999999996,-0.73999999999996,-0.73999999999996,-0.73999999999996,-0.73999999999996,-0.73999999999996,-0.73999999999996 +2.745,-0.74499999999996,-0.74499999999996,-0.74499999999996,-0.74499999999996,-0.74499999999996,-0.74499999999996,-0.74499999999996 +2.75,-0.74999999999996,-0.74999999999996,-0.74999999999996,-0.74999999999996,-0.74999999999996,-0.74999999999996,-0.74999999999996 +2.755,-0.75499999999996,-0.75499999999996,-0.75499999999996,-0.75499999999996,-0.75499999999996,-0.75499999999996,-0.75499999999996 +2.76,-0.75999999999996,-0.75999999999996,-0.75999999999996,-0.75999999999996,-0.75999999999996,-0.75999999999996,-0.75999999999996 +2.765,-0.76499999999996,-0.76499999999996,-0.76499999999996,-0.76499999999996,-0.76499999999996,-0.76499999999996,-0.76499999999996 +2.77,-0.76999999999996,-0.76999999999996,-0.76999999999996,-0.76999999999996,-0.76999999999996,-0.76999999999996,-0.76999999999996 +2.775,-0.77499999999996,-0.77499999999996,-0.77499999999996,-0.77499999999996,-0.77499999999996,-0.77499999999996,-0.77499999999996 +2.78,-0.77999999999996,-0.77999999999996,-0.77999999999996,-0.77999999999996,-0.77999999999996,-0.77999999999996,-0.77999999999996 +2.785,-0.78499999999996,-0.78499999999996,-0.78499999999996,-0.78499999999996,-0.78499999999996,-0.78499999999996,-0.78499999999996 +2.79,-0.78999999999996,-0.78999999999996,-0.78999999999996,-0.78999999999996,-0.78999999999996,-0.78999999999996,-0.78999999999996 +2.795,-0.79499999999996,-0.79499999999996,-0.79499999999996,-0.79499999999996,-0.79499999999996,-0.79499999999996,-0.79499999999996 +2.8,-0.79999999999996,-0.79999999999996,-0.79999999999996,-0.79999999999996,-0.79999999999996,-0.79999999999996,-0.79999999999996 +2.805,-0.80499999999996,-0.80499999999996,-0.80499999999996,-0.80499999999996,-0.80499999999996,-0.80499999999996,-0.80499999999996 +2.81,-0.80999999999996,-0.80999999999996,-0.80999999999996,-0.80999999999996,-0.80999999999996,-0.80999999999996,-0.80999999999996 +2.815,-0.81499999999996,-0.81499999999996,-0.81499999999996,-0.81499999999996,-0.81499999999996,-0.81499999999996,-0.81499999999996 +2.82,-0.81999999999996,-0.81999999999996,-0.81999999999996,-0.81999999999996,-0.81999999999996,-0.81999999999996,-0.81999999999996 +2.825,-0.82499999999996,-0.82499999999996,-0.82499999999996,-0.82499999999996,-0.82499999999996,-0.82499999999996,-0.82499999999996 +2.83,-0.82999999999996,-0.82999999999996,-0.82999999999996,-0.82999999999996,-0.82999999999996,-0.82999999999996,-0.82999999999996 +2.835,-0.83499999999996,-0.83499999999996,-0.83499999999996,-0.83499999999996,-0.83499999999996,-0.83499999999996,-0.83499999999996 +2.84,-0.83999999999996,-0.83999999999996,-0.83999999999996,-0.83999999999996,-0.83999999999996,-0.83999999999996,-0.83999999999996 +2.845,-0.84499999999996,-0.84499999999996,-0.84499999999996,-0.84499999999996,-0.84499999999996,-0.84499999999996,-0.84499999999996 +2.85,-0.84999999999996,-0.84999999999996,-0.84999999999996,-0.84999999999996,-0.84999999999996,-0.84999999999996,-0.84999999999996 +2.855,-0.85499999999996,-0.85499999999996,-0.85499999999996,-0.85499999999996,-0.85499999999996,-0.85499999999996,-0.85499999999996 +2.86,-0.85999999999996,-0.85999999999996,-0.85999999999996,-0.85999999999996,-0.85999999999996,-0.85999999999996,-0.85999999999996 +2.865,-0.86499999999996,-0.86499999999996,-0.86499999999996,-0.86499999999996,-0.86499999999996,-0.86499999999996,-0.86499999999996 +2.87,-0.86999999999996,-0.86999999999996,-0.86999999999996,-0.86999999999996,-0.86999999999996,-0.86999999999996,-0.86999999999996 +2.875,-0.87499999999996,-0.87499999999996,-0.87499999999996,-0.87499999999996,-0.87499999999996,-0.87499999999996,-0.87499999999996 +2.88,-0.87999999999996,-0.87999999999996,-0.87999999999996,-0.87999999999996,-0.87999999999996,-0.87999999999996,-0.87999999999996 +2.885,-0.88499999999996,-0.88499999999996,-0.88499999999996,-0.88499999999996,-0.88499999999996,-0.88499999999996,-0.88499999999996 +2.89,-0.88999999999996,-0.88999999999996,-0.88999999999996,-0.88999999999996,-0.88999999999996,-0.88999999999996,-0.88999999999996 +2.895,-0.89499999999996,-0.89499999999996,-0.89499999999996,-0.89499999999996,-0.89499999999996,-0.89499999999996,-0.89499999999996 +2.9,-0.89999999999996,-0.89999999999996,-0.89999999999996,-0.89999999999996,-0.89999999999996,-0.89999999999996,-0.89999999999996 +2.905,-0.90499999999996,-0.90499999999996,-0.90499999999996,-0.90499999999996,-0.90499999999996,-0.90499999999996,-0.90499999999996 +2.91,-0.90999999999996,-0.90999999999996,-0.90999999999996,-0.90999999999996,-0.90999999999996,-0.90999999999996,-0.90999999999996 +2.915,-0.91499999999996,-0.91499999999996,-0.91499999999996,-0.91499999999996,-0.91499999999996,-0.91499999999996,-0.91499999999996 +2.92,-0.91999999999996,-0.91999999999996,-0.91999999999996,-0.91999999999996,-0.91999999999996,-0.91999999999996,-0.91999999999996 +2.925,-0.92499999999996,-0.92499999999996,-0.92499999999996,-0.92499999999996,-0.92499999999996,-0.92499999999996,-0.92499999999996 +2.93,-0.92999999999996,-0.92999999999996,-0.92999999999996,-0.92999999999996,-0.92999999999996,-0.92999999999996,-0.92999999999996 +2.935,-0.93499999999996,-0.93499999999996,-0.93499999999996,-0.93499999999996,-0.93499999999996,-0.93499999999996,-0.93499999999996 +2.94,-0.93999999999996,-0.93999999999996,-0.93999999999996,-0.93999999999996,-0.93999999999996,-0.93999999999996,-0.93999999999996 +2.945,-0.94499999999996,-0.94499999999996,-0.94499999999996,-0.94499999999996,-0.94499999999996,-0.94499999999996,-0.94499999999996 +2.95,-0.94999999999996,-0.94999999999996,-0.94999999999996,-0.94999999999996,-0.94999999999996,-0.94999999999996,-0.94999999999996 +2.955,-0.95499999999996,-0.95499999999996,-0.95499999999996,-0.95499999999996,-0.95499999999996,-0.95499999999996,-0.95499999999996 +2.96,-0.95999999999996,-0.95999999999996,-0.95999999999996,-0.95999999999996,-0.95999999999996,-0.95999999999996,-0.95999999999996 +2.965,-0.96499999999996,-0.96499999999996,-0.96499999999996,-0.96499999999996,-0.96499999999996,-0.96499999999996,-0.96499999999996 +2.97,-0.96999999999996,-0.96999999999996,-0.96999999999996,-0.96999999999996,-0.96999999999996,-0.96999999999996,-0.96999999999996 +2.975,-0.97499999999996,-0.97499999999996,-0.97499999999996,-0.97499999999996,-0.97499999999996,-0.97499999999996,-0.97499999999996 +2.98,-0.97999999999996,-0.97999999999996,-0.97999999999996,-0.97999999999996,-0.97999999999996,-0.97999999999996,-0.97999999999996 +2.985,-0.98499999999996,-0.98499999999996,-0.98499999999996,-0.98499999999996,-0.98499999999996,-0.98499999999996,-0.98499999999996 +2.99,-0.98999999999996,-0.98999999999996,-0.98999999999996,-0.98999999999996,-0.98999999999996,-0.98999999999996,-0.98999999999996 +2.995,-0.99499999999996,-0.99499999999996,-0.99499999999996,-0.99499999999996,-0.99499999999996,-0.99499999999996,-0.99499999999996 +3,-0.99999999999996,-0.99999999999996,-0.99999999999996,-0.99999999999996,-0.99999999999996,-0.99999999999996,-0.99999999999996 +3.005,-0.99500000000004,-0.99500000000004,-0.99500000000004,-0.99500000000004,-0.99500000000004,-0.99500000000004,-0.99500000000004 +3.01,-0.99000000000004,-0.99000000000004,-0.99000000000004,-0.99000000000004,-0.99000000000004,-0.99000000000004,-0.99000000000004 +3.015,-0.98500000000004,-0.98500000000004,-0.98500000000004,-0.98500000000004,-0.98500000000004,-0.98500000000004,-0.98500000000004 +3.02,-0.98000000000004,-0.98000000000004,-0.98000000000004,-0.98000000000004,-0.98000000000004,-0.98000000000004,-0.98000000000004 +3.025,-0.97500000000004,-0.97500000000004,-0.97500000000004,-0.97500000000004,-0.97500000000004,-0.97500000000004,-0.97500000000004 +3.03,-0.97000000000004,-0.97000000000004,-0.97000000000004,-0.97000000000004,-0.97000000000004,-0.97000000000004,-0.97000000000004 +3.035,-0.96500000000004,-0.96500000000004,-0.96500000000004,-0.96500000000004,-0.96500000000004,-0.96500000000004,-0.96500000000004 +3.04,-0.96000000000004,-0.96000000000004,-0.96000000000004,-0.96000000000004,-0.96000000000004,-0.96000000000004,-0.96000000000004 +3.045,-0.95500000000004,-0.95500000000004,-0.95500000000004,-0.95500000000004,-0.95500000000004,-0.95500000000004,-0.95500000000004 +3.05,-0.95000000000004,-0.95000000000004,-0.95000000000004,-0.95000000000004,-0.95000000000004,-0.95000000000004,-0.95000000000004 +3.055,-0.94500000000004,-0.94500000000004,-0.94500000000004,-0.94500000000004,-0.94500000000004,-0.94500000000004,-0.94500000000004 +3.06,-0.94000000000004,-0.94000000000004,-0.94000000000004,-0.94000000000004,-0.94000000000004,-0.94000000000004,-0.94000000000004 +3.065,-0.93500000000004,-0.93500000000004,-0.93500000000004,-0.93500000000004,-0.93500000000004,-0.93500000000004,-0.93500000000004 +3.07,-0.93000000000004,-0.93000000000004,-0.93000000000004,-0.93000000000004,-0.93000000000004,-0.93000000000004,-0.93000000000004 +3.075,-0.92500000000004,-0.92500000000004,-0.92500000000004,-0.92500000000004,-0.92500000000004,-0.92500000000004,-0.92500000000004 +3.08,-0.92000000000004,-0.92000000000004,-0.92000000000004,-0.92000000000004,-0.92000000000004,-0.92000000000004,-0.92000000000004 +3.085,-0.91500000000004,-0.91500000000004,-0.91500000000004,-0.91500000000004,-0.91500000000004,-0.91500000000004,-0.91500000000004 +3.09,-0.91000000000004,-0.91000000000004,-0.91000000000004,-0.91000000000004,-0.91000000000004,-0.91000000000004,-0.91000000000004 +3.095,-0.90500000000004,-0.90500000000004,-0.90500000000004,-0.90500000000004,-0.90500000000004,-0.90500000000004,-0.90500000000004 +3.1,-0.90000000000004,-0.90000000000004,-0.90000000000004,-0.90000000000004,-0.90000000000004,-0.90000000000004,-0.90000000000004 +3.105,-0.89500000000004,-0.89500000000004,-0.89500000000004,-0.89500000000004,-0.89500000000004,-0.89500000000004,-0.89500000000004 +3.11,-0.89000000000004,-0.89000000000004,-0.89000000000004,-0.89000000000004,-0.89000000000004,-0.89000000000004,-0.89000000000004 +3.115,-0.88500000000004,-0.88500000000004,-0.88500000000004,-0.88500000000004,-0.88500000000004,-0.88500000000004,-0.88500000000004 +3.12,-0.88000000000004,-0.88000000000004,-0.88000000000004,-0.88000000000004,-0.88000000000004,-0.88000000000004,-0.88000000000004 +3.125,-0.87500000000004,-0.87500000000004,-0.87500000000004,-0.87500000000004,-0.87500000000004,-0.87500000000004,-0.87500000000004 +3.13,-0.87000000000004,-0.87000000000004,-0.87000000000004,-0.87000000000004,-0.87000000000004,-0.87000000000004,-0.87000000000004 +3.135,-0.86500000000004,-0.86500000000004,-0.86500000000004,-0.86500000000004,-0.86500000000004,-0.86500000000004,-0.86500000000004 +3.14,-0.86000000000004,-0.86000000000004,-0.86000000000004,-0.86000000000004,-0.86000000000004,-0.86000000000004,-0.86000000000004 +3.145,-0.85500000000004,-0.85500000000004,-0.85500000000004,-0.85500000000004,-0.85500000000004,-0.85500000000004,-0.85500000000004 +3.15,-0.85000000000004,-0.85000000000004,-0.85000000000004,-0.85000000000004,-0.85000000000004,-0.85000000000004,-0.85000000000005 +3.155,-0.84500000000005,-0.84500000000005,-0.84500000000005,-0.84500000000005,-0.84500000000005,-0.84500000000005,-0.84500000000005 +3.16,-0.84000000000005,-0.84000000000005,-0.84000000000005,-0.84000000000005,-0.84000000000005,-0.84000000000005,-0.84000000000005 +3.165,-0.83500000000005,-0.83500000000005,-0.83500000000005,-0.83500000000005,-0.83500000000005,-0.83500000000005,-0.83500000000005 +3.17,-0.83000000000005,-0.83000000000005,-0.83000000000005,-0.83000000000005,-0.83000000000005,-0.83000000000005,-0.83000000000005 +3.175,-0.82500000000005,-0.82500000000005,-0.82500000000005,-0.82500000000005,-0.82500000000005,-0.82500000000005,-0.82500000000005 +3.18,-0.82000000000005,-0.82000000000005,-0.82000000000005,-0.82000000000005,-0.82000000000005,-0.82000000000005,-0.82000000000005 +3.185,-0.81500000000005,-0.81500000000005,-0.81500000000005,-0.81500000000005,-0.81500000000005,-0.81500000000005,-0.81500000000005 +3.19,-0.81000000000005,-0.81000000000005,-0.81000000000005,-0.81000000000005,-0.81000000000005,-0.81000000000005,-0.81000000000005 +3.195,-0.80500000000005,-0.80500000000005,-0.80500000000005,-0.80500000000005,-0.80500000000005,-0.80500000000005,-0.80500000000005 +3.2,-0.80000000000005,-0.80000000000005,-0.80000000000005,-0.80000000000005,-0.80000000000005,-0.80000000000005,-0.80000000000005 +3.205,-0.79500000000005,-0.79500000000005,-0.79500000000005,-0.79500000000005,-0.79500000000005,-0.79500000000005,-0.79500000000005 +3.21,-0.79000000000005,-0.79000000000005,-0.79000000000005,-0.79000000000005,-0.79000000000005,-0.79000000000005,-0.79000000000005 +3.215,-0.78500000000005,-0.78500000000005,-0.78500000000005,-0.78500000000005,-0.78500000000005,-0.78500000000005,-0.78500000000005 +3.22,-0.78000000000005,-0.78000000000005,-0.78000000000005,-0.78000000000005,-0.78000000000005,-0.78000000000005,-0.78000000000005 +3.225,-0.77500000000005,-0.77500000000005,-0.77500000000005,-0.77500000000005,-0.77500000000005,-0.77500000000005,-0.77500000000005 +3.23,-0.77000000000005,-0.77000000000005,-0.77000000000005,-0.77000000000005,-0.77000000000005,-0.77000000000005,-0.77000000000005 +3.235,-0.76500000000005,-0.76500000000005,-0.76500000000005,-0.76500000000005,-0.76500000000005,-0.76500000000005,-0.76500000000005 +3.24,-0.76000000000005,-0.76000000000005,-0.76000000000005,-0.76000000000005,-0.76000000000005,-0.76000000000005,-0.76000000000005 +3.245,-0.75500000000005,-0.75500000000005,-0.75500000000005,-0.75500000000005,-0.75500000000005,-0.75500000000005,-0.75500000000005 +3.25,-0.75000000000005,-0.75000000000005,-0.75000000000005,-0.75000000000005,-0.75000000000005,-0.75000000000005,-0.75000000000005 +3.255,-0.74500000000005,-0.74500000000005,-0.74500000000005,-0.74500000000005,-0.74500000000005,-0.74500000000005,-0.74500000000005 +3.26,-0.74000000000005,-0.74000000000005,-0.74000000000005,-0.74000000000005,-0.74000000000005,-0.74000000000005,-0.74000000000005 +3.265,-0.73500000000005,-0.73500000000005,-0.73500000000005,-0.73500000000005,-0.73500000000005,-0.73500000000005,-0.73500000000005 +3.27,-0.73000000000005,-0.73000000000005,-0.73000000000005,-0.73000000000005,-0.73000000000005,-0.73000000000005,-0.73000000000005 +3.275,-0.72500000000005,-0.72500000000005,-0.72500000000005,-0.72500000000005,-0.72500000000005,-0.72500000000005,-0.72500000000005 +3.28,-0.72000000000005,-0.72000000000005,-0.72000000000005,-0.72000000000005,-0.72000000000005,-0.72000000000005,-0.72000000000005 +3.285,-0.71500000000005,-0.71500000000005,-0.71500000000005,-0.71500000000005,-0.71500000000005,-0.71500000000005,-0.71500000000005 +3.29,-0.71000000000005,-0.71000000000005,-0.71000000000005,-0.71000000000005,-0.71000000000005,-0.71000000000005,-0.71000000000005 +3.295,-0.70500000000005,-0.70500000000005,-0.70500000000005,-0.70500000000005,-0.70500000000005,-0.70500000000005,-0.70500000000005 +3.3,-0.70000000000005,-0.70000000000005,-0.70000000000005,-0.70000000000005,-0.70000000000005,-0.70000000000005,-0.70000000000005 +3.305,-0.69500000000005,-0.69500000000005,-0.69500000000005,-0.69500000000005,-0.69500000000005,-0.69500000000005,-0.69500000000005 +3.31,-0.69000000000005,-0.69000000000005,-0.69000000000005,-0.69000000000005,-0.69000000000005,-0.69000000000005,-0.69000000000005 +3.315,-0.68500000000005,-0.68500000000005,-0.68500000000005,-0.68500000000005,-0.68500000000005,-0.68500000000005,-0.68500000000005 +3.32,-0.68000000000005,-0.68000000000005,-0.68000000000005,-0.68000000000005,-0.68000000000005,-0.68000000000005,-0.68000000000005 +3.325,-0.67500000000005,-0.67500000000005,-0.67500000000005,-0.67500000000005,-0.67500000000005,-0.67500000000005,-0.67500000000005 +3.33,-0.67000000000005,-0.67000000000005,-0.67000000000005,-0.67000000000005,-0.67000000000005,-0.67000000000005,-0.67000000000005 +3.335,-0.66500000000005,-0.66500000000005,-0.66500000000005,-0.66500000000005,-0.66500000000005,-0.66500000000005,-0.66500000000005 +3.34,-0.66000000000005,-0.66000000000005,-0.66000000000005,-0.66000000000005,-0.66000000000005,-0.66000000000005,-0.66000000000005 +3.345,-0.65500000000005,-0.65500000000005,-0.65500000000005,-0.65500000000005,-0.65500000000005,-0.65500000000005,-0.65500000000005 +3.35,-0.65000000000005,-0.65000000000005,-0.65000000000005,-0.65000000000005,-0.65000000000005,-0.65000000000005,-0.65000000000005 +3.355,-0.64500000000005,-0.64500000000005,-0.64500000000005,-0.64500000000005,-0.64500000000005,-0.64500000000005,-0.64500000000005 +3.36,-0.64000000000005,-0.64000000000005,-0.64000000000005,-0.64000000000005,-0.64000000000005,-0.64000000000005,-0.64000000000005 +3.365,-0.63500000000005,-0.63500000000005,-0.63500000000005,-0.63500000000005,-0.63500000000005,-0.63500000000005,-0.63500000000005 +3.37,-0.63000000000005,-0.63000000000005,-0.63000000000005,-0.63000000000005,-0.63000000000005,-0.63000000000005,-0.63000000000005 +3.375,-0.62500000000005,-0.62500000000005,-0.62500000000005,-0.62500000000005,-0.62500000000005,-0.62500000000005,-0.62500000000005 +3.38,-0.62000000000005,-0.62000000000005,-0.62000000000005,-0.62000000000005,-0.62000000000005,-0.62000000000005,-0.62000000000005 +3.385,-0.61500000000005,-0.61500000000005,-0.61500000000005,-0.61500000000005,-0.61500000000005,-0.61500000000005,-0.61500000000005 +3.3899999999999,-0.61000000000005,-0.61000000000005,-0.61000000000005,-0.61000000000005,-0.61000000000005,-0.61000000000005,-0.61000000000005 +3.3949999999999,-0.60500000000005,-0.60500000000005,-0.60500000000005,-0.60500000000005,-0.60500000000005,-0.60500000000005,-0.60500000000005 +3.3999999999999,-0.60000000000005,-0.60000000000005,-0.60000000000005,-0.60000000000005,-0.60000000000005,-0.60000000000005,-0.60000000000005 +3.4049999999999,-0.59500000000005,-0.59500000000005,-0.59500000000005,-0.59500000000005,-0.59500000000005,-0.59500000000005,-0.59500000000005 +3.4099999999999,-0.59000000000005,-0.59000000000005,-0.59000000000005,-0.59000000000005,-0.59000000000005,-0.59000000000005,-0.59000000000005 +3.4149999999999,-0.58500000000005,-0.58500000000005,-0.58500000000005,-0.58500000000005,-0.58500000000005,-0.58500000000005,-0.58500000000005 +3.4199999999999,-0.58000000000005,-0.58000000000005,-0.58000000000005,-0.58000000000005,-0.58000000000005,-0.58000000000005,-0.58000000000005 +3.4249999999999,-0.57500000000005,-0.57500000000005,-0.57500000000005,-0.57500000000005,-0.57500000000005,-0.57500000000005,-0.57500000000005 +3.4299999999999,-0.57000000000005,-0.57000000000005,-0.57000000000005,-0.57000000000005,-0.57000000000005,-0.57000000000005,-0.57000000000005 +3.4349999999999,-0.56500000000005,-0.56500000000005,-0.56500000000005,-0.56500000000005,-0.56500000000005,-0.56500000000005,-0.56500000000005 +3.4399999999999,-0.56000000000005,-0.56000000000005,-0.56000000000005,-0.56000000000005,-0.56000000000005,-0.56000000000005,-0.56000000000005 +3.4449999999999,-0.55500000000005,-0.55500000000005,-0.55500000000005,-0.55500000000005,-0.55500000000005,-0.55500000000005,-0.55500000000005 +3.4499999999999,-0.55000000000005,-0.55000000000005,-0.55000000000005,-0.55000000000005,-0.55000000000005,-0.55000000000005,-0.55000000000005 +3.4549999999999,-0.54500000000005,-0.54500000000005,-0.54500000000005,-0.54500000000005,-0.54500000000005,-0.54500000000005,-0.54500000000005 +3.4599999999999,-0.54000000000005,-0.54000000000005,-0.54000000000005,-0.54000000000005,-0.54000000000005,-0.54000000000005,-0.54000000000005 +3.4649999999999,-0.53500000000005,-0.53500000000005,-0.53500000000005,-0.53500000000005,-0.53500000000005,-0.53500000000005,-0.53500000000005 +3.4699999999999,-0.53000000000005,-0.53000000000005,-0.53000000000005,-0.53000000000005,-0.53000000000005,-0.53000000000005,-0.53000000000005 +3.4749999999999,-0.52500000000005,-0.52500000000005,-0.52500000000005,-0.52500000000005,-0.52500000000005,-0.52500000000005,-0.52500000000005 +3.4799999999999,-0.52000000000005,-0.52000000000005,-0.52000000000005,-0.52000000000005,-0.52000000000005,-0.52000000000005,-0.52000000000005 +3.4849999999999,-0.51500000000005,-0.51500000000005,-0.51500000000005,-0.51500000000005,-0.51500000000005,-0.51500000000005,-0.51500000000005 +3.4899999999999,-0.51000000000005,-0.51000000000005,-0.51000000000005,-0.51000000000005,-0.51000000000005,-0.51000000000005,-0.51000000000005 +3.4949999999999,-0.50500000000005,-0.50500000000005,-0.50500000000005,-0.50500000000005,-0.50500000000005,-0.50500000000005,-0.50500000000005 +3.4999999999999,-0.50000000000005,-0.50000000000005,-0.50000000000005,-0.50000000000005,-0.50000000000005,-0.50000000000005,-0.50000000000005 +3.5049999999999,-0.49500000000005,-0.49500000000005,-0.49500000000005,-0.49500000000005,-0.49500000000005,-0.49500000000005,-0.49500000000005 +3.5099999999999,-0.49000000000005,-0.49000000000005,-0.49000000000005,-0.49000000000005,-0.49000000000005,-0.49000000000005,-0.49000000000005 +3.5149999999999,-0.48500000000005,-0.48500000000005,-0.48500000000005,-0.48500000000005,-0.48500000000005,-0.48500000000005,-0.48500000000005 +3.5199999999999,-0.48000000000005,-0.48000000000005,-0.48000000000005,-0.48000000000005,-0.48000000000005,-0.48000000000005,-0.48000000000005 +3.5249999999999,-0.47500000000005,-0.47500000000005,-0.47500000000005,-0.47500000000005,-0.47500000000005,-0.47500000000005,-0.47500000000005 +3.5299999999999,-0.47000000000005,-0.47000000000005,-0.47000000000005,-0.47000000000005,-0.47000000000005,-0.47000000000005,-0.47000000000005 +3.5349999999999,-0.46500000000005,-0.46500000000005,-0.46500000000005,-0.46500000000005,-0.46500000000005,-0.46500000000005,-0.46500000000005 +3.5399999999999,-0.46000000000005,-0.46000000000005,-0.46000000000005,-0.46000000000005,-0.46000000000005,-0.46000000000005,-0.46000000000005 +3.5449999999999,-0.45500000000005,-0.45500000000005,-0.45500000000005,-0.45500000000005,-0.45500000000005,-0.45500000000005,-0.45500000000005 +3.5499999999999,-0.45000000000005,-0.45000000000005,-0.45000000000005,-0.45000000000005,-0.45000000000005,-0.45000000000005,-0.45000000000005 +3.5549999999999,-0.44500000000005,-0.44500000000005,-0.44500000000005,-0.44500000000005,-0.44500000000005,-0.44500000000005,-0.44500000000005 +3.5599999999999,-0.44000000000005,-0.44000000000005,-0.44000000000005,-0.44000000000005,-0.44000000000005,-0.44000000000005,-0.44000000000005 +3.5649999999999,-0.43500000000005,-0.43500000000005,-0.43500000000005,-0.43500000000005,-0.43500000000005,-0.43500000000005,-0.43500000000005 +3.5699999999999,-0.43000000000005,-0.43000000000005,-0.43000000000005,-0.43000000000005,-0.43000000000005,-0.43000000000005,-0.43000000000005 +3.5749999999999,-0.42500000000005,-0.42500000000005,-0.42500000000005,-0.42500000000005,-0.42500000000005,-0.42500000000005,-0.42500000000005 +3.5799999999999,-0.42000000000005,-0.42000000000005,-0.42000000000005,-0.42000000000005,-0.42000000000005,-0.42000000000005,-0.42000000000005 +3.5849999999999,-0.41500000000005,-0.41500000000005,-0.41500000000005,-0.41500000000005,-0.41500000000005,-0.41500000000005,-0.41500000000005 +3.5899999999999,-0.41000000000005,-0.41000000000005,-0.41000000000005,-0.41000000000005,-0.41000000000005,-0.41000000000005,-0.41000000000005 +3.5949999999999,-0.40500000000005,-0.40500000000005,-0.40500000000005,-0.40500000000005,-0.40500000000005,-0.40500000000005,-0.40500000000005 +3.5999999999999,-0.40000000000005,-0.40000000000005,-0.40000000000005,-0.40000000000005,-0.40000000000005,-0.40000000000005,-0.40000000000005 +3.6049999999999,-0.39500000000005,-0.39500000000005,-0.39500000000005,-0.39500000000005,-0.39500000000005,-0.39500000000005,-0.39500000000005 +3.6099999999999,-0.39000000000005,-0.39000000000005,-0.39000000000005,-0.39000000000005,-0.39000000000005,-0.39000000000005,-0.39000000000005 +3.6149999999999,-0.38500000000005,-0.38500000000005,-0.38500000000005,-0.38500000000005,-0.38500000000005,-0.38500000000005,-0.38500000000005 +3.6199999999999,-0.38000000000005,-0.38000000000005,-0.38000000000005,-0.38000000000005,-0.38000000000005,-0.38000000000005,-0.38000000000005 +3.6249999999999,-0.37500000000006,-0.37500000000006,-0.37500000000006,-0.37500000000006,-0.37500000000006,-0.37500000000006,-0.37500000000006 +3.6299999999999,-0.37000000000006,-0.37000000000006,-0.37000000000006,-0.37000000000006,-0.37000000000006,-0.37000000000006,-0.37000000000006 +3.6349999999999,-0.36500000000006,-0.36500000000006,-0.36500000000006,-0.36500000000006,-0.36500000000006,-0.36500000000006,-0.36500000000006 +3.6399999999999,-0.36000000000006,-0.36000000000006,-0.36000000000006,-0.36000000000006,-0.36000000000006,-0.36000000000006,-0.36000000000006 +3.6449999999999,-0.35500000000006,-0.35500000000006,-0.35500000000006,-0.35500000000006,-0.35500000000006,-0.35500000000006,-0.35500000000006 +3.6499999999999,-0.35000000000006,-0.35000000000006,-0.35000000000006,-0.35000000000006,-0.35000000000006,-0.35000000000006,-0.35000000000006 +3.6549999999999,-0.34500000000006,-0.34500000000006,-0.34500000000006,-0.34500000000006,-0.34500000000006,-0.34500000000006,-0.34500000000006 +3.6599999999999,-0.34000000000006,-0.34000000000006,-0.34000000000006,-0.34000000000006,-0.34000000000006,-0.34000000000006,-0.34000000000006 +3.6649999999999,-0.33500000000006,-0.33500000000006,-0.33500000000006,-0.33500000000006,-0.33500000000006,-0.33500000000006,-0.33500000000006 +3.6699999999999,-0.33000000000006,-0.33000000000006,-0.33000000000006,-0.33000000000006,-0.33000000000006,-0.33000000000006,-0.33000000000006 +3.6749999999999,-0.32500000000006,-0.32500000000006,-0.32500000000006,-0.32500000000006,-0.32500000000006,-0.32500000000006,-0.32500000000006 +3.6799999999999,-0.32000000000006,-0.32000000000006,-0.32000000000006,-0.32000000000006,-0.32000000000006,-0.32000000000006,-0.32000000000006 +3.6849999999999,-0.31500000000006,-0.31500000000006,-0.31500000000006,-0.31500000000006,-0.31500000000006,-0.31500000000006,-0.31500000000006 +3.6899999999999,-0.31000000000006,-0.31000000000006,-0.31000000000006,-0.31000000000006,-0.31000000000006,-0.31000000000006,-0.31000000000006 +3.6949999999999,-0.30500000000006,-0.30500000000006,-0.30500000000006,-0.30500000000006,-0.30500000000006,-0.30500000000006,-0.30500000000006 +3.6999999999999,-0.30000000000006,-0.30000000000006,-0.30000000000006,-0.30000000000006,-0.30000000000006,-0.30000000000006,-0.30000000000006 +3.7049999999999,-0.29500000000006,-0.29500000000006,-0.29500000000006,-0.29500000000006,-0.29500000000006,-0.29500000000006,-0.29500000000006 +3.7099999999999,-0.29000000000006,-0.29000000000006,-0.29000000000006,-0.29000000000006,-0.29000000000006,-0.29000000000006,-0.29000000000006 +3.7149999999999,-0.28500000000006,-0.28500000000006,-0.28500000000006,-0.28500000000006,-0.28500000000006,-0.28500000000006,-0.28500000000006 +3.7199999999999,-0.28000000000006,-0.28000000000006,-0.28000000000006,-0.28000000000006,-0.28000000000006,-0.28000000000006,-0.28000000000006 +3.7249999999999,-0.27500000000006,-0.27500000000006,-0.27500000000006,-0.27500000000006,-0.27500000000006,-0.27500000000006,-0.27500000000006 +3.7299999999999,-0.27000000000006,-0.27000000000006,-0.27000000000006,-0.27000000000006,-0.27000000000006,-0.27000000000006,-0.27000000000006 +3.7349999999999,-0.26500000000006,-0.26500000000006,-0.26500000000006,-0.26500000000006,-0.26500000000006,-0.26500000000006,-0.26500000000006 +3.7399999999999,-0.26000000000006,-0.26000000000006,-0.26000000000006,-0.26000000000006,-0.26000000000006,-0.26000000000006,-0.26000000000006 +3.7449999999999,-0.25500000000006,-0.25500000000006,-0.25500000000006,-0.25500000000006,-0.25500000000006,-0.25500000000006,-0.25500000000006 +3.7499999999999,-0.25000000000006,-0.25000000000006,-0.25000000000006,-0.25000000000006,-0.25000000000006,-0.25000000000006,-0.25000000000006 +3.7549999999999,-0.24500000000006,-0.24500000000006,-0.24500000000006,-0.24500000000006,-0.24500000000006,-0.24500000000006,-0.24500000000006 +3.7599999999999,-0.24000000000006,-0.24000000000006,-0.24000000000006,-0.24000000000006,-0.24000000000006,-0.24000000000006,-0.24000000000006 +3.7649999999999,-0.23500000000006,-0.23500000000006,-0.23500000000006,-0.23500000000006,-0.23500000000006,-0.23500000000006,-0.23500000000006 +3.7699999999999,-0.23000000000006,-0.23000000000006,-0.23000000000006,-0.23000000000006,-0.23000000000006,-0.23000000000006,-0.23000000000006 +3.7749999999999,-0.22500000000006,-0.22500000000006,-0.22500000000006,-0.22500000000006,-0.22500000000006,-0.22500000000006,-0.22500000000006 +3.7799999999999,-0.22000000000006,-0.22000000000006,-0.22000000000006,-0.22000000000006,-0.22000000000006,-0.22000000000006,-0.22000000000006 +3.7849999999999,-0.21500000000006,-0.21500000000006,-0.21500000000006,-0.21500000000006,-0.21500000000006,-0.21500000000006,-0.21500000000006 +3.7899999999999,-0.21000000000006,-0.21000000000006,-0.21000000000006,-0.21000000000006,-0.21000000000006,-0.21000000000006,-0.21000000000006 +3.7949999999999,-0.20500000000006,-0.20500000000006,-0.20500000000006,-0.20500000000006,-0.20500000000006,-0.20500000000006,-0.20500000000006 +3.7999999999999,-0.20000000000006,-0.20000000000006,-0.20000000000006,-0.20000000000006,-0.20000000000006,-0.20000000000006,-0.20000000000006 +3.8049999999999,-0.19500000000006,-0.19500000000006,-0.19500000000006,-0.19500000000006,-0.19500000000006,-0.19500000000006,-0.19500000000006 +3.8099999999999,-0.19000000000006,-0.19000000000006,-0.19000000000006,-0.19000000000006,-0.19000000000006,-0.19000000000006,-0.19000000000006 +3.8149999999999,-0.18500000000006,-0.18500000000006,-0.18500000000006,-0.18500000000006,-0.18500000000006,-0.18500000000006,-0.18500000000006 +3.8199999999999,-0.18000000000006,-0.18000000000006,-0.18000000000006,-0.18000000000006,-0.18000000000006,-0.18000000000006,-0.18000000000006 +3.8249999999999,-0.17500000000006,-0.17500000000006,-0.17500000000006,-0.17500000000006,-0.17500000000006,-0.17500000000006,-0.17500000000006 +3.8299999999999,-0.17000000000006,-0.17000000000006,-0.17000000000006,-0.17000000000006,-0.17000000000006,-0.17000000000006,-0.17000000000006 +3.8349999999999,-0.16500000000006,-0.16500000000006,-0.16500000000006,-0.16500000000006,-0.16500000000006,-0.16500000000006,-0.16500000000006 +3.8399999999999,-0.16000000000006,-0.16000000000006,-0.16000000000006,-0.16000000000006,-0.16000000000006,-0.16000000000006,-0.16000000000006 +3.8449999999999,-0.15500000000006,-0.15500000000006,-0.15500000000006,-0.15500000000006,-0.15500000000006,-0.15500000000006,-0.15500000000006 +3.8499999999999,-0.15000000000006,-0.15000000000006,-0.15000000000006,-0.15000000000006,-0.15000000000006,-0.15000000000006,-0.15000000000006 +3.8549999999999,-0.14500000000006,-0.14500000000006,-0.14500000000006,-0.14500000000006,-0.14500000000006,-0.14500000000006,-0.14500000000006 +3.8599999999999,-0.14000000000006,-0.14000000000006,-0.14000000000006,-0.14000000000006,-0.14000000000006,-0.14000000000006,-0.14000000000006 +3.8649999999999,-0.13500000000006,-0.13500000000006,-0.13500000000006,-0.13500000000006,-0.13500000000006,-0.13500000000006,-0.13500000000006 +3.8699999999999,-0.13000000000006,-0.13000000000006,-0.13000000000006,-0.13000000000006,-0.13000000000006,-0.13000000000006,-0.13000000000006 +3.8749999999999,-0.12500000000006,-0.12500000000006,-0.12500000000006,-0.12500000000006,-0.12500000000006,-0.12500000000006,-0.12500000000006 +3.8799999999999,-0.12000000000006,-0.12000000000006,-0.12000000000006,-0.12000000000006,-0.12000000000006,-0.12000000000006,-0.12000000000006 +3.8849999999999,-0.11500000000006,-0.11500000000006,-0.11500000000006,-0.11500000000006,-0.11500000000006,-0.11500000000006,-0.11500000000006 +3.8899999999999,-0.11000000000006,-0.11000000000006,-0.11000000000006,-0.11000000000006,-0.11000000000006,-0.11000000000006,-0.11000000000006 +3.8949999999999,-0.10500000000006,-0.10500000000006,-0.10500000000006,-0.10500000000006,-0.10500000000006,-0.10500000000006,-0.10500000000006 +3.8999999999999,-0.10000000000006,-0.10000000000006,-0.10000000000006,-0.10000000000006,-0.10000000000006,-0.10000000000006,-0.10000000000006 +3.9049999999999,-0.095000000000061,-0.095000000000061,-0.095000000000061,-0.095000000000061,-0.095000000000061,-0.095000000000061,-0.095000000000061 +3.9099999999999,-0.090000000000061,-0.090000000000061,-0.090000000000061,-0.090000000000061,-0.090000000000061,-0.090000000000061,-0.090000000000061 +3.9149999999999,-0.085000000000061,-0.085000000000061,-0.085000000000061,-0.085000000000061,-0.085000000000061,-0.085000000000061,-0.085000000000061 +3.9199999999999,-0.080000000000061,-0.080000000000061,-0.080000000000061,-0.080000000000061,-0.080000000000061,-0.080000000000061,-0.080000000000061 +3.9249999999999,-0.075000000000061,-0.075000000000061,-0.075000000000061,-0.075000000000061,-0.075000000000061,-0.075000000000061,-0.075000000000061 +3.9299999999999,-0.070000000000062,-0.070000000000062,-0.070000000000062,-0.070000000000062,-0.070000000000062,-0.070000000000062,-0.070000000000062 +3.9349999999999,-0.065000000000062,-0.065000000000062,-0.065000000000062,-0.065000000000062,-0.065000000000062,-0.065000000000062,-0.065000000000062 +3.9399999999999,-0.060000000000062,-0.060000000000062,-0.060000000000062,-0.060000000000062,-0.060000000000062,-0.060000000000062,-0.060000000000062 +3.9449999999999,-0.055000000000062,-0.055000000000062,-0.055000000000062,-0.055000000000062,-0.055000000000062,-0.055000000000062,-0.055000000000062 +3.9499999999999,-0.050000000000062,-0.050000000000062,-0.050000000000062,-0.050000000000062,-0.050000000000062,-0.050000000000062,-0.050000000000062 +3.9549999999999,-0.045000000000062,-0.045000000000062,-0.045000000000062,-0.045000000000062,-0.045000000000062,-0.045000000000062,-0.045000000000062 +3.9599999999999,-0.040000000000062,-0.040000000000062,-0.040000000000062,-0.040000000000062,-0.040000000000062,-0.040000000000062,-0.040000000000062 +3.9649999999999,-0.035000000000062,-0.035000000000062,-0.035000000000062,-0.035000000000062,-0.035000000000062,-0.035000000000062,-0.035000000000062 +3.9699999999999,-0.030000000000062,-0.030000000000062,-0.030000000000062,-0.030000000000062,-0.030000000000062,-0.030000000000062,-0.030000000000062 +3.9749999999999,-0.025000000000063,-0.025000000000063,-0.025000000000063,-0.025000000000063,-0.025000000000063,-0.025000000000063,-0.025000000000063 +3.9799999999999,-0.020000000000063,-0.020000000000063,-0.020000000000063,-0.020000000000063,-0.020000000000063,-0.020000000000063,-0.020000000000063 +3.9849999999999,-0.015000000000063,-0.015000000000063,-0.015000000000063,-0.015000000000063,-0.015000000000063,-0.015000000000063,-0.015000000000063 +3.9899999999999,-0.010000000000063,-0.010000000000063,-0.010000000000063,-0.010000000000063,-0.010000000000063,-0.010000000000063,-0.010000000000063 +3.9949999999999,-0.005000000000063,-0.005000000000063,-0.005000000000063,-0.005000000000063,-0.005000000000063,-0.005000000000063,-0.005000000000063 +3.9999999999999,-6.3058933075233e-14,-6.3058933075233e-14,-6.3058933075233e-14,-6.3058933075233e-14,-6.3058933075233e-14,-6.3058933075233e-14,-6.299995247705e-14 diff --git a/test/tests/auxkernels/materialrealcmm/spring_static.i b/test/tests/auxkernels/materialrealcmm/spring_static.i new file mode 100644 index 0000000000..92de11dc4a --- /dev/null +++ b/test/tests/auxkernels/materialrealcmm/spring_static.i @@ -0,0 +1,303 @@ +# Test for MaterialRealCMMAux auxkernel +# A sawtooth force history is applied at one node of a spring element in each +# of the 6 degrees of freedom. The displacements and rotations are calculated +# at this node. The other node is fixed in all directions. +# This is a static analysis and therefore, the inertia kernels are ommitted. +# The global axial force in the spring is also calculated using the MaterialRealCMMAux +# auxkernel. + +[Mesh] + type = GeneratedMesh + xmin = 0 + xmax = 1 + nx = 1 + dim = 1 + displacements = 'disp_x disp_y disp_z' +[] + +[Variables] + [disp_x] + order = FIRST + family = LAGRANGE + [] + [disp_y] + order = FIRST + family = LAGRANGE + [] + [disp_z] + order = FIRST + family = LAGRANGE + [] + [rot_x] + order = FIRST + family = LAGRANGE + [] + [rot_y] + order = FIRST + family = LAGRANGE + [] + [rot_z] + order = FIRST + family = LAGRANGE + [] +[] + +[AuxVariables] + [global_force_x] + block = '0' + order = CONSTANT + family = MONOMIAL + [] +[] + +[AuxKernels] + [global_force_x] + type = MaterialRealCMMAux + property = global_forces + row = 0 + column = 0 + variable = global_force_x + block = '0' + [] +[] + +[Kernels] + [spring_disp_x] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 0 + variable = disp_x + [] + [spring_disp_y] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 1 + variable = disp_y + [] + [spring_disp_z] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 2 + variable = disp_z + [] + [spring_rot_x] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 3 + variable = rot_x + [] + [spring_rot_y] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 4 + variable = rot_y + [] + [spring_rot_z] + type = StressDivergenceSpring + block = '0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + component = 5 + variable = rot_z + [] +[] + +[BCs] + [fixx1] + type = DirichletBC + variable = disp_x + boundary = left + value = 0.0 + [] + [fixy1] + type = DirichletBC + variable = disp_y + boundary = left + value = 0.0 + [] + [fixz1] + type = DirichletBC + variable = disp_z + boundary = left + value = 0.0 + [] + [fixr1] + type = DirichletBC + variable = rot_x + boundary = left + value = 0.0 + [] + [fixr2] + type = DirichletBC + variable = rot_y + boundary = left + value = 0.0 + [] + [fixr3] + type = DirichletBC + variable = rot_z + boundary = left + value = 0.0 + [] +[] + +[NodalKernels] + [force_x] + type = UserForcingFunctionNodalKernel + variable = disp_x + boundary = right + function = force_x + [] + [force_y] + type = UserForcingFunctionNodalKernel + variable = disp_y + boundary = right + function = force_y + [] + [force_z] + type = UserForcingFunctionNodalKernel + variable = disp_z + boundary = right + function = force_z + [] + [moment_x] + type = UserForcingFunctionNodalKernel + variable = rot_x + boundary = right + function = moment_x + [] + [moment_y] + type = UserForcingFunctionNodalKernel + variable = rot_y + boundary = right + function = moment_y + [] + [moment_z] + type = UserForcingFunctionNodalKernel + variable = rot_z + boundary = right + function = moment_z + [] +[] + +[Functions] + [force_x] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 1.0 0.0 -1.0 0.0' # force + [] + [force_y] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 2.0 0.0 -2.0 0.0' # force + [] + [force_z] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 3.0 0.0 -3.0 0.0' # force + [] + [moment_x] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 10.0 0.0 -10.0 0.0' # moment + [] + [moment_y] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 20.0 0.0 -20.0 0.0' # moment + [] + [moment_z] + type = PiecewiseLinear + x = '0.0 1.0 2.0 3.0 4.0' # time + y = '0.0 30.0 0.0 -30.0 0.0' # force + [] +[] + +[Materials] + [linear_spring_test] + type = LinearSpring + block = 0 + y_orientation = '0.0 1.0 0.0' + displacements = 'disp_x disp_y disp_z' + rotations = 'rot_x rot_y rot_z' + kx = 1.0 + ky = 2.0 + kz = 3.0 + krx = 10.0 + kry = 20.0 + krz = 30.0 + [] +[] + +[Preconditioning] + [smp] + type = SMP + full = true + [] +[] + +[Executioner] + type = Transient + solve_type = NEWTON + line_search = none + nl_rel_tol = 1e-8 + nl_abs_tol = 1e-8 + start_time = 0.0 + end_time = 4.0 + dt = 0.005 + dtmin = 1e-4 + timestep_tolerance = 1e-6 +[] + +[Postprocessors] + [disp_x] + type = PointValue + point = '1.0 0.0 0.0' + variable = disp_x + [] + [disp_y] + type = PointValue + point = '1.0 0.0 0.0' + variable = disp_y + [] + [disp_z] + type = PointValue + point = '1.0 0.0 0.0' + variable = disp_z + [] + [rot_x] + type = PointValue + point = '1.0 0.0 0.0' + variable = rot_x + [] + [rot_y] + type = PointValue + point = '1.0 0.0 0.0' + variable = rot_y + [] + [rot_z] + type = PointValue + point = '1.0 0.0 0.0' + variable = rot_z + [] + [gf_x] + type = PointValue + point = '0.5 0 0' + variable = global_force_x + [] +[] + +[Outputs] + exodus = false + csv = true +[] diff --git a/test/tests/auxkernels/materialrealcmm/tests b/test/tests/auxkernels/materialrealcmm/tests new file mode 100644 index 0000000000..6d8f45de9d --- /dev/null +++ b/test/tests/auxkernels/materialrealcmm/tests @@ -0,0 +1,11 @@ +[Tests] + [./cmm_spring] + type = CSVDiff + input = spring_static.i + csvdiff = spring_static_out.csv + + design = "MaterialRealCMMAux.md" + issues = "#374" + requirement = "The MaterialRealCMMAux AuxKernel shall accurately calculate the AuxVariable from a ColumnMajorMatrix material property." + [../] +[] From e66bf463307d91020af4e6aad737102f48a3a9ad Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 17:41:14 -0600 Subject: [PATCH 4/6] Adding doco for MaterialRealCMMAux. Closes #376. --- .../source/auxkernels/MaterialRealCMMAux.md | 32 +++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100644 doc/content/source/auxkernels/MaterialRealCMMAux.md diff --git a/doc/content/source/auxkernels/MaterialRealCMMAux.md b/doc/content/source/auxkernels/MaterialRealCMMAux.md new file mode 100644 index 0000000000..26e9041d60 --- /dev/null +++ b/doc/content/source/auxkernels/MaterialRealCMMAux.md @@ -0,0 +1,32 @@ +# MaterialRealCMMAux + +!syntax description /AuxKernels/MaterialRealCMMAux + +## Description + +`MaterialRealCMMAux` can be used to populate an `AuxVariable` with an element of +a material property that is of type, ColumnMajorMatrix (CMM). In MASTODON, a +common use case is accessing the deformation or force CMMs from link elements +when the [`LinearSpring`](source/materials/LinearSpring.md) material or any of +the seismic isolator materials (e.g., [`ComputeFPIsolatorElasticity`](source/materials/ComputeFPIsolatorElasticity.md)) is used. This `AuxKernel` and the corresponding +`AuxVariable` should be block-restricted to blocks where the material property +given by the `property` input parameter below is defined. + +For example, the following input file syntax will create the AuxVariable field +called `global_force_x` with the X direction force in the spring element modeled +using the LinearSpring material. The AuxKernel is defined as follows: + +!listing test/tests/auxkernels/materialrealcmm/spring_static.i block=AuxKernels + +In this input, the AuxKernel takes the element (0, 0) from the CMM material +property named, "ls_global_forces". This material property name refers to the +global forces in the linear spring material. (See LinearSpring.C for other +material properties.) In this example, the linear spring material is only defined +in block '0' and therefore this AuxKernel and AuxVariable are also restricted to +this block. + +!syntax parameters /AuxKernels/MaterialRealCMMAux + +!syntax inputs /AuxKernels/MaterialRealCMMAux + +!syntax children /AuxKernels/MaterialRealCMMAux From 1fd37bdcc8a9edac481470e7ab1295b82d01074d Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 17:51:45 -0600 Subject: [PATCH 5/6] Addressing @somu15 comments. --- .../source/auxkernels/MaterialRealCMMAux.md | 21 ++++--------------- 1 file changed, 4 insertions(+), 17 deletions(-) diff --git a/doc/content/source/auxkernels/MaterialRealCMMAux.md b/doc/content/source/auxkernels/MaterialRealCMMAux.md index 26e9041d60..34ca04fad6 100644 --- a/doc/content/source/auxkernels/MaterialRealCMMAux.md +++ b/doc/content/source/auxkernels/MaterialRealCMMAux.md @@ -4,26 +4,13 @@ ## Description -`MaterialRealCMMAux` can be used to populate an `AuxVariable` with an element of -a material property that is of type, ColumnMajorMatrix (CMM). In MASTODON, a -common use case is accessing the deformation or force CMMs from link elements -when the [`LinearSpring`](source/materials/LinearSpring.md) material or any of -the seismic isolator materials (e.g., [`ComputeFPIsolatorElasticity`](source/materials/ComputeFPIsolatorElasticity.md)) is used. This `AuxKernel` and the corresponding -`AuxVariable` should be block-restricted to blocks where the material property -given by the `property` input parameter below is defined. - -For example, the following input file syntax will create the AuxVariable field -called `global_force_x` with the X direction force in the spring element modeled -using the LinearSpring material. The AuxKernel is defined as follows: +`MaterialRealCMMAux` can be used to populate an `AuxVariable` with an element of a material property that is of type, ColumnMajorMatrix (CMM). In MASTODON, a common use case is accessing the deformation or force CMMs from link elements when the [`LinearSpring`](source/materials/LinearSpring.md) material or any of the seismic isolator materials (e.g., [`ComputeFPIsolatorElasticity`](source/materials/ComputeFPIsolatorElasticity.md)) is used. This `AuxKernel` and the corresponding `AuxVariable` should be block-restricted to blocks where the material property given by the `property` input parameter below is defined. + +For example, the following input file syntax will create the AuxVariable field called `global_force_x` with the X direction force in the spring element modeled using the LinearSpring material. Users can look at the LinearSpring.C file in the source code to see that this material has several material properties such as "global_forces", "global_moments", etc. In the case of the linear spring material, the "global_forces" material property is a CMM of size 6x1 and stores the forces in the spring in the global co-ordinate system. The AuxKernel to make this calculation is listed below. !listing test/tests/auxkernels/materialrealcmm/spring_static.i block=AuxKernels -In this input, the AuxKernel takes the element (0, 0) from the CMM material -property named, "ls_global_forces". This material property name refers to the -global forces in the linear spring material. (See LinearSpring.C for other -material properties.) In this example, the linear spring material is only defined -in block '0' and therefore this AuxKernel and AuxVariable are also restricted to -this block. +In this input, the AuxKernel takes the element (0, 0) from the CMM material property named, "global_forces" in the linear spring material. (See LinearSpring.C for other material properties.) In this example, the linear spring material is only defined in block '0' and therefore this AuxKernel and AuxVariable are also restricted to this block. The element (0, 0) of the material property "global_forces" corresponds to the axial force in the spring. Currently, in order to know what each of the elements in a material property CMM correspond to, users will either have to go through the source code or reach out to one of the MASTODON developers [here](help/contact_us.md). !syntax parameters /AuxKernels/MaterialRealCMMAux From 217598fcb8823cacee0802aab4ecf2fdd173b756 Mon Sep 17 00:00:00 2001 From: Chandu Bolisetti Date: Fri, 11 Jun 2021 18:08:51 -0600 Subject: [PATCH 6/6] Fixing material parameter description in computeLRIslatorElasticity. Closes #326. --- src/materials/ComputeLRIsolatorElasticity.C | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/materials/ComputeLRIsolatorElasticity.C b/src/materials/ComputeLRIsolatorElasticity.C index b4a6939878..bb7da1c409 100644 --- a/src/materials/ComputeLRIsolatorElasticity.C +++ b/src/materials/ComputeLRIsolatorElasticity.C @@ -52,8 +52,10 @@ ComputeLRIsolatorElasticity::validParams() "Switch for modeling strength degradation due to lead core heating."); // Strength degradation // due to heating // Material properties - params.addRequiredParam("fy", "Yield stress of the bearing."); - params.addRequiredParam("alpha", "Yield displacement of the bearing."); + params.addRequiredParam("fy", "Yield strength of the bearing."); + params.addRequiredParam("alpha", + "Ratio of post-yield shear stiffness to the initial elastic shear " + "stiffness of the bearing. This is dimensionless"); params.addRequiredParam("G_rubber", "Shear modulus of rubber."); params.addRequiredParam("K_rubber", "Bulk modulus of rubber."); params.addRequiredParam("D1", "Diameter of lead core.");