zlacker

[parent] [thread] 32 comments
1. digdug+(OP)[view] [source] 2023-10-20 05:07:12
The world doesn't need Text-to-CAD. The world needs a fully capable open source parametric 3D geometric CAD kernel.

Solidworks, Creo, AutoCAD, Fusion, etc., can all take their bug ridden unoptimized single threaded rent-seeking monstrosities and stick em where the sun don't shine.

Seriously - if anyone wants to create an absolutely world-changing piece of software, start working on a new CAD kernel that takes the last 50 years of computer science advances into account, because none of the entrenched industry standards have done so. Don't worry about having to provide customer service, because none of the entrenched industry standards worry about that either.

And no - while openCascade and solvespace are impressive, they aren't fully capable, nor do they start from a modern foundation.

replies(8): >>auggie+P >>dymk+G1 >>pca006+ha >>jf___+Vc >>wizzle+af >>virtua+8z >>mhb+4N >>jessfr+Mx1
2. auggie+P[view] [source] 2023-10-20 05:23:33
>>digdug+(OP)
So what's so difficult about such a kernel? What does a modern foundation look like?
replies(1): >>dymk+q1
◧◩
3. dymk+q1[view] [source] [discussion] 2023-10-20 05:33:24
>>auggie+P
The math is quite difficult to do right, and there's a billion corner cases to make a kernel useful for real world designs. Take a fillet: It needs to handle inside corners, outside corners, compound angles coming in from arbitrary numbers of directions, it probably needs the ability to vary along its distance, create more geometry when adjacent faces don't leave enough room, etc etc.

That's just the start of a single feature type. Now you need a bunch more feature types, and they all need to interact well with each other. The kernel also needs some way of solving the topological naming problem to be useful (FreeCAD might get a basic version of this after a decade(?) of work).

It's probably tantamount to writing a modern-day browser in terms of complexity.

replies(1): >>auggie+b4
4. dymk+G1[view] [source] 2023-10-20 05:35:51
>>digdug+(OP)
https://www.fornjot.app/ is an attempt, and also good at showing what a large undertaking such a goal is
◧◩◪
5. auggie+b4[view] [source] [discussion] 2023-10-20 06:05:13
>>dymk+q1
I've written some custom code for computational geometry (like computing the offset for a geometric object using Apple Metal). It was a lot of fun, but also quite hard, and a lot of edge cases I just didn't deal with because I had a particular use case and speed was paramount.

Maybe the idea of a "kernel" is the problem here. A kernel the size of a browser is not a kernel.

I think what's really needed is a full-blown integration with a theorem proving system (which has an easier to define kernel of its own).

replies(1): >>pca006+qd
6. pca006+ha[view] [source] 2023-10-20 07:07:43
>>digdug+(OP)
I think the interesting thing about CAD kernel is that there are different representations and limitations to each representation. You have triangular (or polygonal) mesh, BREP which uses NURBS, SDF which is based on functional representation. I have experience working with triangular meshes and SDF so here are my opinions about them, please correct me if I am wrong:

Triangular mesh is conceptually simple, but requires many faces to approximate curved surfaces with high precision (you may be able to use subdivision surface in some cases but intersection/union in those cases are more challenging). Also, for more complicated models, floating point errors really add up and you either have to use an exact representation (which is really slow) or try some other approaches which can be robust w.r.t. errors (e.g. https://github.com/elalish/manifold but it is really hard to get right). Another disadvantage comparing with BREP is the lack of constraint solving, which I will write about it below.

SDF is nice for mathematically defined objects. They are computationally intensive, so some sdf libraries use GPU to speedup the computation. There are approaches that can speed up the evaluation, but doesn't work well if the function is not really the distance (https://github.com/curv3d/curv/blob/master/docs/shapes/Shape...).

-----

Constraints solving: This is a big problem with mesh-based CAD. Traditional CAD usually allows you to have under-defined constraints, and users can iteratively set constraints until the model is fully defined. There is no such a thing (yet) with mesh-based CAD. Also, we don't really have nice ways to represent constraints relative to curved surfaces because there is no curved surface in our mesh...

Also, one particular challenge with text-based (or code-based) CAD is how to select the surfaces with an ergonomic API. GUI can solve this problem but writing a good GUI is a complicated task (that I am not willing to touch).

7. jf___+Vc[view] [source] 2023-10-20 07:42:18
>>digdug+(OP)
It seems a bizarre statement to state that OpenCASCADE isn't fully capable. Its the only OS licensed kernel that'll read a STEP file. Also "modern foundation" is a misleading statement, any CAD kernel bearing any kind of relevance seems implying a codebase that's been around for a quarter century. Like it or not OpenCASCADE is the hand that was dealt. I've worked with the technology [1] extensively and it provided the underpinnings for a startup I've founded [2]. pythonocc is the bees knees, it allows you to develop a proper CAD app. I'd consider CGAL a modern kernel, but it doesn't cover CAD since there is no BRep support [4] Don't take my word for it, but see also the many publications that have built on the tech [3]

[1] https://github.com/tpaviot/pythonocc-core [2] odico.dk/ [3] https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=pyth... [4] https://www.cgal.org/

replies(4): >>pca006+Ie >>wizzle+lf >>RobotT+9D >>digdug+GN
◧◩◪◨
8. pca006+qd[view] [source] [discussion] 2023-10-20 07:48:19
>>auggie+b4
I doubt if it is possible to have a practical small kernel. Take mesh processing in https://github.com/elalish/manifold for example, we encountered a lot of problems when trying to deal with inexact floating point arithmetic. Using exact arithmetic can probably result in much simpler code (doesn't make sense in our case because the point of the library is about dealing with inexact arithmetic), but exact means slow. Also, a lot of code is about fast-paths that occurs frequently, data structures and algorithms that cut down complexity but are difficult to implement. Can we remove those fast-paths and complex algorithms? Probably yes, but this will slow things down a lot.
replies(1): >>auggie+2r
◧◩
9. pca006+Ie[view] [source] [discussion] 2023-10-20 08:04:32
>>jf___+Vc
Wondering why no BRep support means it is not a CAD. It seems that OpenVSP doesn't use BRep but uses parametric surface [1]. I wonder if mesh-based modeling + some constraints solving will get you a CAD, or are there other requirements that I don't know? I only work on mesh processing library on my free-time and I don't know much about BRep.

[1] https://github.com/elalish/manifold/discussions/549#discussi...

10. wizzle+af[view] [source] 2023-10-20 08:08:57
>>digdug+(OP)
Have you worked with Open Cascade recently? As someone who works with it every day for developing a CAD application, it would be interesting to know what people see as its limitations. It's the only geometry kernel I've had access to, and it seems like an absolute gift
replies(2): >>Kosiri+Ik >>westur+sl
◧◩
11. wizzle+lf[view] [source] [discussion] 2023-10-20 08:11:38
>>jf___+Vc
It does seem a bit of a throwaway statement regarding OCCT - I also work with it every day and, for the most part, it has all the same eccentricities and limitations of any large heritage-listed C++ library. There's a lot it can do!
◧◩
12. Kosiri+Ik[view] [source] [discussion] 2023-10-20 09:04:38
>>wizzle+af
Can you share what you do with it? What kind of projects or industry is your application targeting?
◧◩
13. westur+sl[view] [source] [discussion] 2023-10-20 09:11:56
>>wizzle+af
cadquery wraps OCC (OpenCascades) but used to wrap freeCAD.

Here's a LEGO interlocking block brick in cadquery: https://cadquery.readthedocs.io/en/latest/examples.html#lego... .

awesome-cadquery: https://github.com/CadQuery/awesome-cadquery

cadquery and thus also jupyter-cadquery now have support for build123d.

gumyr/build123d https://github.com/gumyr/build123d :

> Build123d is a python-based, parametric, boundary representation (BREP) modeling framework for 2D and 3D CAD. It's built on the Open Cascade geometric kernel and allows for the creation of complex models using a simple and intuitive python syntax. Build123d can be used to create models for 3D printing, CNC machining, laser cutting, and other manufacturing processes. Models can be exported to a wide variety of popular CAD tools such as FreeCAD and SolidWorks.

> Build123d could be considered as an evolution of CadQuery where the somewhat restrictive Fluent API (method chaining) is replaced with stateful context managers* - e.g. with blocks - thus enabling the full python toolbox: for loops, references to objects, object sorting and filtering, etc.*

"Build123d: A Python CAD programming library" (2023) >>37576296

build123d docs > Tips & Best Practices: https://build123d.readthedocs.io/en/latest/tips.html

BREP: Boundary representation: https://en.wikipedia.org/wiki/Boundary_representation

Manim, Blender, ipyblender, PhysX, o3de, [FEM, CFD, [thermal, fluidic,] engineering]: https://github.com/ManimCommunity/manim/issues/3362

NURBS: Non-Uniform Rational B-Splines: https://en.wikipedia.org/wiki/Non-uniform_rational_B-spline

NURBS for COMPAS: test_curve.py, test_surface.py: https://github.com/gramaziokohler/compas_nurbs :

> This package is inspired by the NURBS-Python package, however uses a NumPy-based backend for better performance.

> Curve, and Surface are non-uniform non-rational B-Spline geometries (NUBS), RationalCurve, and RationalSurface are non-uniform rational B-Spline Geometries (NURBS). They all built upon the class BSpline. Coordinates have to be in 3D space (x, y, z)

https://github.com/compas-dev

compas_rhino, compas_blender,

- [ ] compas_o3de

Blender docs > Modeling Surfaces; NURBs implementation, limits, challenges: https://docs.blender.org/manual/en/latest/modeling/surfaces/...

/? "NURBS" opencascade https://www.google.com/search?q=%22nurbs%22+%22opencascade%2...

OCCT (OCC) Open Cascade Technology: https://en.wikipedia.org/wiki/Open_Cascade_Technology

OCCT > Standard Transient _ MMtg_TShared > Geom_Geometry > Geom_Curve > Geom_BoundedCurve > Geom_BSplineCurve https://dev.opencascade.org/doc/occt-6.9.1/refman/html/class...

OCC > Standard Transient _ MMtg_TShared > Geom_Geometry > Geom_Surface > Geom_BoundedSurface > Geom_BSplineSurface: https://dev.opencascade.org/doc/occt-6.9.1/refman/html/class...

Cadquery.Shape.toSplines(degree: int = 3, tolerance: float = 0.001, nurbs: bool = False)→ T https://cadquery.readthedocs.io/en/latest/classreference.htm...

◧◩◪◨⬒
14. auggie+2r[view] [source] [discussion] 2023-10-20 10:07:15
>>pca006+qd
I completely agree. You want to be able to be fast and correct, where correct also depends on some assumptions you can make depending on what you apply your algorithm to. That's why I think a theorem prover kernel is needed, which connects the various algorithms, and the theory behind them that is needed to prove correctness. Until you get there, you will have various different algorithms that don't really fit together, and where it is not clear how to make them work together and still be confident about the result.
replies(1): >>pca006+lI
15. virtua+8z[view] [source] 2023-10-20 11:26:33
>>digdug+(OP)
I agree 100%.

Truck[1] and Fornjot[2] are recent attempts in the Rust space, both are WIP.

But both seem to be going the traditional way. I.e. B-Rep that can be converted to (trimmed) NURBS.

I think if one wanted to incorporate the last 50 years of computer science, particularly computer graphics, one needed to broaden the feature set considerably.

You need support for precision subdivision surface modeling with variable radius creases (either via reverse subdivision where you make sure the limit surface pass through given constraints or using an interpolating subivision scheme that but has the same perks as e.g. Catmull-Clark).

Then you need to have SDF modeling ofc.

Possibly point based representations. If only as inputs.

And traditional B-Rep.

Finally, the kernel should be able to go back and forth lossless between these representations wherever possible.

And everything must be node-based, like e.g. Houdini. Completely non-destructive.

[1] https://github.com/ricosjp/truck

[2] https://github.com/hannobraun/fornjot

replies(2): >>ameliu+xA >>skippa+oQ
◧◩
16. ameliu+xA[view] [source] [discussion] 2023-10-20 11:38:17
>>virtua+8z
The main problem I have with most CAD software is that I can't do:

    1. create an object, e.g. a cube
    2. perform lots of random transformations
    3. perform the inverse of the above transformations
    4. subtract the object from the original object
    5. end up with exactly nothing
How are modern approaches solving this (robustness), if at all?
replies(3): >>ahmetr+fE >>abrouw+nJ >>iancmc+wS2
◧◩
17. RobotT+9D[view] [source] [discussion] 2023-10-20 12:01:16
>>jf___+Vc
>Its the only OS licensed kernel that'll read a STEP file.

BRL-CAD does that too.

◧◩◪
18. ahmetr+fE[view] [source] [discussion] 2023-10-20 12:12:17
>>ameliu+xA
What is the inverse operation of, say, beveling?
replies(2): >>ameliu+Pg1 >>iancmc+F73
◧◩◪◨⬒⬓
19. pca006+lI[view] [source] [discussion] 2023-10-20 12:45:24
>>auggie+2r
While this would be nice, I don't think enough people working on these algorithms are trained to use formal theorem proving tools.
replies(1): >>auggie+0N
◧◩◪
20. abrouw+nJ[view] [source] [discussion] 2023-10-20 12:52:15
>>ameliu+xA
As someone who does CAD work (product design) every day, I'm curious about the use case? Starting with Surfaces (in Creo, at least) can be provide such flexibility.
replies(1): >>ameliu+8Y
◧◩◪◨⬒⬓⬔
21. auggie+0N[view] [source] [discussion] 2023-10-20 13:20:15
>>pca006+lI
Yes, you are right. Time for better theorem proving tools! A theorem proving tool shouldn't add to your time spent designing and implementing, it should help you to get your stuff done faster and with (much) higher quality.
replies(1): >>namibj+es4
22. mhb+4N[view] [source] 2023-10-20 13:20:30
>>digdug+(OP)
Maybe there's some way to throw money at Rhino to get them into this space. Love Rhino and McNeel.
◧◩
23. digdug+GN[view] [source] [discussion] 2023-10-20 13:24:06
>>jf___+Vc
To be clear - I do think OpenCascade is impressive. Incredibly so, once one becomes aware of the magnitude of the problem it is trying to solve. I will also admit I haven't used it in the past couple of years, but when I did it's limitations in filleting and chamfering alone were enough to make it a non-starter for industry use.

My broader point was that there is a need to start from a new paradigm that leverages the possibilities of modern, highly parallel computing hardware. The hardware requirements for performant and reliable CAD software are incredibly high, and their reliance on high clock speed single core processors is quickly being left behind by modern processing hardware.

◧◩
24. skippa+oQ[view] [source] [discussion] 2023-10-20 13:38:34
>>virtua+8z
Where can one read about these? Are there modern textbooks/overview papers for modern methods?
◧◩◪◨
25. ameliu+8Y[view] [source] [discussion] 2023-10-20 14:11:58
>>abrouw+nJ
I don't have a practical usecase, it's just the fact that most CAD software is not mathematically correct that bothers me.
replies(1): >>dymk+Za4
◧◩◪◨
26. ameliu+Pg1[view] [source] [discussion] 2023-10-20 15:26:29
>>ahmetr+fE
Beveling is a lossy operation, so there is no inverse in general.
27. jessfr+Mx1[view] [source] 2023-10-20 16:48:19
>>digdug+(OP)
Our kernel isnt open source but we have an API, and we rebuilt our own CAD engine to fully utilize GPUs

https://kittycad.io/modeling-app https://kittycad.io/geometry-engine

While we are not fully capable yet, that's the end goal by next year. Years from now we cannot still be using these single-threaded nightmares from 30 years ago.

replies(1): >>digdug+Vpd
◧◩◪
28. iancmc+wS2[view] [source] [discussion] 2023-10-21 02:41:57
>>ameliu+xA
Solidworks would do this no problem
replies(1): >>iancmc+D73
◧◩◪◨
29. iancmc+D73[view] [source] [discussion] 2023-10-21 06:26:07
>>iancmc+wS2
You can even use python or solidworks' own api to do it programmatically
◧◩◪◨
30. iancmc+F73[view] [source] [discussion] 2023-10-21 06:26:54
>>ahmetr+fE
Deleting the bevel
◧◩◪◨⬒
31. dymk+Za4[view] [source] [discussion] 2023-10-21 17:09:15
>>ameliu+8Y
As far as I'm aware, we do not have the math needed to do what you're imagining. CSG requires a lot of numerical methods even for relatively simple operations.
◧◩◪◨⬒⬓⬔⧯
32. namibj+es4[view] [source] [discussion] 2023-10-21 19:12:43
>>auggie+0N
I can't wait for Idris2 (or such? It's quite good for programming and having e.g. core data structure implementations verified to uphold their API invariants, without making you prove everything around it.) to get a good fine tuned LLM attached to replace the human in the interactive theorem proving steps/parts.

I'd expect this to be an LLM application that's unusually suited to automatic-feedback reinforcement learning.

Especially because much of the interactive proving task/process is quite straightforward nondeterministic (Turing) machine stuff: You try (with enough randomization/temperature to not get stuck in non-creative stupidity) a prove attempt/step/hint, and get feedback after a moment of calculation from the proof assistant. Then you try further, until eventually hopefully getting "success" as the assistant's feedback.

Once you got it to succeed once, or seeded with originally human-made step sequences to teach some basic sense into the model, you know an upper bound on the _required_ step count and assistant calculation time to prove the theorem at hand. Thus you can let the LLM auto-play/train with a step limit and computation timeout close to the known upper bound, rewarding for expected total runtime/effort to combine "spamming cheap proof tactics and hope something sticks" and "elaborate careful proof process likely to succeed but always expensive/(semi-exhaustive) to go through".

Perhaps even with a GPT-4 like multi-agent LLM to specialize into the various approaches and have some way of rating/predicting each agent's expected efficiency/cost each "chat message":

Turns out, interactive theorem proving is literally a (beyond-)NP heuristics-guiding sampler (traditionally a human with trained gut feeling based problem-solving brainstorming creativity) chatting with a non-creative algorithmic oracle. At the start, if initiated not by the "human", the oracle would info-dump the theorem and appropriate context along a description of what "today's" task is:

This may be: 1) "I expect the theorem to be False because: [reason for what caused the expectation]." 2) "I expect the theorem to be True because: [reason for what caused the expectation]."

3) "I expect a weaker (in it's implications, so less general) form of that theorem to be sufficient in the proof for this situation here. It could make proving them for data types we'd like to use with this implementation/specialization cheaper than the general API's demands, potentially (though not the reason for today's task!) allowing weaker data types to be used here than in the general API. In particular, there (maybe) are cheaper (to implement) yet weaker (in effects) variants of the functions we call on the supplied data type, that still suffice in our algorithm (like weaker demands on a comparison function when only stable sorting is used; dropping distinctiveness demands for either hash function or comparison function if only the other is relied on for Set/Map data structure efficiency). Confidence of it to be True is at X (confidence measure/score the LLM can feel natively) and here's what certain outcomes would be worth (either enumerated weaker forms or a scoring function (in a form the LLM agent comparer can utilize to plan what to aim for, when to pivot, and when to declare defeat))."

4) "I expect a weaker (in it's demands, so more general) form of that theorem to be sufficient in the proof for this situation here. It would allow us to use this implementation/specialization with more data types than the general API promises. Confidence of it to be True is at X (confidence measure/score the LLM can feel natively) and here's what certain outcomes would be worth (either enumerated weaker requirements for the theorem or a scoring function (in a form the LLM agent comparer can utilize to plan what to aim for, when to pivot, and when to declare defeat))."

5) "Find (and codify!) sufficient invariants demanded from (functions on) data types so we can uphold our API's promised invariants. [Optionally, aim for something that matches this natural language description of what (parts/aspects of) the data type's function's purpose[s]/semantics are supposed to be.]"

6) "Prove this is constant-time/constant-memory-access-pattern w.r.t. that part of data (even in ways where the data in question may be in chunks behind some memcpy-reducing indirection), or e.g. that this key here affects nothing persistent other than that ciphertext/plaintext/signature/hash/success-flag."

7) "Prove time/memory complexity of this implementation. Here's what certain bounds are worth: split between various kinda of bounds, e.g. lower bounds, upper bounds, average (w.r.t. that data) complexity, best/worst case (w.r.t. that data), handle those parameters by enumerating over those and proving with their values fixed (because a general equation may be too complicated/weak, or even just too hard to derive)."

8) A classic: "Prove these two implementations produce identical (under that comparison/test-sampling method) results."

9) "Find bounds on when (conditions) and/or how much (some appropriate supplied measure) these two implementations differ."

10) "Find a faster implementation along with proven limits on it's inaccuracy. Combining the two (+) dimensions of candidate quality from the obvious pareto frontier into a single number score is according to this: [formula in useful format for directing search/exploration]."

11) "Cough up an implementation limited to those numerical primitives there, along with proof of it complying with these accuracy requirements, for this implementation that uses (inherently computationally-unsuitable) real numbers. Speed/memory performance importance: [scoring function suitable for directing where to aim, when to stop, and when to give up]."

And afterwards, the "human" would ask/explore the oracle about context, suggest/try proof tactics, and in some cases write/transform code in both LLM-style and by commanding ("textbook"/library/archive, or even freshly written) rewrite rules.

That process could then be trained with reinforcement learning, even if intermediate states have no useful score function defined, as the presence of certain results after certain amounts of expended effort is directly useful as a score for/of the solver/agent itself. The multi-agent suggestions applicability/efficiency predictor/arbiter should be amenable to more normal (stochastic) backpropagation at a completed-chat granularity, as the final efficiency/score will be known, and if it were a perfect predictor, it'd have predicted that exact score for the entire time along the path that was taken. The intermediate predictions on how much effort the chosen agent's suggestion actually took to complete is also easily recorded for training the per-step cost predictions as a more fine-grained aspect of the final-score-when-taking-this-branching-path-now machinery.

◧◩
33. digdug+Vpd[view] [source] [discussion] 2023-10-24 19:29:40
>>jessfr+Mx1
Very interesting. Do you have information/documentation on the engine itself? I'm not seeing anything technical that explains what it is/isn't capable of, or how it relates to existing CAD kernels.

Super cool work either way! I'll be wishing you luck.

[go to top]