search

svalinn / mcnp2cad

Convert MCNP input files to a general CAD format
BSD 2-Clause "Simplified" License
34 stars 31 forks source link

Explore matrix capability of CGM #28

Closed gonuke closed 8 years ago

gonuke commented 8 years ago

It may be possible to replace Armadillo with methods in AnalyticGeometryTool.[ch]pp and CubitMatrix.[ch]pp

pshriwise commented 8 years ago

After looking through those files, I don't see the capability to calculate the pseudo-inverse directly or via SVD :(

gonuke commented 8 years ago

CubitMatrix has an inverse() method for 3x3 matrices - isn't that all you need?

pshriwise commented 8 years ago

I specifically need the pseudo-inverse (aka Moore-Penrose inverse) to find the inverse that has a minimal translation magnitude (which manifests as an inverse with the solution to a linear least squares problem) for the GQ in order to get the translation correct for GQs with an infinite dimension (e.g. elliptical cylinders).

This could also be accomplished by using the singular value decomposition and then some matrix multiplication, but I don't see that in there either.

gonuke commented 8 years ago

OK - got it...

gonuke commented 8 years ago

BTW, I think Armadillo's license (Mozilla) will allow us to vendor it into the plugin code base. It's CMake build system should further support that.