Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Nonlinear elasticity projection rom example #17

Merged
merged 9 commits into from
Oct 13, 2023
Merged

Conversation

dreamer2368
Copy link
Collaborator

@dreamer2368 dreamer2368 commented Sep 20, 2023

Example commands with 3d mesh (as in librom.net)

python3 nonlinear_elasticity_global_rom.py --mesh ../data/beam-hex-nurbs.mesh -offline -dt 0.01 -tf 5.0 -s 14 -vs 10 -sc 3.90 -id 0
python3 nonlinear_elasticity_global_rom.py --mesh ../data/beam-hex-nurbs.mesh -offline -dt 0.01 -tf 5.0 -s 14 -vs 10 -sc 4.10 -id 1
python3 nonlinear_elasticity_global_rom.py --mesh ../data/beam-hex-nurbs.mesh -merge -ns 2 -dt 0.01 -tf 5.0
python3 nonlinear_elasticity_global_rom.py --mesh ../data/beam-hex-nurbs.mesh -offline -dt 0.01 -tf 5.0 -s 14 -vs 5 -sc 3.92 -id 2
python3 nonlinear_elasticity_global_rom.py --mesh ../data/beam-hex-nurbs.mesh -online -dt 0.01 -tf 5.0 -s 14 -vs 5 -hyp -rvdim 40 -rxdim 10 -hdim 71 -nsr 200 -sc 3.92

FOM result

tep 500, t = 5.000000, EE = 1.41649e-03, KE = 4.71724e-04, ΔTE = -2.91350e-03
Elapsed time for time integration loop 1.86668e+02
Elapsed time for entire simulation 2.10924e+02

ROM result

Elapsed time for time integration loop 3.47888e+01
Relative error of ROM position (x) at t_final: 5.000000 is 5.08981700e-03
Relative error of ROM velocity (v) at t_final: 5.000000 is 9.11743695e+00
Elapsed time for entire simulation 6.63521e+01

@dreamer2368 dreamer2368 marked this pull request as ready for review September 22, 2023 02:05
@dreamer2368 dreamer2368 added the RFR Ready for review label Sep 22, 2023
@dreamer2368
Copy link
Collaborator Author

@chldkdtn @siuwuncheung @bumpwy @pravija12345 , this PR is now ready for review.

@bumpwy
Copy link
Collaborator

bumpwy commented Oct 11, 2023

everything runs fine now with the docker (albeit slow). Although, the relative error for the ROM predicted velocity seems a bit high. Is that of concern?

@dreamer2368
Copy link
Collaborator Author

everything runs fine now with the docker (albeit slow). Although, the relative error for the ROM predicted velocity seems a bit high. Is that of concern?

I also noticed that. Although, this is consistent with the c++ libROM result, so it should be rather considered as the performance of the method itself.

# # delete A; # TODO: find a good solution for this.
# return S;

# TODO: move this to the library?
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think this is general enough to be used in other examples, so it may be a good idea to move it to the library.

@dreamer2368 dreamer2368 merged commit 258f178 into main Oct 13, 2023
5 checks passed
@ckendrick ckendrick mentioned this pull request Oct 18, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
RFR Ready for review
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants