Commit f6ac8e7d authored by Fabian's avatar Fabian
Browse files

Example 3 now computes the energy norm

parent 4f5b4f50
......@@ -276,7 +276,12 @@ int main(int argc, char **argv) {
std::cout << "\\> Problem solved. \n";
std::cout << std::setprecision(16);
std::cout << " Energy-norm of interior = " << ui.transpose()*Kh*ui << std::endl;
std::cout << " Energy-norm = " << ui.transpose()*Kh*ui
+2*ui.transpose()*Krhs*ub
+ub.transpose() * assembler.assembleMatrix(
boundaryFESpace, boundaryFESpace,
BuilderData<gridView_t, febasis_t>(gridView, febasis, w)
)*ub << std::endl;
return 0;
}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment