Dear Kyusik, So if I understand you question correctly, you're not sure how to compute the "grad_p" part of the RHS during assembly? This you would do using the get_function_gradients() <https://www.dealii.org/8.5.0/doxygen/deal.II/classFEValuesBase.html#a69c63f8be0311c8970253276ac7df138> function in FEValues or, even better (and by using Extractors), the get_function_gradients() <https://www.dealii.org/8.5.0/doxygen/deal.II/classFEValuesViews_1_1Scalar.html#aac98b8207ba9c750a30165fd24fe1e6e> in the FEValuesViews. There are a number of tutorials that demonstrate this functionality - step-9 <https://www.dealii.org/8.5.0/doxygen/deal.II/step_9.html> and step-15 <https://www.dealii.org/8.5.0/doxygen/deal.II/step_15.html>, for example, get this information without Extractors while step-44 <https://www.dealii.org/8.5.0/doxygen/deal.II/step_44.html> does the same using Extractors.
I hope this helps! Regards, Jean-Paul On Monday, August 21, 2017 at 5:58:36 AM UTC+2, [email protected] wrote: > > 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 x component of the equation > (Inv_K*u)_x+(grad_p)_x=x * y * (grad_p)_y Here subscript > means component of vector (vector_x : x component of this vector) > > In the y component of the equation > (Inv_K*u)_y+(grad_p)_y=x * y * (grad_p)_x > > where grad_p comes from previous solution(Let say I use Picard iteration.) > > Could you please let me know How I can do these things? > > 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
