Thank you, David.
I also thought about changing the return value of the “Real RBEvaluation::rb_solve(unsigned int N),” but I did not edit it because of the fear of an error. Since you said there is no other way, I will try to add the code in “rb_construction.C” or “rb_evaluation.C” instead of changing the return value. I always appreciate your help. Best regards, SKang From: David Knezevic <david.kneze...@akselos.com> Sent: Monday, January 21, 2019 10:05 PM To: 강신성 <ss.k...@pusan.ac.kr> Cc: libmesh-users <libmesh-users@lists.sourceforge.net> Subject: Re: [Libmesh-users] [RB] How to print output error bounds in the offline stage (RB training) On Sun, Jan 20, 2019 at 10:46 PM <ss.k...@pusan.ac.kr <mailto:ss.k...@pusan.ac.kr> > wrote: Thank you for a reply, David. I know the code, "libMesh::out << "Maximum error bound is " << training_greedy_error << std::endl" in rb_construction.C in order to print error bounds during the RB training. However, the "training_greedy_error" is a "Maximum solution error bound," not a "Maximum output error bound." If you look at the "Real RBEvaluation::rb_solve(unsigned int N)," variables associated with "Maximum solution error bound" and "Maximum output error bound" are as follows. ======================================================================================== * "abs_error_bound" at line 272 >> This is the "solution error bound," which is already in print at every step >> of RB training. * "RB_output_error_bounds" at line 277 >> This is the "output error bound," the error I want to print at every step of >> RB training. ===================================================================================== In other words, my question is as follows: Can I print a "Maximum output solution error bound" in each basis dimension? I look forward to hearing from you. Oh OK, I see, you want to print the output error bounds. This should be possible, but the code isn't set up to do that at the moment. I think you'd have to edit rb_construction.C in order to achieve this. You should be able to call the same code that generates the output error indicator during an online solve, and use that to print out the output indicator during training, but I'm afraid that I don't have time to look into the details myself at the moment. Best regards, David From: David Knezevic <david.kneze...@akselos.com <mailto:david.kneze...@akselos.com> > Sent: Friday, January 18, 2019 11:57 PM To: 강신성 <ss.k...@pusan.ac.kr <mailto:ss.k...@pusan.ac.kr> > Cc: libmesh-users <libmesh-users@lists.sourceforge.net <mailto:libmesh-users@lists.sourceforge.net> > Subject: Re: [Libmesh-users] [RB] How to print output error bounds in the offline stage (RB training) On Fri, Jan 18, 2019 at 1:37 AM <ss.k...@pusan.ac.kr <mailto:ss.k...@pusan.ac.kr> > wrote: Hello, all. I have a question about printing output error bounds in the offline stage. We can see that maximum solution error bounds of each step in the RB training process, such as ================================================================== . ---- Basis dimension: 1 ---- Performing RB solves on training set Maximum error bound is 1.07082 Performing truth solve at parameter: . ================================================================== Besides, using the code "rb_eval.RB_output_error_bounds[0]" in the main program, I printed the output error bound after the RB training. However, I cannot find any way to print all output error bounds during an RB training. In other words, I want to see output error bound corresponding the basis dimension similar to the above code. Could you please tell me how to print them? I always thank you for your help. I'm not really sure what you want to print... but the print statement you referred to above is in rb_construction.C, where it does "libMesh::out << "Maximum error bound is " << training_greedy_error << std::endl << std::endl;". I guess you could do something like that in your code, or add libMesh::out calls in your local libMesh in order to print more info? Best, David _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
