Re: [deal.II] Reuse of the system_matrix in the Picard iteration for Inhomogenous Dirichlet boundary conditions

Oh, thank you very much!!! I was so stupid!! I thought that "system_matrix_tmp = system_matrix" is the same as "system_matrix_tmp.copy_from(system_matrix)" Thank you. Kyusik. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.goog

Re: [deal.II] segmentation error in refine_grid()

> > OK, thank you for your help. > 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 are subscribed to the Google Groups "deal.II User Group" group.

Re: [deal.II] segmentation error in refine_grid()

Dr. Bangerth, Thank you for your answer. But, could I have one more question...? Do you think this problem is related to "min_grid_level=0" I solve my problem just by changing from "min_grid_level=0 & max_grid_level=3" to min_grid_level=1 & max_grid_level=4" But I can't understand why thi

[deal.II] segmentation error in refine_grid()

Dear all, I am trying to apply deal.II for my research. But, the segmentation fault happens in the refine_grid() function that comes from step-26.cc tutorial example. The reason why I think the segmentation error occurs in this function are as follows... template void Step6::refine_grid (c

[deal.II] Re: Acess to Each component of RHS

Thank you very much for your help!! It seems that I can solve my problem now. 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 are subscribed to the

[deal.II] Re: Acess to Each component of RHS

Dear Jean-Paul Thank you for your kind answer. But, could I have one more question...? In step-9, "RightHandSide::value" and "RightHandSide::value_list" have to get "component" as input. I want to know the role of it exactly so, please check if what I am thinking is right or not. Th

[deal.II] Acess to Each component of RHS

Dear All, I'm studying step-20 tutorial example. In step-20, RHS of the first equation (that is Inv_K*u+grad_p=0) is zero. But, I want to make it non-zero function. So, I wonder whether it is possible for me to give the non-zero value to each component of RHS vector. For instance In the

[deal.II] Re: "operator+=" error in the assemble_system()

Dear Jean-Paul, > Firstly (and just so that we're all clear on what you're trying to > implement), was I correct that your comment that I linked contained the > description of the integral that you are trying to compute? > > Yes it is. that is something I have to do. > No, this is not true

[deal.II] Re: "operator+=" error in the assemble_system()

Dear Bruuno and Jean-Paul First of all thank you for your reply I'm sorry, I should have been more specific about my question As far as I know fe_values2.shape_value(i,q_index) is vector-valued shape function(fe2 is distributed by "dof_handler2.distribute_dofs (fe2)") and sol_grad is also

[deal.II] "operator+=" error in the assemble_system()

Dear All, I'm trying to use two FE_Q in my assemble_system() to project the gradient of scalar solution onto the vector-valued FE space (that is, one is for scalar solution and the other is for vector) my assemble_q() is as follow... template void Step6::assemble_q () { const QGauss quadra

Re: [deal.II] Error in the second derivative of solution...?

Thank you very much!! It seems that I can solve my problem. 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 are subscribed to the Goog

Re: [deal.II] Error in the second derivative of solution...?

Dear Jean-Paul, Thank you again for a detailed explanation. However, please excuse my ignorance... > > 1. Create a vector-valued finite element system (maybe an > FESystem(FE_Q(1),dim) ) and use it to distribute DoFs onto a new > DoFHandler. > 2. Create a mass matrix "M" with the new DoFHandl

Re: [deal.II] Error in the second derivative of solution...?

Dear Jean-Paul Pelteret, Thank you very much for your help. It is very helpful. > However, you'll note that this function only makes sense (and its > therefore only defined) for a vector field. Your solution field is scalar > though, which presents a bit of an issue. On the face of it you'd ha

Re: [deal.II] Error in the second derivative of solution...?

Dr. Wolfgang Bangerth, Thank you very much for your kind reply. But, could I have one more question? > > I don't know whether the problem really comes from that, but you can't > take > the second derivative of a finite element solution without problem. The > issue > is that it exists, of

Re: [deal.II] How to make G vector in lecture 21.65

Dr. Bangerth, Thank you for the kind explanation. Now, I understand what is wrong... I'm really sorry but, could I have one more question? > > The approach you use, making things look like an assembly step, is not > appropriate for what you want to do. I would just call > VectorTools::int

Re: [deal.II] How to make G vector in lecture 21.65

Dr. Bangerth, Thank you very much for your kind reply. > What is the content of 'constraints' in your case? > > DoFTools::make_hanging_node_constraints (dof_handler, constraints); VectorTools::interpolate_boundary_values (dof_handler,

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

Dear Jean-Paul Sorry for my belated reply. I was solving other problems... Anyway, thank you for your reply. I will study what you mention. thank you Kyusik. 2017년 3월 3일 금요일 오후 3시 15분 18초 UTC+9, Jean-Paul Pelteret 님의 말: > > Dear Kyusik, > > >> /user2/hanks318/dealii/dealii_pre/examples/ste

Re: [deal.II] Re: Lecture 31.65(Picard iteration)

Thank you for your relpy > But you could try to just set matrices and vectors to zero, and see > whether that makes a difference. > > As you mentioned, in "assemble_system()", I replace "system_rhs.reinit" and "solution.reinit" with for (unsigned int i=0;ihttp://www.dealii.org/ For maili

[deal.II] Re: Lecture 31.65(Picard iteration)

I'm sorry, please let me ask you one more question > Why is that expected? It should be the same. The new sparsity_pattern > should be the same as the old one > > > According to your comment, I remove system_matrix.reinit(sparsity_pattern) only (system_rhs.reinit and solution.reinit are no

[deal.II] Re: Lecture 31.65(Picard iteration)

Thanks for your reply please let me check that I understand what you mean. > The new sparsity_pattern should be the same as the old one and you still > have the same number of degrees of freedom. That's why there is no need to > do the reinit. > As I mentioned, the results from two cases ar

[deal.II] Lecture 31.65(Picard iteration)

Hi all, I have a question in the lecture 31.65 that is about Picard iteration. In the lecture, the professor doesn't include system_matrix.reinit (sparsity_pattern), system_rhs.reinit (dof_handler.n_dofs()), and solution.reinit (dof_handler.n_dofs()) in assemble_system(). However, I was wor

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

Thanks for your reply again. > Yes, of course. Via VectorTools::interpolate_boundary_values and > VectorTools::project_boundary_values. > > Due to your help, I can apply ZeroBoundaryCondition for pressure by adding //test for strong boundary condtion std::map boundary_values_strong;

Re: [deal.II] ZeroBoundaryCondition when FESystem is used

Thank you very much your kind reply. Could I have one more question..? You mention that > > The difference that in step-3, you are imposing *strong* boundary > conditions, > whereas in step-20, you are imposing *weak* boundary conditions. > > Then, is there anyway to use strong boundary cond

[deal.II] Re: ZeroBoundaryCondition when FESystem is used

I'm sorry, I should have include "part for boundary of pressure" What I mean by "part for boundary pressure in assemble_system" is (in step-20) for (unsigned int face_no=0; face_no::faces_per_cell; ++face_no) if (cell->at_boundary(face_no))

[deal.II] ZeroBoundaryCondition when FESystem is used

Hi all, I'm studying step-20 tutorial example. But, I can't understand it, because it seems to have different things compared to step-3 in terms of boundary condition. In the introduction of step-3, it is said that ∫Ω∇φ⋅∇u−∫∂Ωφn⋅∇u=∫Ωφf. The test function φ has to satisfy the same kind of bou

Re: [deal.II] second derivatives of solution

OK Thank you very much! Kyusik. 2017년 2월 13일 월요일 오후 1시 3분 27초 UTC+9, Wolfgang Bangerth 님의 말: > > On 02/12/2017 06:15 PM, hank...@gmail.com wrote: > > > > Do you mean that it is right way to get second derivatives but it could > not be > > continuous? > > Yes, you are using the correct way

[deal.II] second derivatives of solution

Hi all, It seems that my question is very easy one, but I can't find the similar thing... I know how I can use solution and first derivative solution that comes from fe_values.get_function_values and fe_values.get_function_gradients So I thought that second derivative can be obtained by

Re: [deal.II] Duplication access to the cell?

Thank you very much!!! The problem is solved, due to your help!! Kyusik. 2016년 12월 14일 수요일 오전 1시 55분 53초 UTC+9, Wolfgang Bangerth 님의 말: > > Kyusik > Material ids are only stored as unsigned char objects, which can only > represent values up to 255. You may want to use the user_index instead. >

[deal.II] Duplication access to the cell?

Dear All, I'm trying to calculate something differently depending on material id. So, I have to set material id before next iteration.(I'm going to use picard iteration) The way I have done is 1.the middle point of mesh would be selected. (it is indicated as max_cell in my code) 2.the neighb

[deal.II] Re: Distorted cell when it is crated by gmsh

Oh, I didn't know that sorry. Thank you all for your replies. 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 are subscribed to the Google Groups

Re: [deal.II] Accesing cell from cell index

Dear Pelteret, I'm sorry to have bothered you... I should have been more considerate before asking something. 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 t

Re: [deal.II] GridTools::laplace_transform and interior points

OK, I see. Anyway, as always, Thanks for educating me free of charge :-) 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 are subscribed to the Googl

Re: [deal.II] GridTools::laplace_transform and interior points

Dr. Bangerth, Thanks for your reply again. > Yes, but that's because you have a small bug in your code: you fix *all* > vertices before you call laplace_transform. You should only fix the ones > that > are are at the boundary of the domain. > > As you said above, I only touched the points on

Re: [deal.II] GridTools::laplace_transform and interior points

Thank you very much for your quick reply. > > I'm not quite sure what you're getting at here. The mesh you show looks > correct to me. What is it that you're surprised about? > > > I'm sorry, I should have been more specific. What I mean that the result has no difference compared to when I didn

[deal.II] Re: How to change material ids during iteration.

Thanks for your reply. you are right, I shouldn't have been in a hurry to show some results to my boss... I should have asked, after looking through several examples especially step49 that includes GridTools::laplace_transform. I guess, using GridTools::laplace_transform, I can solve my proble

[deal.II] Re: How to change material ids during iteration.

First of all, Thank you very much for your kind reply. > you mean physically move the triangulation to change the spatial position > of the interface between two different materials, then this is also > possible. Look at step-18 >

[deal.II] How to change material ids during iteration.

Hi all, I'm trying to solve 2 materials problem(Vacuum and Plasma). And the boundary between these two materials must be defined as contour of minimum value of previous solution(I'm going to use picard iteration) I can input initial guess for this boundary(and material ids) by using GMSH. But,

Re: [deal.II] maximum value of solution

Oh, I see. Thank you very much for your kind reply. 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 are subscribed to the Google Groups "deal.II Us

[deal.II] Re: Visit error...?

Oh, I see I'm sorry I thought about it in the wrong way. Thank you very much your kind reply. 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

[deal.II] Re: Visit error...?

Arndt, First of all, Thank you very much for your reply. > What do you mean by "average value" here? You probably just wnat to > compute the error in the L2-norm, isn't it? > I'm sorry this is not average, it is RMS value. In deal.ii, from step7, L2 norm is calculated like this... L2_n

[deal.II] Re: Visit error...?

I'm sorry, I missed the sqrt(#cell) in the calculation of L2norm I did std::cout<<"l2norm(error)="

Re: [deal.II] Absolute values of deal.II vector

Dr. Bangerth, Thank you very much. It works. 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 are subscribed to the Google Groups "deal.II User Gr

[deal.II] Absolute values of deal.II vector

Hi All, First of all, my question is very simple. But I can't find similar question. Anyway, what I'm trying to do is to get error = |exact_solution - numerical_solution| (absolute value of two Deal.II Vector) and then plot this one using Visit. But I don't know how to change all elements of e

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

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.

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.

> > > 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] mesh generation for arbitrary boundary

Dr. Bangerth, I'm sorry... I found the similar question about this problem. That's why I eliminate my question... Anyway Thank you very much for your reply. Kyusik. 2016년 11월 1일 화요일 오전 2시 22분 58초 UTC+9, Wolfgang Bangerth 님의 말: > > On 10/31/2016 04:00 AM, hank...@gmail.com wrote: > > > > I

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

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

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

Daniel, > If you know that your FiniteElement is interpolating (which I assume given > the way you compute the maximum) and that J_tor= x*solution+solution_max/y > is in your ansatz space (probably not) or an interpolation is sufficient, > you can use something like > > J_tor(local_dof_indices

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

> First of all, Thank you very much for your reply > > This can only work if you know that the maximum is attained at one of the > support_points > of your FiniteElement. In general, you can use > VectorTools::integrate_difference to find a good approximation > to the maximum value of your

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

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

Re: [deal.II] how to use solution.

2016년 10월 5일 수요일 오후 6시 55분 19초 UTC+9, Daniel Arndt 님의 말: > > Kyusik, > > As Wolfgang already mentioned, you want to use > FEValues::get_function_values() > to extract from a global FE Vector the local values at the quadrature > points in the cell you are on. > These quadrature points normally

Re: [deal.II] how to use solution.

2016년 10월 5일 수요일 오전 1시 17분 31초 UTC+9, Wolfgang Bangerth 님의 말: > > On 10/04/2016 09:12 AM, hank...@gmail.com wrote: > > > > template > > void Step6::get_profile () > > { > > const QGauss quadrature_formula(3); > > FEValues fe_values (fe, quadrature_formula, > >

Re: [deal.II] how to use solution.

2016년 10월 5일 수요일 오전 1시 17분 31초 UTC+9, Wolfgang Bangerth 님의 말: > > On 10/04/2016 09:12 AM, hank...@gmail.com wrote: > > > > template > > void Step6::get_profile () > > { > > const QGauss quadrature_formula(3); > > FEValues fe_values (fe, quadrature_formula, > >

[deal.II] how to use solution.

Hi, I'm trying to study how to use solution after the calculation is done. So, What I'm trying to do in step6 tutorial program is to calculate f=x*solution+y; So I add the following class.(that almost comes from assemble_system()) template void Step6::get_profile () { const QGauss quadratu

Re: [deal.II] Re: How can I have access to the all element of previous solution?

> > Thank you very much for your kind reply!!! Best, 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 are subscribed to the Google Groups "deal.I

[deal.II] Re: How can I have access to the all element of previous solution?

First of all, Thank you very much for your reply And, I'm sorry, I should have been more specific about my problem In the lecture 31.65, the professor change the step-6.cc a little to apply Picard iteration. Doing this, he uses previous_solution_gradients[q_index]. So, like this, I just

[deal.II] How can I have access to the all element of previous solution?

Hi, I'm studying Wolfgang's lecture 31.65 that deals with picard iteration. In the lecture, previous_solution_gradients[] is used. But, I also want to know How I can use just previous_solution[]. At first, I used previous_solution[q_index] in assemble_system as in previous_solution_gradients[

[deal.II] Re: After the calculation, Can the result be used??

Thank you very much!!! Your answers would be very helpful I will study what you mention and come back with more specific questions Thanks Best regards -- 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

[deal.II] Re: After the calculation, Can the result be used??

I -- 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 are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop rec

[deal.II] Re: After the calculation, Can the result be used??

I -- 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 are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop rec

[deal.II] Re: After the calculation, Can the result be used??

I -- 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 are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop rec

[deal.II] Re: After the calculation, Can the result be used??

Thank you very much!!! I will study what you mention. Thanks Best regards 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 are subscribed to the Googl

[deal.II] After the calculation, Can the result be used??

Hi I'm Kyusik. I just wonder if it is possible that The result of the calculation can be used for the other calculation. For example, In step5, The Eq is -div(alp(x)*grad(u(x)))=f with appropriate Boundary Condition. After this calculation, is it possible for me to use the solution that is u

Re: [deal.II] singularity error due to 1/0

Thanks for your answer. But, I still don't know how I can deal with this error. Anyway, Thank you very much. Kyusik. 2016년 9월 5일 월요일 오후 9시 47분 21초 UTC+9, Wolfgang Bangerth 님의 말: > > On 09/05/2016 06:04 AM, hank...@gmail.com wrote: > > > > K_Inv[0][0]=2*r/(cos(theta)*cos(theta)), > > K_Inv[

[deal.II] Re: singularity error due to 1/0

I have tried to solve this problem including what you suggested. But,I'm sorry to tell you, Everything didn't work out well... It seems that csc(theta) or sec(theta) can't be used in cylindrical coord system But, I solve GradShaf Eq that I was originally trying to solve using step-6 Anyway, T

[deal.II] Re: singularity error due to 1/0

Oh, sorry my mistake Not theta=x/y in my code but theta=atan(x/y) 2016년 9월 5일 월요일 오후 9시 22분 49초 UTC+9, Jean-Paul Pelteret 님의 말: > > Hi Kyusik, > > r is calculated by "r=sqrt(x^2+y^2) ", and theta is calculated by >> "theta=x/y" > > > At the very least, your definition of theta in the (radial)

Re: [deal.II] Re: applying 3D mesh in step-20

Thank you very much!! 2016년 9월 5일 월요일 오전 4시 38분 52초 UTC+9, Wolfgang Bangerth 님의 말: > > On 09/04/2016 05:38 AM, hank...@gmail.com wrote: > > > > I have to show my result to my boss(variable p). And, To show the inside > of > > torus, It seems that I have to use contour plot using Visit.(In the

[deal.II] Re: applying 3D mesh in step-20

Thanks for your Answer Bruno. But, It seems that I am using the old version. When I try to use Torus manifold, The Error message occurs. And, the problem is solved by making new mesh. Anyway, Thank you very much for your Answering my Question. Best. Kyusik 2016년 8월 31일 수요일 오후 11시 55분 34초 UTC

[deal.II] Re: applying 3D mesh in step-20

Oh, Thank you very much! The problem is solved!!! you are right. My mesh was the problem. I was stupid!!! But, Could I have the last question? I have to show my result to my boss(variable p). And, To show the inside of torus, It seems that I have to use contour plot using Visit.(In the

[deal.II] Re: applying 3D mesh in step-20

First of all, Thank you very much for your answering my Question. Could I have one more question...? Then, Is it okay If I use Prism cell instead of the Hexahedron cell in UCD format?? Can the deal.II treat the Prism cell?? And, Do you also agree with that "SphericalManifold boundary" is not r

[deal.II] applying 3D mesh in step-20

Hi I'm studying step-20 tutorial example, watching Wolfgang's Lectures In step20, I would like to solve that Eq in the 3D mesh that I made by myself(I attached B_coord3D.inp that is UCD format) But, It didn't work out. The error message is like this [ 50%] Built target step-20 make[3]: W

Re: [deal.II] Question step5 in the lectures(about Coefficient...)

Thank you very much indeed!! 2016년 7월 6일 수요일 오후 3시 57분 51초 UTC+9, Praveen C 님의 말: > > Hi > > In value_list function, you can get (x,y) like this > > double x = points[i][0]; > double y = points[i][1]; > > Best > praveen > > On Wed, Jul 6, 2016 at 12:13 P

[deal.II] Question step5 in the lectures(about Coefficient...)

Hi, I am trying to study step-5.cc. But I have some problem... I would like to make Coefficient 1/x (As you know, In step-5.cc, The Coefficient is originally 20(if r<0.5) or 1(if r>0.5)) of course not on the circle grid(circle-grid.inp), but on the other grid not including the point that x=0