kindred/solver
1
0
Fork 0
Code Issues 10 Pull Requests 1 Actions Packages Projects Releases Wiki Activity

works locally, but imagine it breaks elsewhere

Browse Source
This commit is contained in:
John Dupuy
2023-11-06 16:31:57 -06:00
parent 84ec3cf5e1
commit 87bd168b50
18 changed files with 83 additions and 59 deletions
Download Patch File Download Diff File Expand all files Collapse all files

47
OndselSolver/MBDynNode.cpp
View File

@@ -20,28 +20,29 @@ void MbD::MBDynNode::outputLine(int i, std::ostream& os)
auto x = asmtPart->xs->at(i);
auto y = asmtPart->ys->at(i);
auto z = asmtPart->zs->at(i);
auto bryantAngles = std::make_shared<EulerAngles<double>>();
bryantAngles->setRotOrder(1, 2, 3);
bryantAngles->at(0) = asmtPart->bryxs->at(i);
bryantAngles->at(1) = asmtPart->bryys->at(i);
bryantAngles->at(2) = asmtPart->bryzs->at(i);
bryantAngles->calc();
auto aA = bryantAngles->aA;
auto vx = asmtPart->vxs->at(i);
auto vy = asmtPart->vys->at(i);
auto vz = asmtPart->vzs->at(i);
auto omex = asmtPart->omexs->at(i);
auto omey = asmtPart->omeys->at(i);
auto omez = asmtPart->omezs->at(i);
os << id << " ";
os << x << " " << y << " " << z << " ";
auto row = aA->at(0);
os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
row = aA->at(1);
os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
row = aA->at(2);
os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
os << vx << " " << vy << " " << vz << " ";
os << omex << " " << omey << " " << omez << " ";
// TODO: undo the breaking I just did on purpose
// auto bryantAngles = std::make_shared<EulerAngles<double>>();
// bryantAngles->setRotOrder(1, 2, 3);
// bryantAngles->at(0) = asmtPart->bryxs->at(i);
// bryantAngles->at(1) = asmtPart->bryys->at(i);
// bryantAngles->at(2) = asmtPart->bryzs->at(i);
// bryantAngles->calc();
// auto aA = bryantAngles->aA;
// auto vx = asmtPart->vxs->at(i);
// auto vy = asmtPart->vys->at(i);
// auto vz = asmtPart->vzs->at(i);
// auto omex = asmtPart->omexs->at(i);
// auto omey = asmtPart->omeys->at(i);
// auto omez = asmtPart->omezs->at(i);
// os << id << " ";
// os << x << " " << y << " " << z << " ";
// auto row = aA->at(0);
// os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
// row = aA->at(1);
// os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
// row = aA->at(2);
// os << row->at(0) << " " << row->at(1) << " " << row->at(2) << " ";
// os << vx << " " << vy << " " << vz << " ";
// os << omex << " " << omey << " " << omez << " ";
os << std::endl;
}