[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] r4609 - /trunk/getfem/interface/src/gf_compute.cc
From: |
mathieu . fabre |
Subject: |
[Getfem-commits] r4609 - /trunk/getfem/interface/src/gf_compute.cc |
Date: |
Mon, 14 Apr 2014 12:38:09 -0000 |
Author: fabremathieu
Date: Mon Apr 14 14:38:09 2014
New Revision: 4609
URL: http://svn.gna.org/viewcvs/getfem?rev=4609&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=4609&r1=4608&r2=4609&view=diff
==============================================================================
--- trunk/getfem/interface/src/gf_compute.cc (original)
+++ trunk/getfem/interface/src/gf_compute.cc Mon Apr 14 14:38:09 2014
@@ -514,6 +514,33 @@
getfem::error_estimate(mim, *mf, U.cplx(), err, mim.convex_index());
);
+
+#ifdef EXPERIMENTAL_PURPOSE_ONLY
+
+ /address@hidden E = ('error estimate nitsche', @tmim mim, @scalar GAMMAC,
@scalar GAMMAN)
+ 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 address@hidden/
+ sub_command
+ ("error_estimate_nitsche", 3, 3, 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();
+ darray err =
+ out.pop().create_darray_h
+ (unsigned(mim.linked_mesh().convex_index().last_true()+1));
+ if (!U.is_complex())
+ getfem::error_estimate_nitsche(mim, *mf, U.real(), err,
mim.convex_index(), GAMMAC, GAMMAN);
+ else
+ getfem::error_estimate_nitsche(mim, *mf, U.cplx(), err,
mim.convex_index(), GAMMAC, GAMMAN);
+ );
+
+#endif
+
+
+
+
/address@hidden E = ('convect', @tmf mf_v, @dvec V, @scalar dt, @int nt[,
@str option[, @dvec per_min, @dvec per_max]])
Compute a convection of `U` with regards to a steady state velocity
field `V` with a Characteristic-Galerkin method. This
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-commits] r4609 - /trunk/getfem/interface/src/gf_compute.cc,
mathieu . fabre <=