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::interpolate_boundary_values to get the pairs of dof_index > and value, and then convert the result of that function into a vector. > > At first, the reason why I would like to make vector G by myself (without using VectorTools::interpolate_boundary_values) is that I would like to input calculated values (not just function of position) on the boundary. I would like to solve free boundary problem. So, the boundary condition would be calculated from sources(current in coils) and line integral of green function and so on.... Especially, the calculation of boundary values includes the line integral, that is why I want to use an assembly step. So, do you think it is hard to apply boundary condition without using VectorTools::interpolate_boundary_values...? 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 Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.