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

2346 Commits

Author SHA1 Message Date
Arne Morten Kvarving
dbc2202b5a ASMu3D::evaluateBasis: pass vector instead of finite element and basis index 2023-09-05 11:51:08 +02:00
Arne Morten Kvarving
82ccdc49e3 ASMu3D::constrainFace 
use midpoint()
 2023-09-05 11:51:08 +02:00
Arne Morten Kvarving
3ad0b872c8 changed: rename ASMmxBase::geoBasis to ASMmxBase::itgBasis 
more precise as this represents the basis holding the integration elements
 2023-09-05 10:12:26 +02:00
Arne Morten Kvarving
6aeefcd0f0 changed: use PETSC_NULLPTR 
for PETSc < 3.19 we define this to PETSC_NULL.
use of PETSC_NULL causes a deprecation warning in 3.19
 2023-09-04 20:02:53 +02:00
Arne Morten Kvarving
9014cc6335 changed: use INSTANTIATE_TEST_SUITE_P 
INSTANTIATE_TEST_CASE_P causes deprecation warnings with
newer gtest versions
 2023-09-04 20:02:53 +02:00
Arne Morten Kvarving
ab100c290d changed: use getBasis in ASMxxD::assembleL2matrices 
make it clear that we want the integration and projection bases.
 2023-09-01 17:20:42 +02:00
Arne Morten Kvarving
ad7ff2b19e ASMuxDmx: use getBasis(ASM::REFINEMENT_BASIS) where appropriate 2023-09-01 17:20:42 +02:00
Arne Morten Kvarving
d9056a92b2 ASMuxD: use getBasis(ASM::PROJECTION_BASIS) where appropriate 2023-09-01 17:20:42 +02:00
Arne Morten Kvarving
00486b30e8 use getBasis(ASM::PROJECTION) in ASMuxD::refine 
and restructure a little to make code more readable.
 2023-09-01 17:20:42 +02:00
Arne Morten Kvarving
c3bb2447cd use getBasis(ASM::INTEGRATION_BASIS) in ASMsxDmx::generateFEMTopology() 
this makes it a bit clearer what is going on.
introduce a lambda to avoid duplicated code and some generic cleanup
 2023-09-01 17:20:42 +02:00
Arne Morten Kvarving
f91a62e372 add tests for unstructured ASM write 
also serves as tests for the getBasis implementation
 2023-09-01 15:41:44 +02:00
Arne Morten Kvarving
9cf613cbca changed: move ASMuSquare and ASMuCube to separate header 
for reuse like in the structured ASMs tests
 2023-09-01 15:41:44 +02:00
Arne Morten Kvarving
76a6139bfc changed: no reason to override write in ASMuxDmx 
we can just use the ASMuxD implementation
 2023-09-01 15:41:44 +02:00
Arne Morten Kvarving
d57ebd0d0c add tests for structured ASM write 
also serves as tests for the getBasis implementation
 2023-09-01 14:29:10 +02:00
Arne Morten Kvarving
dc7a2e444a changed: unify ASM(Cube|Square) and ASMmx(Cube|Square) 
in particular do not call generateFEMTopology in the latter's
constructor.
 2023-09-01 14:29:10 +02:00
Arne Morten Kvarving
53f5d1f830 changed: no reason to override write in ASMsxDmx 
we can just use the ASMsxD implementation
 2023-09-01 14:29:10 +02:00
Arne Morten Kvarving
005f8f24b1 ASMuxDmx: restructure basis establishment 
make it (almost) similar to code in structured ASMs
 2023-09-01 14:09:00 +02:00
Arne Morten Kvarving
809a650bc4 fixed: avoid cloning in ASMsxDmx 
instead point to the shared ptrs. make sure to null ptrs
before ASMstruct wants to delete them.
 2023-09-01 14:09:00 +02:00
Arne Morten Kvarving
fdd7b189b1 use a lambda in ASMu3D::getBoundary1Nodes 
switch is preferable to ifs
 2023-09-01 13:46:14 +02:00
Arne Morten Kvarving
3d6d950474 changed: remove ASMxxD::get(Surface|Volume) 
use getBasis instead
 2023-09-01 08:51:04 +02:00
Arne Morten Kvarving
957af79433 added: extend getBasis method 
allow obtaining additional bases, not just the FE bases
in particular geometry basis and projection basis.
add an enumeration to make things readable
 2023-09-01 08:51:04 +02:00
Arne Morten Kvarving
4abf4c64a7 changed: override separateProjectionBasis in ASMxxD(mx) 
in preparation for having a separate geometry basis. in that case
we cannot use the geomB pointer to check for equivalence.
also fixes two field mixed formulations with the higher
order basis being the second basis. in that case we are
not projecting on basis 1, but rather on basis 2.
 2023-08-29 14:36:03 +02:00
Arne Morten Kvarving
0c2e160d90 rename some variables to make code more readable 2023-08-29 14:34:33 +02:00
Arne Morten Kvarving
9d133d8be1 move refinement basis into ASMLRSpline 2023-08-29 14:05:20 +02:00
Arne Morten Kvarving
b3b3e1fa7a ASMu2D: make createLRNurbs return a shared_ptr directly 
no reason to use bare pointers here, they are always put into
a shared_ptr
 2023-08-29 14:05:20 +02:00
Arne Morten Kvarving
59ad91c245 rename altProjThreadGroups to proj2ThreadGroups 2023-08-29 14:05:20 +02:00
Arne Morten Kvarving
505d9d9a6c move projection bases into ASMLRSpline 2023-08-29 14:05:20 +02:00
Arne Morten Kvarving
dce989f1c6 changed: remove proj members in ASMs2D and ASMs3D 
as these are stored as projB in ASMstruct.
less confusion as there are less member pointers.
 2023-08-29 10:52:23 +02:00
Arne Morten Kvarving
086592c2f9 changed: move secondary projection basis into ASMstruct 2023-08-29 10:52:23 +02:00
Arne Morten Kvarving
1d32025671 ASMuxDmx: simplify remapErrors by using elm->midpoint() 2023-08-29 10:12:45 +02:00
Arne Morten Kvarving
bd5609035b changed: rename ASMLRSpline::geo to geomB 
to be consistent with ASMstruct. we store it
as a shared_ptr. too many pointers to keep manually
track of. it can be (soon) be a separate object,
one of the FE bases, the projection basis, ...
 2023-08-29 09:11:09 +02:00
Knut Morten Okstad
3a4fc9e081 Changed: Second (mostly dummy) argument in read method is removed. 
Instead, introducing a separate readBasis method.
 2023-08-25 16:43:46 +02:00
Arne Morten Kvarving
0a33cedfc3 changed: introduce a method to obtain first node on integration basis 
use this in ASMsxD::getElementCoordinates. we can then
reuse this method in ASMsxDmx
 2023-08-25 12:26:26 +02:00
Arne Morten Kvarving
c0b5911549 changed: introduce a method to obtain last node on integration basis 
use this in ASMsxD::getParametric[Area|Length|Volume]. we can then
reuse these method in ASMsxDmx
 2023-08-25 12:26:26 +02:00
Arne Morten Kvarving
deb2bd53ea fixed: make sure to always set geoBasis 
if not we risk getting the wrong setting if we instance
a SIM with e.g. subgrid followed by instancing a SIM with Taylor-Hood.
 2023-08-25 11:16:05 +02:00
Knut Morten Okstad
bc066b9336 Changed: Unnecessary copying of bfs.N removed. 
Fixed: Debug print in ASMu2D::integrate().
 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
f2f876c77c added: BasisFunctionCache in ASMu3Dmx 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
730fdca306 added: BasisFunctionCache in ASMu3D 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
2f676bf6b3 added: BasisFunctionCache in ASMu2Dmx 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
f2924c1cbd added: BasisFunctionCache in ASMu2D 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
bffcfa640f added: BasisFunctionCache in ASMs3DmxLag 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
e595bd3556 added: BasisFunctionCache in ASMs3Dmx 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
003d567cd3 added: BasisFunctionCache in ASMs3DLag 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
a247271241 added: BasisFunctionCache in ASMs3D 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
a6debc90af added: BasisFunctionCache in ASMs2DmxLag 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
f268e533a8 added: BasisFunctionCache in ASMs2Dmx 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
eff2b08e03 added: BasisFunctionCache in ASMs2DLag 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
6ce589a6b5 added: BasisFunctionCache in ASMs2D 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
c168894ebf added: input file option for defining cache policy 2023-08-15 11:47:15 +02:00
Arne Morten Kvarving
44bbec509a added: base class for basis function cache 2023-08-15 11:47:15 +02:00
Knut Morten Okstad
4fc972737a Fixed: G_MATRIX and AVERAGE calculation for 3D mixed LR 2023-08-01 14:42:06 +02:00
Knut Morten Okstad
58a1bf78d1 Added: Error exit if integrand with reduced integration for mixed LR. 
Fixed: Solution transfer for mixed in 3D.
 2023-08-01 12:23:22 +02:00
Knut Morten Okstad
d8ac636fee Added: Convenience method in utl namespace for calculation of (1st and 2nd) 
derivatives of mixed basis functions w.r.t. Cartesian coordinates,
including the Jacobian and Hessian matrices
 2023-08-01 10:26:02 +02:00
Knut Morten Okstad
91b9a07230 Removed unused and potentially dangerous evaluateBasis method, 
storing basis function derivatives w.r.t. parameter domain coordinates
in the array intended for derivatives w.r.t. Cartesian coordinates
 2023-08-01 10:26:02 +02:00
Knut Morten Okstad
e63aebcc03 Added: Convenience typedef SplinePtr in ASMu2Dmx and ASMu3Dmx 
and use it instead of auto in range-based loops for readability.
Changed: Use getElement(iel) instead of elementBegin()+iel.
Fixed: Don't increment element counter when skipping element
not in current threading group when doing interface integral.
 2023-08-01 10:26:02 +02:00
Knut Morten Okstad
3d3816130c Added: Convenience method getElementsAt() for ASMu2Dmx and ASMu3Dmx. 2023-08-01 10:26:02 +02:00
Knut Morten Okstad
78635f56d3 Fixed: Integration point counter for mixed integrand. 
Changed: Return false instead of break in loop without threading.
 2023-08-01 10:26:02 +02:00
Arne Morten Kvarving
733c9c51a9 added: set parameters in SplineUtils::project 2023-06-08 11:00:43 +02:00
Arne Morten Kvarving
29d70c1262 fixed: initial conditions with partitioning 
calculate these on root node and broadcast instead of calculating
separately on each process. this can cause desync between processes,
for instance if the initial condition uses random numbers
 2023-06-08 10:24:58 +02:00
Arne Morten Kvarving
6736f7829f added: broadcast wrapper to ProcessAdm 2023-06-08 10:24:58 +02:00
Arne Morten Kvarving
933bc34944 ASMu2Dmx: add rational support 2023-05-26 09:40:42 +02:00
Arne Morten Kvarving
6a227741d6 added: NURBS support in LRSplineField(s)2D(mx) 2023-05-25 08:36:08 +02:00
Arne Morten Kvarving
84badbb8c8 changed: make NURBS basis evaluation methods static and public 
for reuse in fields classes
 2023-05-25 08:36:08 +02:00
Arne Morten Kvarving
be6dd45918 changed: make getCoordinates a static and public method 
for reuse in field classes
 2023-05-25 08:36:08 +02:00
Arne Morten Kvarving
046b242920 added: function to query whether a ASMu2D is rational or not 2023-05-25 08:36:08 +02:00
Arne Morten Kvarving
835384655c fixed: handle rational projection basis 2023-05-23 20:47:09 +02:00
Arne Morten Kvarving
373cdf281e changed: convert lambda to static member function 2023-05-23 20:47:09 +02:00
Arne Morten Kvarving
ada3102a76 ASMu2Dmx: use computeBasis wrapper 
in preparation for rational support
 2023-05-23 13:09:28 +02:00
Arne Morten Kvarving
a107664301 added: eps output functions for LR nurbs 2023-05-23 12:03:40 +02:00
Arne Morten Kvarving
acb0c471eb changed: unify computeBasis signatures 
always take a spline pointer parameter. use reference in the NURBS interfaces.
 2023-05-23 11:38:54 +02:00
Arne Morten Kvarving
403cc925bc fixed: use internal basis evaluation in jump integration method 
to handle NURBS
 2023-05-23 11:13:58 +02:00
Knut Morten Okstad
f9c1be46ed Fixed a valgrind issue in test problem (missing initialization) 2023-05-23 09:30:15 +02:00
Arne Morten Kvarving
9ec1b0d6d9 fixed: have to regenerate bezier basis when switching projection basis 
we also have to switch geometry as we need elements without the
additional refinement when projecting the pressure
 2023-05-22 14:16:35 +02:00
Arne Morten Kvarving
bb7d87efb9 changed: fold ASMu2Dnurbs into ASMu2D 
this so we get access to NURBS in ASMu2Dmx.
remove ASM::LRNurbs and option -LRn as we now detect
rational from the underlying tensor spline
 2023-05-16 22:02:36 +02:00
Knut Morten Okstad
81b341d03f Fixed: External forces where added twice for HHTSIM 2023-05-03 06:12:02 +02:00
Knut Morten Okstad
d775740068 Added: Methods for debug print of element matrices, with optional names 2023-05-03 06:11:18 +02:00
Knut Morten Okstad
234f4c152b Added: Unit test for NewmarkNLSIM 2023-05-02 13:20:47 +02:00
Arne Morten Kvarving
fb34513b76 added: set parameters in ASMs3DLag::integrate(B) 2023-03-27 14:07:55 +02:00
Knut Morten Okstad
6f0d8d56a6 Changed: ASMbase::globalL2projection() does not need to be virtual. 
Added: Enforcement of end point values in the global L2-projection.
 2023-02-15 07:38:04 +01:00
Knut Morten Okstad
d2c496c920 Fixed: Use actual spline parameters when evaluating functions - not the normalized ones. 
Changed: Use matrix::fillColumn() instead of loop.
 2023-02-14 13:37:18 +01:00
Knut Morten Okstad
9e35a86598 Cosmetic doxy fixes 2023-02-14 12:07:08 +01:00
Knut Morten Okstad
747b21b655 Changed: Simply VariationDiminishingSplineApproximation() signature 2023-01-10 13:47:06 +01:00
Knut Morten Okstad
90f7129729 Changed: Use && instead of bit-wise &, to avoid warnings with newer compilers. 
The result is that the mesh will not be saved if the refinement fails,
which probably makes more sense anyway.
 2023-01-10 13:39:52 +01:00
Knut Morten Okstad
d89c31558f Changed: Return int instead of bool such the bitwise or 
can be used on the result without warning for newer compilers
 2023-01-10 13:39:52 +01:00
Arne Morten Kvarving
4ea903be4f Fixed: Wrong return value 2023-01-10 13:39:52 +01:00
Knut Morten Okstad
4618eab9d0 Added: Optional element shrinkage factor for LR meshes to VTF 2022-10-29 11:50:45 +02:00
Knut Morten Okstad
69e43205ca Fixed: Potential infinite recursion 2022-10-04 14:54:14 +02:00
Arne Morten Kvarving
a8194b9500 changed: avoid use of PetscPrintf 
the format specifier %D is deprecated and non-working in newer
petsc. rather than figuring out how to do this properly across
petsc versions, use our own printing facilities.
 2022-10-04 13:45:02 +02:00
Arne Morten Kvarving
080fe1c2bc fixed: graph based partitioning with MPC constraints 2022-10-04 13:45:02 +02:00
Arne Morten Kvarving
bebf8c1313 added: SAM::getUniqueEqns 
this returns the unique equations for a given element.
in particular it resolves MPC couplings
 2022-10-04 13:45:02 +02:00
Knut Morten Okstad
6d42ff4e22 Added: AlgEqSystem methods addScalar() and getScalars(), and the 
SIMbase method extractScalars(). The latter is a generalization and
replacement of the virtual method getExtLoad() which is no longer needed.
 2022-10-04 10:13:19 +02:00
Knut Morten Okstad
f9370d3e67 Fixed: isConstant() for some function types 2022-10-04 10:13:19 +02:00
Knut Morten Okstad
e2c76eaa08 Added: Time-derivative of time-dependent spatial functions 2022-10-04 10:13:19 +02:00
Arne Morten Kvarving
cfa6a8b4a2 fixed: periodic single patch models with multiple bases 2022-10-04 09:51:05 +02:00
Arne Morten Kvarving
d748bdc0eb added: utl::getDirs 
this creates a direction integer like '123'
 2022-10-04 09:51:05 +02:00
Arne Morten Kvarving
0c66ef79ee changed: facilitate reuse of ForceBase derived integrands 
previously there was no way to clear the buffers used during
integration. now initBuffer resizes only if necessary, but
always zeros the underlying buffer
 2022-08-01 13:45:57 +02:00
Arne Morten Kvarving
ad4264505e use INSTANTIATE_TEST_SUITE_P instead of INSTANTIATE_TEST_CASE_P 
latter is deprecated
 2022-05-19 11:26:24 +02:00
Arne Morten Kvarving
063fb58b26 added: proper preallocation for block matrices with graph partitioning 2022-05-19 10:59:43 +02:00
Arne Morten Kvarving
042a70bf65 fixed: fix SplineUtils::project for 3D functions 2022-05-13 11:58:09 +02:00
Arne Morten Kvarving
bbb90d468e changed: suppress profiler output in unit tests 2022-05-11 15:16:13 +02:00
Arne Morten Kvarving
a4acd0efac added: support dune-istl >= 2.7 2022-05-11 10:52:02 +02:00
Arne Morten Kvarving
6bbe09a312 fixed: add missing initializers 2022-05-10 11:32:16 +02:00
Arne Morten Kvarving
b2707a8137 changed: make SIMbase::project virtual 
needs to allow overriding in subclasses for when
implementing support for alternative projection bases
 2022-05-10 11:32:16 +02:00
Arne Morten Kvarving
459463878a introduce support for an alternative projection basis 2022-05-10 11:32:16 +02:00
Arne Morten Kvarving
1d88adf5af added: getProjectedField 
this returns a scalar field over the projection basis
with the given coefficients.
 2022-05-10 11:32:16 +02:00
Arne Morten Kvarving
fa73096ebe changed: remove support for subgridH 
no longer needed, pressure and velocity are
projected on separate bases now
 2022-05-10 10:15:10 +02:00
Arne Morten Kvarving
73d6e85745 fixed: wrong projection basis for tensor subgrid 2022-05-10 10:15:10 +02:00
Arne Morten Kvarving
3e6edab293 added: allow instancing a FieldFunction from an existing field 2022-05-10 09:59:38 +02:00
Arne Morten Kvarving
f27b3a4efb added: allow instancing scalar fields directly from a surface/volume 2022-05-10 09:59:38 +02:00
Arne Morten Kvarving
6b861e4789 added: initPatch member in Integrand 
this tells the integrand the index of the patch currently
being processed.
 2022-05-10 09:06:02 +02:00
Arne Morten Kvarving
0fb14fe00d fixed: handle 'default' preconditioner for ISTL 2022-05-09 11:02:37 +02:00
Arne Morten Kvarving
5391cfbc6b added: support for block matrices for graph partitioned models 2022-05-09 11:02:37 +02:00
Arne Morten Kvarving
646bf62b0f changed: refactor setup of sparsity pattern into separate methods 
makes the code more readable
 2022-05-09 11:02:37 +02:00
Knut Morten Okstad
048e37ffcb Fixup for commit 3ce4b7f36904303a01f0af1ef61ac5f5e402024a: valgrind complaint 2022-04-28 08:45:57 +02:00
Knut Morten Okstad
6991c28f31 Added: Extra (optional) int argument in writeGlvA 2022-04-26 16:25:20 +02:00
Knut Morten Okstad
5d7270e612 Added: Internal forces debug print 2022-04-26 16:25:20 +02:00
Knut Morten Okstad
3ce4b7f369 Fixed: Clear reaction force component after it being added when a group 
of nodes is specified, such that nodes that are member in more than one
set are not added twice
 2022-04-26 16:25:20 +02:00
Knut Morten Okstad
b672585ecc Changed: Also log multi-threading on faces to IFEM::cout 2022-04-26 16:25:20 +02:00
Knut Morten Okstad
6b322da4b6 Changed: Generalize ReactionsOnly class for internal force calculation. 
Added: SIMbase::assembleForces() and IntegrandBase::setSecondaryInt().
Added: virtual method for calculation/print of interface force resultants.
 2022-04-26 16:25:20 +02:00
Arne Morten Kvarving
37d33f0648 fixed: set parameters when projecting function 
necessary when using a fieldfunction
 2022-04-26 12:45:04 +02:00
Arne Morten Kvarving
e3d3b526e8 fixed: clear additional functions in clearProperties 
this is necessary to not register multiple functions if
used in adaptive simulations
 2022-04-26 11:39:18 +02:00
Arne Morten Kvarving
cadcefb454 fixed: memory leak in mixed LR integrands 
tensorspline was not cleared up after LR bases
had been established
 2022-04-26 11:37:34 +02:00
Arne Morten Kvarving
ddf260795f changed: add 'symmetric' as an alias for 'symmetrized' beta type 2022-04-22 12:15:39 +02:00
Arne Morten Kvarving
71fa5975da fixed: check number of components in function 
in particular for mixed integrands using anasol boundary conditions
nf can be greater than the number of components in the function. this is
because nf is set to the total number of fields and not the number
of fields on a single basis
 2022-04-22 10:57:33 +02:00
Arne Morten Kvarving
e80aef161b re-enable multithreaded mixed 3D assembly 
the problems leading to it being disabled should have been fixed now
 2022-04-20 12:20:10 +02:00
Arne Morten Kvarving
673eec21c7 added: more field function types 
gradient and laplacian of scalar and vectorial fields.
this will be used for deriving source functions for simulations
where we use another simulation result as "exact" solution
 2022-04-20 12:19:20 +02:00
Arne Morten Kvarving
02de121ca3 added: ability to use a fixed time level in FieldFunctions 2022-04-20 12:19:20 +02:00
Arne Morten Kvarving
c5c43beca2 fixed/add: hessianFE in SplineFields 
previously the derivatives on the reference domain was returned
 2022-04-20 12:19:20 +02:00
Arne Morten Kvarving
9bfe8af84e added: support FieldFunctions for mixed spline fields 
If no dash is found in the basis name, we assume the field
names given are on separate bases, one per basis

this is used to instance fieldfunctions for div-compatible
solution fields
 2022-04-20 12:19:20 +02:00
Arne Morten Kvarving
3dc01f07a4 added: implement valueCoor in mixed spline fields 2022-04-20 12:19:20 +02:00
Arne Morten Kvarving
e479951031 added: ability to read mixed bases from streams 2022-04-20 12:19:20 +02:00
Knut Morten Okstad
421fbdf198 Changed: Enable MathJax for equation rendering 2022-03-08 15:31:27 +01:00
Knut Morten Okstad
d015090e9f Fixed: Calculation of reaction and interface forces. 
Boundary contributions due to in-homogeneous neumann conditions were missing
due to the flawed check on the presence of myEqSys->getVector().
We now rely on the hasBoundaryTerms() method of the integrand instead.
Use property type OTHER when internal forces are calculated,
such that dirichlet boundaries in the same patch that does not belong
to the interface we are calculating interface forces on are skipped.
 2022-03-02 12:30:16 +01:00
Knut Morten Okstad
58e5969947 Added: Convenience method SAM::printVector (for debugging) 2022-03-02 12:30:16 +01:00
Knut Morten Okstad
7736de71ca Added: SIMgeneric::getInterfaceForces() 2022-03-02 12:30:16 +01:00
Knut Morten Okstad
88ee21ffe0 Changed: Use SIMbase::getBoundaryNodes in SIM::initBoundaryNodeMap() 2022-03-02 12:30:16 +01:00
Knut Morten Okstad
d984430fca Added: Extract internal forces from ReactionsOnly integral 2022-03-02 12:30:16 +01:00
Knut Morten Okstad
a5fd95b0ed Fixed: Don't flag Neumann-order or local tangent direction 
for edge boundaries in 3D through the Property::lindx member.
It is not needed and also makes it easier to address volume
edges (in range 1 to 12) correctly without the mod'ing.
 2022-03-02 11:53:07 +01:00
Knut Morten Okstad
83ae2c73fa Fixed: Ensure unique boundary nodes nodes also for ldim=nsd and ldim=4. 
Added: Also consider 1D boundaries (edges) in 3D and vertices in 2D.
 2022-03-02 11:53:07 +01:00
Knut Morten Okstad
3712524a0e Replaced virtual method getEdge() by getBoundary1Nodes() 
with a similar interface as getBoundaryNodes()
 2022-03-02 11:53:07 +01:00
Knut Morten Okstad
5c3e3d2bb2 Fixed: The *sol vector is considered the control point values of the 
projected secondary solution if the description field equals "projected",
and the results field contains DataExporter::SECONDARY.
 2022-03-02 11:42:20 +01:00
Arne Morten Kvarving
7787ed15d9 fixed: these are pointers 
use nullptr, not 0
 2022-02-17 09:19:25 +01:00
Arne Morten Kvarving
2e6002d024 fixed: we have to set the option here to avoid trouble 
if we have previously reused the matrix factorization,
the FACTORED value is stored. this leads to segfaults.
we now properly flag that we only want to reuse the
sparsity pattern
 2022-02-17 09:19:25 +01:00
Arne Morten Kvarving
ab6d03b141 added: extend element matrix cache support to mixed 2022-02-17 08:45:06 +01:00
Arne Morten Kvarving
a3555e8942 changed: do not resize the newton matrix up front 
this way we can use its non-emptyness to signal that
it has already been assembled.
 2022-02-16 22:44:10 +01:00
Arne Morten Kvarving
5618e4b000 changed: print iteration history in NonLinSIM even at msgLevel -1 
allows suppressing norm and solution summary output, while still
keeping the iteration history. useful in split simulators
such as Chorin
 2022-02-14 07:00:47 +01:00
Arne Morten Kvarving
ffe5499917 fixed: wrong variable used in 3D mixed tensor spline fields 2022-02-11 15:32:38 +01:00
Arne Morten Kvarving
22d7167c7d fixed: we need to free the MPI_Comms even in serial 2022-02-09 12:52:53 +01:00
Arne Morten Kvarving
ec6b692c3a fixed: also put multiplier in block 2 with component based blocks 2022-01-27 11:44:00 +01:00
Arne Morten Kvarving
a2cf0932f1 Print number of local equations in each block 2022-01-27 11:44:00 +01:00