search

robertsj / libdetran

deterministic transport utilities
MIT License
2 stars 2 forks source link

Steady state ROM #26

Closed Relzohery closed 3 years ago

Relzohery commented 4 years ago

This PR adds three main classes for the steady state ROM, namely;

  1. ROMBasis: This gets a user basis form a binary file and stores it in a dense matrix.
  2. OperatotProjection : This takes any matrix type (sparse, dense or shell) and basis and produce the reduced matrix(UTAU)
  3. ROMSolver: This takes user input, mesh, material, and operator type and solves the eigenvalue problem of the reduced system.

Each class has its unit test.

In the Solver test, I used 1D slab reactor and tried three operators (Diffusion, EnergyDependent and EnergyIndependent). The Diffusion results seems reasonable. For the EnergyDpendent operator, while the eigenvalues of the ROM is close to the true value, they don't match the one from the diffusion.

Relzohery commented 4 years ago

@robertsj I did not solve all the issues yet. Also, I can break this into three smaller PRs, if you want.

robertsj commented 4 years ago

I'll take a peak today or tomorrow. I've got an NRC conference call this afternoon and won't make a meeting with you. I hope F that I can be in the office.

From: Relzohery notifications@github.com Sent: Wednesday, July 15, 2020 10:06 AM To: robertsj/libdetran libdetran@noreply.github.com Cc: Jeremy Roberts jaroberts@ksu.edu; Mention mention@noreply.github.com Subject: Re: [robertsj/libdetran] Steady state ROM (#26)

@robertsjhttps://github.com/robertsj I did not solve all the issues yet. Also, I can break this into three smaller PRs, if you want.

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/robertsj/libdetran/pull/26#issuecomment-658820741, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AAH66G5Z6RAYWV4WLXEY5KTR3XAVXANCNFSM4O2TVGFQ.