Author: fabremathieu
Date: Tue Apr 22 14:24:33 2014
New Revision: 4621
URL: http://svn.gna.org/viewcvs/getfem?rev=4621&view=rev
Log:
work in progress
Modified:
trunk/getfem/interface/src/gf_compute.cc
Modified: trunk/getfem/interface/src/gf_compute.cc
URL:
http://svn.gna.org/viewcvs/getfem/trunk/getfem/interface/src/gf_compute.cc?rev=4621&r1=4620&r2=4621&view=diff
==============================================================================
--- trunk/getfem/interface/src/gf_compute.cc (original)
+++ trunk/getfem/interface/src/gf_compute.cc Tue Apr 22 14:24:33 2014
@@ -517,23 +517,27 @@
#ifdef EXPERIMENTAL_PURPOSE_ONLY
- /*@FUNC E = ('error estimate nitsche', @tmim mim, @scalar GAMMAC, @scalar
GAMMAN)
+ /*@FUNC E = ('error estimate nitsche', @tmim mim, @scalar GAMMAC, @scalar
GAMMAN, @scalar lambda, @scalar mu, @scalar gamma0,, @scalar f_coeff)
Compute an a posteriori error estimate in the case of Nitsche method.
Currently there is only one which is available: for each convex,
the jump of the normal derivative is integrated on its faces.@*/
sub_command
- ("error_estimate_nitsche", 3, 3, 0, 1,
+ ("error_estimate_nitsche", 5, 5, 0, 1,
const getfem::mesh_im &mim = *in.pop().to_const_mesh_im();
scalar_type GAMMAC = in.pop().to_scalar();
scalar_type GAMMAN = in.pop().to_scalar();
+ scalar_type lambda = in.pop().to_scalar();
+ scalar_type mu = in.pop().to_scalar();
+ scalar_type gamma0 = in.pop().to_scalar();
+ scalar_type f_coeff = in.pop().to_scalar();
unsigned si = unsigned(mim.linked_mesh().convex_index().last_true()+1);
darray err =
out.pop().create_darray_h(si);
getfem::base_vector ERR(si);
getfem::base_vector UU(U.real().size());
gmm::copy(U.real(), UU);
- getfem::error_estimate_nitsche(mim, *mf, UU, ERR);
+ getfem::error_estimate_nitsche(mim, *mf, UU, GAMMAC, GAMMAN, lambda,
mu, ERR);
gmm::copy(ERR, err);
);
_______________________________________________
Getfem-commits mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-commits
