Changeset 9eb7580
- Timestamp:
- Dec 4, 2010, 11:54:32 PM (14 years ago)
- Branches:
- Action_Thermostats, Add_AtomRandomPerturbation, Add_FitFragmentPartialChargesAction, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Added_ParseSaveFragmentResults, AddingActions_SaveParseParticleParameters, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_ParticleName_to_Atom, Adding_StructOpt_integration_tests, AtomFragments, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.5.4, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, CombiningParticlePotentialParsing, Combining_Subpackages, Debian_Package_split, Debian_package_split_molecuildergui_only, Disabling_MemDebug, Docu_Python_wait, EmpiricalPotential_contain_HomologyGraph, EmpiricalPotential_contain_HomologyGraph_documentation, Enable_parallel_make_install, Enhance_userguide, Enhanced_StructuralOptimization, Enhanced_StructuralOptimization_continued, Example_ManyWaysToTranslateAtom, Exclude_Hydrogens_annealWithBondGraph, FitPartialCharges_GlobalError, Fix_BoundInBox_CenterInBox_MoleculeActions, Fix_ChargeSampling_PBC, Fix_ChronosMutex, Fix_FitPartialCharges, Fix_FitPotential_needs_atomicnumbers, Fix_ForceAnnealing, Fix_IndependentFragmentGrids, Fix_ParseParticles, Fix_ParseParticles_split_forward_backward_Actions, Fix_PopActions, Fix_QtFragmentList_sorted_selection, Fix_Restrictedkeyset_FragmentMolecule, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, Fix_fitting_potentials, Fixes, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, FragmentAction_writes_AtomFragments, FragmentMolecule_checks_bonddegrees, GeometryObjects, Gui_Fixes, Gui_displays_atomic_force_velocity, ImplicitCharges, IndependentFragmentGrids, IndependentFragmentGrids_IndividualZeroInstances, IndependentFragmentGrids_IntegrationTest, IndependentFragmentGrids_Sole_NN_Calculation, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, MoreRobust_FragmentAutomation, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PdbParser_setsAtomName, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, Rewrite_FitPartialCharges, RotateToPrincipalAxisSystem_UndoRedo, SaturateAtoms_findBestMatching, SaturateAtoms_singleDegree, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, Switchable_LogView, ThirdParty_MPQC_rebuilt_buildsystem, TrajectoryDependenant_MaxOrder, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, TremoloParser_setsAtomName, Ubuntu_1604_changes, stable
- Children:
- 3bc926
- Parents:
- cfc53b
- git-author:
- Frederik Heber <heber@…> (11/15/10 10:29:25)
- git-committer:
- Frederik Heber <heber@…> (12/04/10 23:54:32)
- Location:
- src
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/AnalysisAction/PrincipalAxisSystemAction.cpp
rcfc53b r9eb7580 48 48 49 49 // reset inertia tensor 50 InertiaTensor. zero();50 InertiaTensor.setZero(); 51 51 52 52 // sum up inertia tensor -
src/Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.cpp
rcfc53b r9eb7580 54 54 55 55 // reset inertia tensor 56 InertiaTensor. zero();56 InertiaTensor.setZero(); 57 57 58 58 // sum up inertia tensor … … 105 105 // summing anew for debugging (resulting matrix has to be diagonal!) 106 106 // reset inertia tensor 107 InertiaTensor. zero();107 InertiaTensor.setZero(); 108 108 109 109 // sum up inertia tensor -
src/Actions/SelectionAction/AllAtomsInsideCuboidAction.cpp
rcfc53b r9eb7580 50 50 getParametersfromValueStorage(); 51 51 52 RotationMatrix. rotation(params.Xangle, params.Yangle, params.Zangle);52 RotationMatrix.setRotation(params.Xangle, params.Yangle, params.Zangle); 53 53 54 54 DoLog(1) && (Log() << Verbose(1) << "Selecting all atoms inside a rotated " << RotationMatrix << " cuboid at " << params.position << " and extension of " << params.extension << "." << endl); … … 72 72 Matrix RotationMatrix; 73 73 74 RotationMatrix. rotation(state->params.Xangle, state->params.Yangle, state->params.Zangle);74 RotationMatrix.setRotation(state->params.Xangle, state->params.Yangle, state->params.Zangle); 75 75 Shape s = translate(transform(stretch(Cuboid(),state->params.extension),RotationMatrix),state->params.position); 76 76 World::getInstance().selectAllAtoms(AtomByShape(s)); -
src/Actions/SelectionAction/NotAllAtomsInsideCuboidAction.cpp
rcfc53b r9eb7580 49 49 // obtain information 50 50 getParametersfromValueStorage(); 51 RotationMatrix. rotation(params.Xangle, params.Yangle, params.Zangle);51 RotationMatrix.setRotation(params.Xangle, params.Yangle, params.Zangle); 52 52 53 53 DoLog(1) && (Log() << Verbose(1) << "Unselecting all atoms inside a rotated " << RotationMatrix << " cuboid at " << params.position << " and extension of " << params.extension << "." << endl); … … 72 72 Matrix RotationMatrix; 73 73 74 RotationMatrix. rotation(state->params.Xangle, state->params.Yangle, state->params.Zangle);74 RotationMatrix.setRotation(state->params.Xangle, state->params.Yangle, state->params.Zangle); 75 75 Shape s = translate(transform(stretch(Cuboid(),state->params.extension),RotationMatrix),state->params.position); 76 76 World::getInstance().unselectAllAtoms(AtomByShape(s)); -
src/Box.cpp
rcfc53b r9eb7580 39 39 { 40 40 M= new Matrix(); 41 M-> one();41 M->setIdentity(); 42 42 Minv = new Matrix(); 43 Minv-> one();43 Minv->setIdentity(); 44 44 conditions.resize(3); 45 45 conditions[0] = conditions[1] = conditions[2] = Wrap; … … 253 253 void Box::setCuboid(const Vector &endpoint){ 254 254 ASSERT(endpoint[0]>0 && endpoint[1]>0 && endpoint[2]>0,"Vector does not define a full cuboid"); 255 M-> one();255 M->setIdentity(); 256 256 M->diagonal()=endpoint; 257 257 Vector &dinv = Minv->diagonal(); -
src/LinearAlgebra/Matrix.cpp
rcfc53b r9eb7580 106 106 } 107 107 108 void Matrix:: one(){108 void Matrix::setIdentity(){ 109 109 for(int i=NDIM;i--;){ 110 110 for(int j=NDIM;j--;){ … … 114 114 } 115 115 116 void Matrix:: zero(){116 void Matrix::setZero(){ 117 117 for(int i=NDIM;i--;){ 118 118 for(int j=NDIM;j--;){ … … 122 122 } 123 123 124 void Matrix:: rotation(const double x, const double y, const double z)124 void Matrix::setRotation(const double x, const double y, const double z) 125 125 { 126 126 set(0,0, cos(y)*cos(z)); -
src/LinearAlgebra/Matrix.hpp
rcfc53b r9eb7580 48 48 * Set the matrix to a unit matrix. 49 49 */ 50 void one();50 void setIdentity(); 51 51 52 52 /** 53 53 * Set all matrix entries to zero. 54 54 */ 55 void zero();55 void setZero(); 56 56 57 57 /** … … 60 60 * with given angles. 61 61 */ 62 void rotation(const double x, const double y, const double z);62 void setRotation(const double x, const double y, const double z); 63 63 64 64 /** -
src/UIElements/Qt4/Pipe/BoxQtQueryPipe.cpp
rcfc53b r9eb7580 34 34 { 35 35 tmpM = new Matrix(); 36 tmpM-> zero();36 tmpM->setZero(); 37 37 } 38 38 -
src/unittests/BoxUnittest.cpp
rcfc53b r9eb7580 41 41 ASSERT_DO(Assert::Throw); 42 42 unit = new Matrix; 43 unit-> one();43 unit->setIdentity(); 44 44 zero = new Matrix; 45 45 invertible = new Matrix; … … 52 52 unitBox = new Box; 53 53 stretchedBox1 = new Box; 54 boxMat. one();54 boxMat.setIdentity(); 55 55 boxMat.diagonal() = Vector(1,2,3); 56 56 stretchedBox1->setM(boxMat); 57 57 58 58 stretchedBox2 = new Box; 59 boxMat. one();59 boxMat.setIdentity(); 60 60 boxMat.diagonal() = Vector(2,3,1); 61 61 stretchedBox2->setM(boxMat); 62 62 63 63 stretchedBox3 = new Box; 64 boxMat. one();64 boxMat.setIdentity(); 65 65 boxMat.diagonal() = Vector(3,1,2); 66 66 stretchedBox3->setM(boxMat); 67 67 68 68 stretchedBox4 = new Box; 69 boxMat. one();69 boxMat.setIdentity(); 70 70 boxMat.diagonal() = Vector(2,2,2); 71 71 stretchedBox4->setM(boxMat); 72 72 73 73 tiltedBox1 = new Box; 74 boxMat. one();74 boxMat.setIdentity(); 75 75 boxMat.column(0) = Vector(1,0,1); 76 76 tiltedBox1->setM(boxMat); 77 77 78 78 tiltedBox2 = new Box; 79 boxMat. one();79 boxMat.setIdentity(); 80 80 boxMat.column(0) = Vector(1,1,1); 81 81 tiltedBox2->setM(boxMat); 82 82 83 83 tiltedBox3 = new Box; 84 boxMat. one();84 boxMat.setIdentity(); 85 85 boxMat.column(1) = Vector(0,1,1); 86 86 tiltedBox3->setM(boxMat); 87 87 88 88 tiltedBox4 = new Box; 89 boxMat. one();89 boxMat.setIdentity(); 90 90 boxMat.column(0) = Vector(1,1,1); 91 91 boxMat.column(1) = Vector(0,1,1); -
src/unittests/MatrixUnittest.cpp
rcfc53b r9eb7580 121 121 CPPUNIT_ASSERT_EQUAL(mat.diagonal(),zeroVec); 122 122 123 mat. one();123 mat.setIdentity(); 124 124 CPPUNIT_ASSERT_EQUAL(mat.row(0),unitVec[0]); 125 125 CPPUNIT_ASSERT_EQUAL(mat.row(1),unitVec[1]); … … 171 171 172 172 // transpose of unit is unit 173 res. one();173 res.setIdentity(); 174 174 (const Matrix)res.transpose(); 175 175 CPPUNIT_ASSERT_EQUAL(res,*one); 176 176 177 177 // transpose of transpose is same matrix 178 res. zero();178 res.setZero(); 179 179 res.set(2,2, 1.); 180 180 CPPUNIT_ASSERT_EQUAL(res.transpose().transpose(),res); … … 258 258 259 259 // zero rotation angles yields unity matrix 260 res. rotation(0,0,0);260 res.setRotation(0,0,0); 261 261 CPPUNIT_ASSERT_EQUAL(*one, res); 262 262 263 263 // arbitrary rotation matrix has det = 1 264 res. rotation(M_PI/3.,1.,M_PI/7.);264 res.setRotation(M_PI/3.,1.,M_PI/7.); 265 265 CPPUNIT_ASSERT(fabs(fabs(res.determinant()) -1.) < MYEPSILON); 266 266 267 267 // inverse is rotation matrix with negative angles 268 res. rotation(M_PI/3.,0.,0.);269 inverse. rotation(-M_PI/3.,0.,0.);268 res.setRotation(M_PI/3.,0.,0.); 269 inverse.setRotation(-M_PI/3.,0.,0.); 270 270 CPPUNIT_ASSERT_EQUAL(*one, res * inverse); 271 271 272 272 // ... or transposed 273 res. rotation(M_PI/3.,0.,0.);273 res.setRotation(M_PI/3.,0.,0.); 274 274 CPPUNIT_ASSERT_EQUAL(inverse, ((const Matrix) res).transpose()); 275 275 }
Note:
See TracChangeset
for help on using the changeset viewer.