//----------------------------------------------------------------------------- // All declarations not grouped specially elsewhere. // // Copyright 2008-2013 Jonathan Westhues. //----------------------------------------------------------------------------- #ifndef SOLVESPACE_H #define SOLVESPACE_H #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include // We declare these in advance instead of simply using FT_Library // (defined as typedef FT_LibraryRec_* FT_Library) because including // freetype.h invokes indescribable horrors and we would like to avoid // doing that every time we include solvespace.h. struct FT_LibraryRec_; struct FT_FaceRec_; typedef struct _cairo cairo_t; typedef struct _cairo_surface cairo_surface_t; // The few floating-point equality comparisons in SolveSpace have been // carefully considered, so we disable the -Wfloat-equal warning for them #ifdef __clang__ # define EXACT(expr) \ (_Pragma("clang diagnostic push") \ _Pragma("clang diagnostic ignored \"-Wfloat-equal\"") \ (expr) \ _Pragma("clang diagnostic pop")) #else # define EXACT(expr) (expr) #endif // Debugging functions #if defined(__GNUC__) #define ssassert(condition, message) \ do { \ if(__builtin_expect((condition), true) == false) { \ SolveSpace::AssertFailure(__FILE__, __LINE__, __func__, #condition, message); \ __builtin_unreachable(); \ } \ } while(0) #else #define ssassert(condition, message) \ do { \ if((condition) == false) { \ SolveSpace::AssertFailure(__FILE__, __LINE__, __func__, #condition, message); \ abort(); \ } \ } while(0) #endif #define dbp SolveSpace::Platform::DebugPrint #define DBPTRI(tri) \ dbp("tri: (%.3f %.3f %.3f) (%.3f %.3f %.3f) (%.3f %.3f %.3f)", \ CO((tri).a), CO((tri).b), CO((tri).c)) namespace SolveSpace { using std::min; using std::max; using std::swap; using std::fabs; [[noreturn]] void AssertFailure(const char *file, unsigned line, const char *function, const char *condition, const char *message); #if defined(__GNUC__) __attribute__((__format__ (__printf__, 1, 2))) #endif std::string ssprintf(const char *fmt, ...); inline bool IsReasonable(double x) { return std::isnan(x) || x > 1e11 || x < -1e11; } inline int WRAP(int v, int n) { // Clamp it to the range [0, n) while(v >= n) v -= n; while(v < 0) v += n; return v; } inline double WRAP_NOT_0(double v, double n) { // Clamp it to the range (0, n] while(v > n) v -= n; while(v <= 0) v += n; return v; } inline double WRAP_SYMMETRIC(double v, double n) { // Clamp it to the range (-n/2, n/2] while(v > n/2) v -= n; while(v <= -n/2) v += n; return v; } #define CO(v) (v).x, (v).y, (v).z static constexpr double ANGLE_COS_EPS = 1e-6; static constexpr double LENGTH_EPS = 1e-6; static constexpr double VERY_POSITIVE = 1e10; static constexpr double VERY_NEGATIVE = -1e10; inline double Random(double vmax) { return (vmax*rand()) / RAND_MAX; } #include "platform/platform.h" #include "platform/gui.h" #include "resource.h" using Platform::AllocTemporary; using Platform::FreeAllTemporary; class Expr; class ExprVector; class ExprQuaternion; class RgbaColor; enum class Command : uint32_t; enum class Unit : uint32_t { MM = 0, INCHES, METERS }; template using handle_map = std::map; class Group; class SSurface; #include "dsc.h" #include "polygon.h" #include "srf/surface.h" #include "render/render.h" class Entity; class hEntity; class Param; class hParam; typedef IdList EntityList; typedef IdList ParamList; enum class SolveResult : uint32_t { OKAY = 0, DIDNT_CONVERGE = 10, REDUNDANT_OKAY = 11, REDUNDANT_DIDNT_CONVERGE = 12, TOO_MANY_UNKNOWNS = 20 }; #include "sketch.h" #include "ui.h" #include "expr.h" // Utility functions that are provided in the platform-independent code. class utf8_iterator : std::iterator { const char *p, *n; public: utf8_iterator(const char *p) : p(p), n(NULL) {} bool operator==(const utf8_iterator &i) const { return p==i.p; } bool operator!=(const utf8_iterator &i) const { return p!=i.p; } ptrdiff_t operator- (const utf8_iterator &i) const { return p -i.p; } utf8_iterator& operator++() { **this; p=n; n=NULL; return *this; } utf8_iterator operator++(int) { utf8_iterator t(*this); operator++(); return t; } char32_t operator*(); const char* ptr() const { return p; } }; class ReadUTF8 { const std::string &str; public: ReadUTF8(const std::string &str) : str(str) {} utf8_iterator begin() const { return utf8_iterator(&str[0]); } utf8_iterator end() const { return utf8_iterator(&str[0] + str.length()); } }; #define arraylen(x) (sizeof((x))/sizeof((x)[0])) #define PI (3.1415926535897931) void MakeMatrix(double *mat, double a11, double a12, double a13, double a14, double a21, double a22, double a23, double a24, double a31, double a32, double a33, double a34, double a41, double a42, double a43, double a44); void MultMatrix(double *mata, double *matb, double *matr); int64_t GetMilliseconds(); void Message(const char *fmt, ...); void MessageAndRun(std::function onDismiss, const char *fmt, ...); void Error(const char *fmt, ...); class System { public: enum { MAX_UNKNOWNS = 1024 }; EntityList entity; ParamList param; IdList eq; // A list of parameters that are being dragged; these are the ones that // we should put as close as possible to their initial positions. List dragged; enum { // In general, the tag indicates the subsys that a variable/equation // has been assigned to; these are exceptions for variables: VAR_SUBSTITUTED = 10000, VAR_DOF_TEST = 10001, // and for equations: EQ_SUBSTITUTED = 20000 }; // The system Jacobian matrix struct { // The corresponding equation for each row hEquation eq[MAX_UNKNOWNS]; // The corresponding parameter for each column hParam param[MAX_UNKNOWNS]; // We're solving AX = B int m, n; struct { Expr *sym[MAX_UNKNOWNS][MAX_UNKNOWNS]; double num[MAX_UNKNOWNS][MAX_UNKNOWNS]; } A; double scale[MAX_UNKNOWNS]; // Some helpers for the least squares solve double AAt[MAX_UNKNOWNS][MAX_UNKNOWNS]; double Z[MAX_UNKNOWNS]; double X[MAX_UNKNOWNS]; struct { Expr *sym[MAX_UNKNOWNS]; double num[MAX_UNKNOWNS]; } B; } mat; static const double RANK_MAG_TOLERANCE, CONVERGE_TOLERANCE; int CalculateRank(); bool TestRank(int *rank = NULL); static bool SolveLinearSystem(double X[], double A[][MAX_UNKNOWNS], double B[], int N); bool SolveLeastSquares(); bool WriteJacobian(int tag); void EvalJacobian(); void WriteEquationsExceptFor(hConstraint hc, Group *g); void FindWhichToRemoveToFixJacobian(Group *g, List *bad, bool forceDofCheck); void SolveBySubstitution(); bool IsDragged(hParam p); bool NewtonSolve(int tag); void MarkParamsFree(bool findFree); int CalculateDof(); SolveResult Solve(Group *g, int *rank = NULL, int *dof = NULL, List *bad = NULL, bool andFindBad = false, bool andFindFree = false, bool forceDofCheck = false); SolveResult SolveRank(Group *g, int *rank = NULL, int *dof = NULL, List *bad = NULL, bool andFindBad = false, bool andFindFree = false); void Clear(); }; #include "ttf.h" class StepFileWriter { public: void ExportSurfacesTo(const Platform::Path &filename); void WriteHeader(); void WriteProductHeader(); int ExportCurve(SBezier *sb); int ExportCurveLoop(SBezierLoop *loop, bool inner); void ExportSurface(SSurface *ss, SBezierList *sbl); void WriteWireframe(); void WriteFooter(); List curves; List advancedFaces; FILE *f; int id; }; class VectorFileWriter { protected: Vector u, v, n, origin; double cameraTan, scale; public: FILE *f; Platform::Path filename; Vector ptMin, ptMax; static double MmToPts(double mm); static VectorFileWriter *ForFile(const Platform::Path &filename); void SetModelviewProjection(const Vector &u, const Vector &v, const Vector &n, const Vector &origin, double cameraTan, double scale); Vector Transform(Vector &pos) const; void OutputLinesAndMesh(SBezierLoopSetSet *sblss, SMesh *sm); void BezierAsPwl(SBezier *sb); void BezierAsNonrationalCubic(SBezier *sb, int depth=0); virtual void StartPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) = 0; virtual void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) = 0; virtual void Bezier(SBezier *sb) = 0; virtual void Triangle(STriangle *tr) = 0; virtual bool OutputConstraints(IdList *) { return false; } virtual void StartFile() = 0; virtual void FinishAndCloseFile() = 0; virtual bool HasCanvasSize() const = 0; virtual bool CanOutputMesh() const = 0; }; class DxfFileWriter : public VectorFileWriter { public: struct BezierPath { std::vector beziers; }; std::vector paths; IdList *constraint; static const char *lineTypeName(StipplePattern stippleType); bool OutputConstraints(IdList *constraint) override; void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return false; } bool CanOutputMesh() const override { return false; } bool NeedToOutput(Constraint *c); }; class EpsFileWriter : public VectorFileWriter { public: Vector prevPt; void MaybeMoveTo(Vector s, Vector f); void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return true; } bool CanOutputMesh() const override { return true; } }; class PdfFileWriter : public VectorFileWriter { public: uint32_t xref[10]; uint32_t bodyStart; Vector prevPt; void MaybeMoveTo(Vector s, Vector f); void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return true; } bool CanOutputMesh() const override { return true; } }; class SvgFileWriter : public VectorFileWriter { public: Vector prevPt; void MaybeMoveTo(Vector s, Vector f); void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return true; } bool CanOutputMesh() const override { return true; } }; class HpglFileWriter : public VectorFileWriter { public: static double MmToHpglUnits(double mm); void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return false; } bool CanOutputMesh() const override { return false; } }; class Step2dFileWriter : public VectorFileWriter { StepFileWriter sfw; void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return false; } bool CanOutputMesh() const override { return false; } }; class GCodeFileWriter : public VectorFileWriter { public: SEdgeList sel; void StartPath( RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void FinishPath(RgbaColor strokeRgb, double lineWidth, bool filled, RgbaColor fillRgb, hStyle hs) override; void Triangle(STriangle *tr) override; void Bezier(SBezier *sb) override; void StartFile() override; void FinishAndCloseFile() override; bool HasCanvasSize() const override { return false; } bool CanOutputMesh() const override { return false; } }; #ifdef LIBRARY # define ENTITY EntityBase # define CONSTRAINT ConstraintBase #else # define ENTITY Entity # define CONSTRAINT Constraint #endif class Sketch { public: // These are user-editable, and define the sketch. IdList group; List
groupOrder; IdList constraint; IdList request; IdList style; // These are generated from the above. IdList entity; IdList param; inline CONSTRAINT *GetConstraint(hConstraint h) { return constraint.FindById(h); } inline ENTITY *GetEntity (hEntity h) { return entity. FindById(h); } inline Param *GetParam (hParam h) { return param. FindById(h); } inline Request *GetRequest(hRequest h) { return request.FindById(h); } inline Group *GetGroup (hGroup h) { return group. FindById(h); } // Styles are handled a bit differently. void Clear(); BBox CalculateEntityBBox(bool includingInvisible); Group *GetRunningMeshGroupFor(hGroup h); }; #undef ENTITY #undef CONSTRAINT class SolveSpaceUI { public: TextWindow *pTW; TextWindow &TW; GraphicsWindow GW; // The state for undo/redo typedef struct { IdList group; List
groupOrder; IdList request; IdList constraint; IdList param; IdList style; hGroup activeGroup; void Clear() { group.Clear(); request.Clear(); constraint.Clear(); param.Clear(); style.Clear(); } } UndoState; enum { MAX_UNDO = 16 }; typedef struct { UndoState d[MAX_UNDO]; int cnt; int write; } UndoStack; UndoStack undo; UndoStack redo; std::map, Platform::PathLess> images; bool ReloadLinkedImage(const Platform::Path &saveFile, Platform::Path *filename, bool canCancel); void UndoEnableMenus(); void UndoRemember(); void UndoUndo(); void UndoRedo(); void PushFromCurrentOnto(UndoStack *uk); void PopOntoCurrentFrom(UndoStack *uk); void UndoClearState(UndoState *ut); void UndoClearStack(UndoStack *uk); // Little bits of extra configuration state enum { MODEL_COLORS = 8 }; RgbaColor modelColor[MODEL_COLORS]; Vector lightDir[2]; double lightIntensity[2]; double ambientIntensity; double chordTol; double chordTolCalculated; int maxSegments; double exportChordTol; int exportMaxSegments; double cameraTangent; double gridSpacing; double exportScale; double exportOffset; bool fixExportColors; bool drawBackFaces; bool showContourAreas; bool checkClosedContour; bool turntableNav; bool immediatelyEditDimension; bool automaticLineConstraints; bool showToolbar; Platform::Path screenshotFile; RgbaColor backgroundColor; bool exportShadedTriangles; bool exportPwlCurves; bool exportCanvasSizeAuto; bool exportMode; struct { double left; double right; double bottom; double top; } exportMargin; struct { double width; double height; double dx; double dy; } exportCanvas; struct { double depth; int passes; double feed; double plungeFeed; } gCode; Unit viewUnits; int afterDecimalMm; int afterDecimalInch; int afterDecimalDegree; bool useSIPrefixes; int autosaveInterval; // in minutes std::string MmToString(double v); std::string MmToStringSI(double v, int dim = 0); std::string DegreeToString(double v); double ExprToMm(Expr *e); double StringToMm(const std::string &s); const char *UnitName(); double MmPerUnit(); int UnitDigitsAfterDecimal(); void SetUnitDigitsAfterDecimal(int v); double ChordTolMm(); double ExportChordTolMm(); int GetMaxSegments(); bool usePerspectiveProj; double CameraTangent(); // Some stuff relating to the tangent arcs created non-parametrically // as special requests. double tangentArcRadius; bool tangentArcManual; bool tangentArcModify; // The platform-dependent code calls this before entering the msg loop void Init(); void Exit(); // File load/save routines, including the additional files that get // loaded when we have link groups. FILE *fh; void AfterNewFile(); void AddToRecentList(const Platform::Path &filename); Platform::Path saveFile; bool fileLoadError; bool unsaved; typedef struct { char type; const char *desc; char fmt; void *ptr; } SaveTable; static const SaveTable SAVED[]; void SaveUsingTable(const Platform::Path &filename, int type); void LoadUsingTable(const Platform::Path &filename, char *key, char *val); struct { Group g; Request r; Entity e; Param p; Constraint c; Style s; } sv; static void MenuFile(Command id); void Autosave(); void RemoveAutosave(); static constexpr size_t MAX_RECENT = 8; static constexpr const char *SKETCH_EXT = "slvs"; static constexpr const char *BACKUP_EXT = "slvs~"; std::vector recentFiles; bool Load(const Platform::Path &filename); bool GetFilenameAndSave(bool saveAs); bool OkayToStartNewFile(); hGroup CreateDefaultDrawingGroup(); void UpdateWindowTitles(); void ClearExisting(); void NewFile(); bool SaveToFile(const Platform::Path &filename); bool LoadAutosaveFor(const Platform::Path &filename); bool LoadFromFile(const Platform::Path &filename, bool canCancel = false); void UpgradeLegacyData(); bool LoadEntitiesFromFile(const Platform::Path &filename, EntityList *le, SMesh *m, SShell *sh); bool ReloadAllLinked(const Platform::Path &filename, bool canCancel = false); // And the various export options void ExportAsPngTo(const Platform::Path &filename); void ExportMeshTo(const Platform::Path &filename); void ExportMeshAsStlTo(FILE *f, SMesh *sm); void ExportMeshAsQ3doTo(FILE *f, SMesh *sm); void ExportMeshAsObjTo(FILE *fObj, FILE *fMtl, SMesh *sm); void ExportMeshAsThreeJsTo(FILE *f, const Platform::Path &filename, SMesh *sm, SOutlineList *sol); void ExportMeshAsVrmlTo(FILE *f, const Platform::Path &filename, SMesh *sm); void ExportViewOrWireframeTo(const Platform::Path &filename, bool exportWireframe); void ExportSectionTo(const Platform::Path &filename); void ExportWireframeCurves(SEdgeList *sel, SBezierList *sbl, VectorFileWriter *out); void ExportLinesAndMesh(SEdgeList *sel, SBezierList *sbl, SMesh *sm, Vector u, Vector v, Vector n, Vector origin, double cameraTan, VectorFileWriter *out); static void MenuAnalyze(Command id); // Additional display stuff struct { SContour path; hEntity point; } traced; SEdgeList nakedEdges; struct { bool draw; Vector ptA; Vector ptB; } extraLine; struct { bool draw, showOrigin; Vector pt, u, v; } justExportedInfo; struct { bool draw; bool dirty; Vector position; } centerOfMass; class Clipboard { public: List r; List c; void Clear(); bool ContainsEntity(hEntity old); hEntity NewEntityFor(hEntity old); }; Clipboard clipboard; void MarkGroupDirty(hGroup hg, bool onlyThis = false); void MarkGroupDirtyByEntity(hEntity he); // Consistency checking on the sketch: stuff with missing dependencies // will get deleted automatically. struct { int requests; int groups; int constraints; int nonTrivialConstraints; } deleted; bool GroupExists(hGroup hg); bool PruneOrphans(); bool EntityExists(hEntity he); bool GroupsInOrder(hGroup before, hGroup after); bool PruneGroups(hGroup hg); bool PruneRequests(hGroup hg); bool PruneConstraints(hGroup hg); static void ShowNakedEdges(bool reportOnlyWhenNotOkay); enum class Generate : uint32_t { DIRTY, ALL, REGEN, UNTIL_ACTIVE, }; void GenerateAll(Generate type = Generate::DIRTY, bool andFindFree = false, bool genForBBox = false); void SolveGroup(hGroup hg, bool andFindFree); void SolveGroupAndReport(hGroup hg, bool andFindFree); SolveResult TestRankForGroup(hGroup hg, int *rank = NULL); void WriteEqSystemForGroup(hGroup hg); void MarkDraggedParams(); void ForceReferences(); void UpdateCenterOfMass(); bool ActiveGroupsOkay(); // The system to be solved. System *pSys; System &sys; // All the TrueType fonts in memory TtfFontList fonts; // Everything has been pruned, so we know there's no dangling references // to entities that don't exist. Before that, we mustn't try to display // the sketch! bool allConsistent; Platform::TimerRef showTWTimer; Platform::TimerRef generateAllTimer; Platform::TimerRef autosaveTimer; void ScheduleShowTW(); void ScheduleGenerateAll(); void ScheduleAutosave(); static void MenuHelp(Command id); void Clear(); // We allocate TW and sys on the heap to work around an MSVC problem // where it puts zero-initialized global data in the binary (~30M of zeroes) // in release builds. SolveSpaceUI() : pTW(new TextWindow()), TW(*pTW), pSys(new System()), sys(*pSys) {} ~SolveSpaceUI() { delete pTW; delete pSys; } }; void ImportDxf(const Platform::Path &file); void ImportDwg(const Platform::Path &file); extern SolveSpaceUI SS; extern Sketch SK; } #ifndef __OBJC__ using namespace SolveSpace; #endif #endif