OilfieldToolsNet/IFEM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity
2,006 Commits 1 Branch 1 Tag 15 MiB
4b341216c4
Commit Graph

2006 Commits

This Branch
This Branch
All Branches
Author SHA1 Message Date
Arne Morten Kvarving
dc48750a66 add hook for post-processing element norms 
use this for entries which are sums norms
and where you have to do the summation after the
square-root has been applied.
 2017-10-05 13:58:44 +02:00
Arne Morten Kvarving
9b13301cae increase maximum number of norms to 30 2017-10-05 13:58:44 +02:00
Arne Morten Kvarving
8a6382155d changed: check if an empty secondary solution is associated with residual norms 
if so, count the norm entries.

also add a flag for externally provided projections in the norm class.
this is necessary for applications where the projection is performed
on a different basis than than the first finite element basis. in particular,
this is used with subgrid and taylor-hood discretizations of stokes
 2017-10-05 13:58:44 +02:00
Arne Morten Kvarving
298f15af78 fixed: lower-dimension connections is not an error 
just ignore them in serial
 2017-10-05 11:30:25 +02:00
Knut Morten Okstad
e56e247492 Added: Implementation of the method to interpolate a function over a domain, 
defaults to L2-projection for ASMu2D and ASMu3D.
Changed: Merged to two edgeL2projection methods for ASMu2D into one method
using the FunctionBase interface. Also, don't use the FiniteElement class
in the L2-projection of functions on boundaries, since no Integrand here.
Other changes are cosmetics, doxygen fixes, etc.
 2017-10-04 13:52:16 +02:00
Knut Morten Okstad
dcb433e1d0 Changed: Not likely you have log and vtf-files in your working directory 
of this repository after the out-splitting of the Apps in separate repos.
Therefore removed. Also removed incomplete ignore of app-directories.
 2017-10-04 13:52:16 +02:00
Arne Morten Kvarving
20a1ff8ec6 update README with umfpack info 2017-10-04 08:50:36 +02:00
Kjetil Andre Johannessen
d3bab37cd1 Bugfix: 3D-LR projection (CGL2) bezier evaluation was wrong 2017-09-29 11:22:27 +02:00
Arne Morten Kvarving
9329ac73f6 add umfpack solver to SparseMatrix 2017-09-28 22:10:56 +02:00
Arne Morten Kvarving
5c9bac9cf9 fixed: initial guess mess with preonly ksp in PETSc (lu solves) 2017-09-28 22:10:56 +02:00
Arne Morten Kvarving
20ecda0d55 changed: print CFL for negative maxCFL value 
this allows to print but not apply time stepping control
 2017-09-22 11:48:17 +02:00
Knut Morten Okstad
27d4538482 Added: Dump of solution matrix/vectors of several steps to the same file. 
Changed: Moved the myInterfaces member to class SIMinput (only used there).
 2017-09-20 08:41:29 +02:00
Knut Morten Okstad
dec5d33e3f Fixed: Seems like the unordered_map container sorts after all (STL bug?) 
So replace by a vector of pairs instead.
 2017-09-16 08:18:35 +02:00
Knut Morten Okstad
5f73052452 Added: unit test for adding explicit node sets 2017-09-16 08:16:30 +02:00
Knut Morten Okstad
7ef84d2350 Changed: Undo of 02d84c07e2 
It feels unneccassary and confusing with two versions of solveStep
and the need for the using statement when one of them is overriden.
 2017-09-12 17:05:40 +02:00
Arne Morten Kvarving
5001025d2a changed: use bezier extraction in evaluation of secondary solution. 2017-09-08 11:13:08 +02:00
Arne Morten Kvarving
50bc13f17d disable multithreading for LR splines 
currently broken
 2017-09-07 14:10:03 +02:00
Knut Morten Okstad
0aff9324e7 Fixed: Copy the norm prefixes to avoid valgrind complaint 2017-09-01 09:11:42 +02:00
Knut Morten Okstad
af7d7b6728 Added: Specification of nodes in a node set on the input file. 
Fixed: Use a unordered_map instead of a map for the node sets,
to avoid that the list is sorted with respect to the names, which
would invalidate the set indices returned by the getNodeSetIdx method.
 2017-08-29 21:27:09 +02:00
Knut Morten Okstad
e19d34bbd7 Changed: utl::parseIntegers can now also read a list of integers. 
Changed: Using EXPECT instead of ASSERT in the test programs, and
put the xml input in a character string instead of a file.
 2017-08-29 21:27:09 +02:00
Knut Morten Okstad
d181ccb048 Fixed: Set nNod correctly for multi-patch models with shared grids 2017-08-29 21:27:09 +02:00
Knut Morten Okstad
ed1bd4ef5e Fixed: SIMbase::getBoundaryNodes now works also for explicit node sets 2017-08-29 21:27:08 +02:00
Knut Morten Okstad
fcf078e8a6 Added: New method ThreadGroups::oneGroup in case of no multi-threading 2017-08-29 21:26:54 +02:00
Knut Morten Okstad
d21cad9774 Added: 2D Lagrange sub-class for FE model from Matlab file 2017-08-29 21:15:11 +02:00
Knut Morten Okstad
b14cec6b63 Added: Protected method setCoord to assign nodal coordinates. 
Changed: Don't set fe.u, fe.v and fe.w if no spline parameters available,
and use iel when calculating fe.iGp such that it also works spline-less.
 2017-08-29 21:15:11 +02:00
Knut Morten Okstad
de1d029f2b Added: New virtual methods getNodeSetIdx and getNodeSet. 
Also added a method to prescribe/fix a list of nodes.
 2017-08-29 21:15:11 +02:00
Knut Morten Okstad
87c6c2ccb8 Added: Member variables for polynomial order. 
Fixed: One evalSolution signature for the Lagrange ASMs.
 2017-08-29 21:15:11 +02:00
Knut Morten Okstad
021c205a63 Fixed: Avoid adding MPC for nodes with identical node numbers. 
It resulted in infinite recursion and stack overflow.
Changed: Use range-based for loops where more appropriate.
 2017-08-29 21:08:03 +02:00
Knut Morten Okstad
d5948ec5d6 Added: Optionally collapse degenerated edges into a single point/node 2017-08-29 21:08:03 +02:00
Kjetil Andre Johannessen
f046ca76bd bumped LRSpline library for multipatch connections 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
6a85046c8c added: multipatch refinement (2D and 3D) 
fixed: we have to refine the mesh with the highest order

we have to remap the element wise errors from the geometry mesh
elements to the refinement mesh basis functions.
 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
364ef7c56b fixed: only use MNPC for first basis for injecting projected solution into integrands 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
e630986c53 changes needed to make projection work with mixed 
- only use nodes on first basis where relevant
- use extractPatchSolution instead of directly calling
  extractNodeVec to tap into additional MADOF functionality
 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
76e26e68bd changed: flag for whether or not to include other bases in additional MADOFs 
recovered quantities are only using basis1 so we need to be able to
create a MADOF which holds dummies for the other bases.

also clear addMADOFs in SIMbase::clear() to make it work for adaptive simulations
 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
88fe4728a0 fix assembleL2Matrices in ASMu2Dmx 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
ddf1511e88 implement ASMu2Dmx::evalSolution 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
dc53797f66 added: hessian support in ASMu2Dmx 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
61a84b33c3 fixed: store element size in FiniteElement for various ASMs 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
b1e6ba8cbb changed: keep all topology information 
not just those to off-process patches. needed for multi-patch
(adaptive) refinements
 2017-08-29 10:22:33 +02:00
Arne Morten Kvarving
e7db893c1f fixed: check for nullptr before dereferencing 2017-08-29 09:11:02 +02:00
Arne Morten Kvarving
47a3080a5e changed: move Interface class out of DomainDecomposition 
will be reused for intra-process connection representations
 2017-08-28 11:18:37 +02:00
Knut Morten Okstad
ce3c1e6a59 Changed: The bool argument of SIMbase::getNoNodes is removed. 
Its default value (false) was not used in practice and it was
thus only a potential source of confusion.
Fixed: Also include the 'X'-nodes (if any) when counting the nodes
associated with the first basis. That means, we assume that only
the first basis can have Dirichlet boundary conditions defined in
local coordinate axes.
 2017-08-28 09:37:10 +02:00
Arne Morten Kvarving
a8bfe42442 added: explicit warning in AdaptiveSIM::adaptMesh with negative norms 
happens if the orientation of your model is wrong
 2017-08-22 15:09:54 +02:00
Knut Morten Okstad
d346f85263 Fixed: Partly undo of commit 89faa82c6f 
Do only one iteration if the first iteration already satisfies the
convergence tolerance when refNopt=ALL. No check on subiter&FIRST here.
 2017-08-22 09:34:44 +02:00
Arne Morten Kvarving
bab3a6aa03 fixed: generate topology and topologysets even with a pre-existing model 
this happens in adaptive simulations, where the model is kept but the
rest of the input files is re-processed.
 2017-08-22 08:32:35 +02:00
Knut Morten Okstad
515f6145cc Fixed: The code inside the #if SP_DEBUG > 2 did not compile. 
B is therefore changed into a single vector (of dimension nnod*ncomp)
instead of an array of vectors of length nnod. The equation solver
can handle this format, no need to loop over the RHS outside here.
Also some indentation fixes here and there.
 2017-08-15 08:11:10 +02:00
Knut Morten Okstad
76134520b0 Added: Dump of FE model and optionally topology sets to Matlab 2017-08-15 07:38:48 +02:00
Knut Morten Okstad
7d586966ff Added: Scaling option for the non-uniform explicit refinement 2017-08-14 22:24:23 +02:00
Arne Morten Kvarving
f2080cdf54 make single-parameter constructors as explicit 
to avoid unexpected implicit casts
 2017-08-14 21:27:48 +02:00
Arne Morten Kvarving
3722992673 added: unstructured patch support in domain decomposition 2017-08-14 13:15:57 +02:00