Files
solver/GNN/OndselSolver/BasicIntegrator.h
forbes-0023 98051ba0c9 feat: add Phase 1 constraint solver addon, move prior content to GNN/
- Move existing OndselSolver, GNN ML layer, and tooling into GNN/
  directory for integration in later phases
- Add Create addon scaffold: package.xml, Init.py
- Add expression DAG with eval, symbolic diff, simplification
- Add parameter table with fixed/free variable tracking
- Add quaternion rotation as polynomial Expr trees
- Add RigidBody entity (7 DOF: position + unit quaternion)
- Add constraint classes: Coincident, DistancePointPoint, Fixed
- Add Newton-Raphson solver with symbolic Jacobian + numpy lstsq
- Add pre-solve passes: substitution + single-equation
- Add DOF counting via Jacobian SVD rank
- Add KindredSolver IKCSolver bridge for kcsolve integration
- Add 82 unit tests covering all modules

Registers as 'kindred' solver via kcsolve.register_solver() when
loaded by Create's addon_loader.
2026-02-20 20:35:47 -06:00

59 lines
2.0 KiB
C++

/***************************************************************************
* Copyright (c) 2023 Ondsel, Inc. *
* *
* This file is part of OndselSolver. *
* *
* See LICENSE file for details about copyright. *
***************************************************************************/
#pragma once
#include <vector>
#include "Integrator.h"
namespace MbD {
class IntegratorInterface;
class DifferenceOperator;
class BasicIntegrator : public Integrator
{
//istep iTry maxTry tpast t tnew h hnew order orderNew orderMax opBDF continue
public:
virtual void calcOperatorMatrix();
virtual void incrementTime();
virtual void incrementTry();
void initialize() override;
void initializeGlobally() override;
void initializeLocally() override;
void iStep(size_t i) override;
void postFirstStep() override;
void postStep() override;
void postRun() override;
void preFirstStep() override;
void preRun() override;
void preStep() override;
void reportStats() override;
void run() override;
void selectOrder() override;
void subsequentSteps() override;
void setSystem(Solver* sys) override;
void logString(const std::string& str) override;
virtual void setorder(size_t o);
virtual void settnew(double t);
virtual void sett(double t);
void settime(double t);
double tprevious();
IntegratorInterface* system;
size_t istep = 0, iTry = 0, maxTry = 0;
std::shared_ptr<std::vector<double>> tpast;
double t = 0.0, tnew = 0.0, h = 0, hnew = 0.0;
size_t order = 0, orderNew = 0, orderMax = 0;
std::shared_ptr<DifferenceOperator> opBDF;
bool _continue = false;
};
}