Re: [deal.II] implement weak formulation of biharmonic equation

2018-02-26 Thread Timo Heister
> "make run" have error, and I don't understand what it mean and how to do: Did you read the error message, namely: > You are requesting information from an FEValues/FEFaceValues/FESubfaceValues > object for which this kind of information has not been computed. What > information these objects co

Re: [deal.II] implement weak formulation of biharmonic equation

2018-02-26 Thread chucui1016
Dear Timo, Thank you very much! And I correct what you say, then "make" can go, but "make run" have error, and I don't understand what it mean and how to do: yucan@ubuntu-csrc-b224:~/boost/deal.II/examples/practice-4order-39$ make run [ 50%] Built target practice4order39 [100%] Run with Debug

Re: [deal.II] implement weak formulation of biharmonic equation

2018-02-26 Thread chucui1016
Dear Timo, Thank you very much! And I correct what you say, then "make" can go, but "make run" have error, and I don't understand what it mean and how to do: yucan@ubuntu-csrc-b224:~/boost/deal.II/examples/practice-4order-39$ make run [ 50%] Built target practice4order39 [100%] Run with Debug

Re: [deal.II] implement weak formulation of biharmonic equation

2018-02-26 Thread Timo Heister
for (unsigned int d=0; d hessian_phi_v = fe.shape_hessian_component(i, k, d); Mii += dx * (hessian_phi_v.trace * hessian_phi_v.trace); You are missing curly braces: the way you wrote it, the Mii assignment is not done inside the loop. -- Timo Heister http

[deal.II] implement weak formulation of biharmonic equation

2018-02-26 Thread chucui1016
Dear all, I want to solve biharmonic equation problem with interior penalty method like step-39, but I need to rewrite weak formulation. In step-39, it use functions like LocalIntegrators::Laplace::cell_matrix