MYSTRAN 18.0.0

MYSTRAN 18.0.0

This major update brings massive performance improvements, a couple new features (and removal of some obsolete ones), and many bug fixes.

General performance improvements

The most noticeable change users will be able to notice is are massive performance improvements. We've seen speedups from 2x to 4x all the way to 10x.

This was achieved through optimisation of several loops and removing unused cruft. Props to @Copper280z and @victorkemp for their huge effort!

Multithreaded SuperLU support

Support has been added for the multithreaded variant of SuperLU. Some larger decks can run over twice as fast!

However, due to the fiddly nature of SuperLU's run-time memory scaling parameters, this has not been made into a default feature.

When building MYSTRAN, one can enable this feature with the -DUSE_SUPERLU_MT=yes. This requires OpenMP in your system, but most Linux and Windows (MSYS2) will already have it. If you do decide to use it, keep an eye on the SPIENV 6/7/8 PARAMs, especially SPIENV7. See PR #182 by @Bruno02468 for details.

F04 and ANS removed

The subroutine timing and log file (F04) is obsolete when debuggers and profilers exist -- and it required every source file to have a header and a footer.

The ANS file contained data similar to F06 but more machine-friendly. In the presence of OP2 and F06 tools (like pyNastran and nastools respectively), it can also be considered obsolete.

Both these formats required a ton of repetitive code for no tangible benefit, and @victorkemp undertook the grueling effort to rid us of them both.

The state of Lanczos...

This update brings several changes to the Lanczos eigenvalue extraction method. Two pieces of good news and one breaking change (warranting the major-version increase).

SuperLU for Lanczos

Thanks to @victorkemp, Lanczos solutions can now use SuperLU for the factorization step, and do so by default. It's a massive performance increase and memory usage decrease: decks that used to take hours and gigabytes now take minutes and megabytes.

Adaptive Lanczos driver for frequency-band mode

Frequency-band Lanczos mode (where you don't specify a fixed number of eigenvalues, just a frequency band) was essentially broken. The inertia-based calculation required a band matrix, which we no longer have due to the removal of BANDIT. It's not viable to run it for non-band matrices, hence the EIGESTL parameter, only allowing the estimator to run for decks with 5000 or less L-set DoF. Increasing this value just allowed the estimator to run for hours.

Instead, an adaptive Lanczos driver was implemented by @Bruno02468. Put very simply, it starts with 10 modes and doubles the amount until the number of in-band modes found becomes stable (or reaches the absurd limit of 10240). There's a bit more to it than that, but that's the gist.

What matters is, frequency-band mode now works for >5k L-set DoF.

Removed EIGRL continuation card

The EIGRL continuation card has been removed. It represented an incompatibility with major Nastran-compatible solvers, its fields were poorly-explained, and none of them should be set by the user.

Other fixes

There's also been many other bug fixes: K6ROT, MPCs, progress counters, timestamps...

Feel free to check the commits since 17.0.0 or the PRs below for details:

• Fixed existing cheap K6ROT that's dof6 springs to ground by @victorkemp in #175
• SuperLU for Lanczos by @victorkemp in #176
• Mass matrices by @victorkemp in #179
• Fix principal angle bug for tiny strains or stresses by @victorkemp in #181
• Multithreaded SuperLU support by @Bruno02468 in #182
• Remove EIGRL continuation card and LANCMETH parameter by @Bruno02468 in #183
• Fixing frequency-band mode: Adaptive Lanczos by @Bruno02468 in #186
• One buckling fix, one eigen fix by @Bruno02468 in #194
• Bugfix for 4-term MPCs by @victorkemp in #190
• Solid stress/strain in material coordinates by @victorkemp in #187
• Change PSOLID IN and ISOP options by @victorkemp in #108
• Buckling MPC fixes by @Bruno02468 in #195
• Lower DARPACK for small models by @Bruno02468 in #197
• reordered array dimensions by @victorkemp in #199
• converted internal grid number to sequence number by @victorkemp in #200
• Redefine get grid num comps by @victorkemp in #201
• gated writes to F06 by @victorkemp in #203
• Remove ans by @victorkemp in #205
• Use reverse lookup instead of nested loop for autospc by @Copper280z in #204
• Small counters update by @Bruno02468 in #207
• UNIX_TIME taking too much time by @Bruno02468 in #208
• Common link message by @victorkemp in #209
• Optimize matadd sss by @victorkemp in #210
• LINK9 speedup by @victorkemp in #212
• Remove f04 by @victorkemp in #213
• Removed SUBR_BEGEND_LEVELS from BLAS files by @victorkemp in #214
• added dev branch to build action by @victorkemp in #215
• Dev to main by @victorkemp in #216

MYSTRAN 17.0.0

MYSTRAN 17.0.0

This update has been in the works for a while, and several new features have been introduced.

Among the additions, the most relevant are the MITC shell and the QUAD8 element -- both courtesy of @victorkemp.

There have also been bug fixes and performance improvements.

For details, check out the list of pull requests below, or the commit-level diff between 16.0.0 and 17.0.0 if you prefer.

Pull requests

• Changed default AUTOSPC NSET from 3 to 1 by @victorkemp in #132
• shell differential stiffness implemented but a few loose ends by @victorkemp in #127
• Better counters, better build script, SuperLU update by @Bruno02468 in #131
• fixed rotated coordinate system bug by @victorkemp in #133
• Condition to skip K6ROT code if 0 to avoid crash by @victorkemp in #134
• Changed ALPVEC use for solid, membrane, composite by @victorkemp in #137
• Composite buckling by @victorkemp in #139
• Bailout fixes by @victorkemp in #140
• Mitc8 basic stiffness by @victorkemp in #144
• Mitc8 with stress/strain/force output by @victorkemp in #147
• MITC4/MITC4+ basic stiffness by @victorkemp in #149
• MITC4+ more features by @victorkemp in #153
• included -1/t^2 factor in MID4 output by @victorkemp in #154
• MITC4+ gravity and eigenvalue by @victorkemp in #155
• Updated release_notes.md in preparation for 16.0.0 release by @zchlrnr in #157
• Fatal error is PSHELL's MID4 is not 0 or blank by @victorkemp in #159
• Pre-v17 memory bug fixes by @Bruno02468 in #160
• Prevents crash but not sure if OP2 is correct by @victorkemp in #163
• Fix MITC4+ corner stress bug by @victorkemp in #168
• Fix missing default SIGMA and set to -1.0 by @victorkemp in #169
• Removed DEF EXT from INI file and made bdf and dat both defaults. by @victorkemp in #171
• Getting rid of fort.N files and the many issues they cause by @Bruno02468 in #174

MYSTRAN 16.0.0

This update contains several fixes, including breaking changes to the output.

Changes

• Runtime performance improvement for large problems by @Copper280z in #87
• Bandit fix by @ceanwang in #91
• fix issues #95 #96 #98 by @victorkemp in #100
• Fix issue #99 shell orthotropic thermal strain by @victorkemp in #101
• fix #102 PCOMP THETA by @victorkemp in #105
• fix TRIA3K nonsymmetric composite error message by @victorkemp in #106
• PARAM QUADTYP MIN4 default by @victorkemp in #107
• fix PARAM STR_CID for solids by @victorkemp in #111
• Fix issue #109 CBUSH long format by @victorkemp in #110
• Solid corner stress by @victorkemp in #113
• Fixed MCID option for CQUAD4 and CTRIA3 by @victorkemp in #117
• Q4SURFIT default value by @victorkemp in #116
• Corrected CHAN's I2,J; and I1's A,I1,I2,J by @victorkemp in #119
• fixed and enabled solid element differential stiffness by @victorkemp in #123
• fix two memory bugs by @Bruno02468 in #124

All commits since 15.2.1: 15.2.1...16.0.0

MYSTRAN 15.2.1

This update contains improvements to the input parser and fixes several crashes.

Changes

• Improvements to the input parser by @SteveDoyle2 in GH-58, GH-68 and GH-71
• Fix for a segfault caused by a closed ANS file by @SteveDoyle2 in GH-64
• Fix for broken UTF-8 in some output files by @Bruno02468 in GH-79
• Other minor bug fixes by @Bruno02468 and @SteveDoyle2 in GH-85

All commits since 15.2.0: 15.2.0...15.2.1

MYSTRAN 15.2.0

This update implements a significant improvement to the RBE3 element (which would hardly be possible without the invaluable assistance we got from Victor from the MecWay project), a TUBE2 option for PBARL props, and adds grid point forces in OP2 output.

If you're interested in technical details, see the pull-request-level information:

• Adding TUBE2 option to PBARL by @SteveDoyle2 in #40
• Adding GPFORCE op2 writing to Mystran by @SteveDoyle2 in #55
• The RBE3 Fix by @Bruno02468 in #59

MYSTRAN 15.1.7

This update fixes a crash found when running some buckling solutions. Check out GH-44 for details.

MYSTRAN 15.1.6

This patch integrates memory bug fixes for buckling decks. See GH-41 for details.

MYSTRAN 15.1.5

This update integrates four PRs by @SteveDoyle2:

• fix PARAM,GRDPNT bug in #21
• case control no longer throws warnings for PLOT unless it's GPFORCE in #23
• fixing op2 subcase id in #32
• op2 gets date correct now instead of faked value #35

All commits compared to previous release: 15.1.4...15.1.5

MYSTRAN 15.1.4

This update fixes a memory bug when writing BAR engineering forces to a .ANS file. See GH-14.

MYSTRAN 15.1.3

This patch-level update fixes a memory bug that could arise in models with PCOMPs. See GH-13.