58 lines
2.3 KiB
C++
58 lines
2.3 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2023 Ondsel, Inc. *
|
|
* *
|
|
* This file is part of OndselSolver. *
|
|
* *
|
|
* See LICENSE file for details about copyright. *
|
|
***************************************************************************/
|
|
|
|
#pragma once
|
|
|
|
#include "EigenDecomposition.h"
|
|
|
|
namespace MbD {
|
|
class MomentOfInertiaSolver : public EigenDecomposition
|
|
{
|
|
//See document 9100moment.fodt
|
|
//m aJPP aJoo rPoP aAPo aJcmP aJcmPcopy rPcmP aJpp aAPp colOrder
|
|
//aJoo == aJpp when rPoP == rPcmP and aAPo == aAPp
|
|
public:
|
|
static void example1();
|
|
void doFullPivoting(size_t p);
|
|
void forwardEliminateWithPivot(size_t p) override;
|
|
void backSubstituteIntoDU() override;
|
|
void postSolve() override;
|
|
FColDsptr basicSolvewithsaveOriginal(FMatDsptr fullMat, FColDsptr fullCol, bool saveOriginal) override;
|
|
FColDsptr basicSolvewithsaveOriginal(SpMatDsptr spMat, FColDsptr fullCol, bool saveOriginal) override;
|
|
void preSolvewithsaveOriginal(FMatDsptr fullMat, FColDsptr fullCol, bool saveOriginal) override;
|
|
void preSolvewithsaveOriginal(SpMatDsptr spMat, FColDsptr fullCol, bool saveOriginal) override;
|
|
double getmatrixArowimaxMagnitude(size_t i) override;
|
|
void doPivoting(size_t p) override;
|
|
|
|
void setm(double mass);
|
|
void setJPP(FMatDsptr mat);
|
|
void setrPoP(FColDsptr col);
|
|
void setAPo(FMatDsptr mat);
|
|
void setrPcmP(FColDsptr col);
|
|
FMatDsptr getJoo();
|
|
DiagMatDsptr getJpp();
|
|
FMatDsptr getAPp();
|
|
void calc();
|
|
void calcJoo();
|
|
void calcJpp();
|
|
void calcAPp();
|
|
FColDsptr eigenvectorFor(double lam);
|
|
void calcJppFromDiagJcmP();
|
|
void calcJppFromFullJcmP();
|
|
double modifiedArcCos(double val);
|
|
|
|
|
|
double m = 0.0;
|
|
FMatDsptr aJPP, aJoo, aAPo, aJcmP, aJcmPcopy, aAPp;
|
|
FColDsptr rPoP, rPcmP;
|
|
DiagMatDsptr aJpp;
|
|
std::shared_ptr<FullRow<size_t>> colOrder;
|
|
|
|
};
|
|
}
|