Re: [deal.II] computing some other output quantity after solution is calculated.

2016-11-01 Thread hanks0227
Thank you very much for your reply It seems that I can solve the problems I have now. Thank you. Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you ar

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-11-01 Thread Daniel Arndt
Kyusik, [...] > 3. I calculated l2 norm > Difference = solution; > Difference -= exact_solution; > norm = Difference.l2norm(); > > But, according to what you mentioned, I calculated l2norm of U_i-U_i_exact > not l2norm of u_h(x)-u_h_exact(x) > Yes, this gives you the l2-norm of U_i-U_i_exact

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread hanks0227
Dr, Bangerth, > >u_h(x) = \sum_i U_i \varphi_i(x) > > First of all, I'm really sorry to bother you continuously... Could I have please one more question? I also wanted to calculate l2 norm of |u_h(x)-u_h_exact(x)| (where u_h_exact(x) is exact solution of u_h(x)). So what I did is fo

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread Wolfgang Bangerth
On 10/31/2016 08:11 PM, hanks0...@gmail.com wrote: But, my question was that in usual tutorial program of deal.II, the vector solution contains values for U_i not u_h(x). We use them interchangeably. The *solution vector* that comes out of the linear solver is U_i, and it corresponds to the *

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread hanks0227
> > > Well, both are the solution -- one is a function of x, the other is a > vector > of nodal values, but they ultimately represent the same object. U_i is the > vector you get out of the linear solver. > I'm sorry to make you confused I should have been more specific. But, my question was

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread Wolfgang Bangerth
This is one way to do it. You are just projecting J_tor to a finite element space and use that for output or further processing. As others have already noted, the only problem I can see is that computing solution_max may be inaccurate. This is because in general, if you have

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread hanks0227
Dr. Bangerth, Thank you very much for your kind reply. > This is one way to do it. You are just projecting J_tor to a finite > element space and use that for output or further processing. As others > have already noted, the only problem I can see is that computing > solution_max may be inacc

Re: [deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread Wolfgang Bangerth
Kyusik, I'm trying to calculate J_tor= x*solution+solution_max/y solution is already calculated vector from 3 function (setup_system, assemble_system, solve) To do this calculation I added the following three function(setup_Jtor, assmeble_Jtor, solve_Jtor) [...] This is one way to do it. Y

[deal.II] computing some other output quantity after solution is calculated.

2016-10-31 Thread hanks0227
Hi All, I'm trying to calculate J_tor= x*solution+solution_max/y solution is already calculated vector from 3 function (setup_system, assemble_system, solve) To do this calculation I added the following three function(setup_Jtor, assmeble_Jtor, solve_Jtor) template void Step6::setup_Jtor