Skip to content

Latest commit

 

History

History
662 lines (574 loc) · 31.3 KB

CHANGELOG.md

File metadata and controls

662 lines (574 loc) · 31.3 KB
Raw

Changelog

All notable changes to this project will be documented in this file.

4.14

Added

  • Finite element BDM2 and BDM2ortho in test, Bug in BDM2ortho corrected the 4 sept 2014 in version: v4.13-130-g1af52457
  • Conversion of matrix or transpose of matrix in int[int][int] array to get the structure of sparse matrix. see tutorial/sparse-matrix.edp example at end 
    matrix A = va(Ph,Vh); 
int[int] a = A, at= A';
  • a meshL finite function can be see as real function with 1, or 2 parameters 
     meshL ThL = segment(10); fespace VhL(ThL,P1); VhL u= x;
 cout << u(0.5)   << endl; 
 cout << u(0.5,0) << endl;
  • Exemple to code convolution of 2 function with one with a small support too be not to expanxive see tutorial/Convolution-Sample.edp example
  • Support for dense blocks in PETSc matrices
  • GenEO for saddle-point examples with PCHPDDM in PETSc
  • Distributed ParaView output on meshS
  • Interface to mmg2d for two-dimensional mesh
  • Support for Mmg parameters localParameter, -nosizreq, -hgradreq`

Changed

  • PETSc 3.20.2

Deprecated

Removed

Fixed

  • bug in P3pnc3d in vectorial case (thank to [email protected] )
  • in segment(10,region=1,label=ll); region is now used..

4.13

Added

  • Composite FE spaces and variational forms for coupled problems:

    • can now define composite FE spaces with different meshes/mesh types as 
      fespace Uh(Th1,[P2,P2]);
fespace Ph(Th2,P1);
fespace Vh=Uh*Ph;
    • can define coupled problems using composite FE spaces, or directly with < > syntax: 
      fespace Uh(Th1,[P2,P2]);
fespace Ph(Th2,P1);
Uh [u1,u2],[v1,v2];
Ph p,q;

solve Pb (<[u1,u2],[p]>, <[v1,v2],[q]>) = ...
      see examples/examples/stokes_composite.edp and examples/examples/stokes_periodic_composite.edp
    • this new type of composite problem can be used for FEM-BEM coupling and also benefits from automatic parallel assembly (in test) and can be easily solved using the distributed solver MUMPS, see examples/bem/Helmholtz-2d-FEM-BEM-coupling-MUMPS-composite.edp
    • composite problems can also be solved using PETSc (in test), see examples/hpddm/Helmholtz-2d-FEM-BEM-coupling-PETSc-composite.edp

  • remove spurious cout in Curve/Line DG definition.

  • add New Finite element 2d on mesh : RT0dc (discontinios RT0 ) in plugin Element_Mixte see example plugin/RT0dc.edp and P1nc (Crouziex-Raviart) + bulle : name P1bnc in plugin Element_P1ncdc and P1nc totaly discontinous + bulle ; name P1bdcnc in plugin Element_P1ncdc see example plugin/example testp1dcnc.edp for [email protected]

  • add New finite element: P4S P^4 on meshS , P3pnc3d in Element_P3pnc_3d (Couziex-Raviart with P3 ) see [email protected] for more information.

  • add new interface for metis (see examples/plugin/metis.edp)

  • Correct jump, mean, otherside of finite element function on mesh3, meshS, meshL (add missing code in method: MeshPoint::SetAdj() thanks to [email protected])

  • try to build dmg install mac version

  • add file script to build meshS from boundary meshL TL if the boundary is the graph of function from mean plane. see example in examples/3dSurf/buildmeshS.edp meshS Ts=buildmeshSminsurf(TL,1);// minimal surface meshS Tsl=buildmeshSLap(TL,1);// Laplace Surface .. meshS Tsl=buildmesh(TL,1,op);// op = 0 Lap and op =1 => minsurf.

  • add sparse block to sparse matrix matrix A = va(Vh,Vh); matrix B(A.n5,A.n5); int i=2; B.add(1.+10i,A,indof,i*ndof);

Changed

  • change isoline to do the job for meshS, see example plugin/isoline.edp
  • change Curve function to be with 3 components to use the isoline data.
  • change Curvature plugin to compatible with new isoline data for 3 d case.
  • change some sprintf in snprint to remove warning

Deprecated

Removed

Fixed

  • bug in all P0face, P0edge, P0VF on mesh3,meshS, MeshL and also discontinous version (missing initialisation)
  • bug in plot function and ffglut with parameter pdf="file.pdf" , because shift in plot named parameter not change in ffglut.
  • genere a bug if zero size element in read MeshL from file.
  • remove mistake when the border is badly defined , remove empty element in buildmeshL function.
  • bug in array quadrature FE.

4.12

Added

  • add new finite Element P2pnc3d of Stokes problem like Crouzeix-Raviard in 3d of P2 pylynome (see G. Allaire or [email protected] for detail )

  • add pdfPLOT form [email protected] [http://www-an.acs.i.kyoto-u.ac.jp/~fujiwara/ff++-programs/] usage: plot( ..., pdf="filename.pdf", svg="filename.svg" );

  • add missing code for Discontinous Galerkin in 3d for RHS ( see [problem-in-3d-discontinuous-galerkin-computation][https://community.freefem.org/t/problem-in-3d-discontinuous-galerkin-computation/2015/6] )

  • add in examples/mpi/chamonix.edp : radiative transfer (use new plugin plugin/mpi/RadiativeTransfer_htool.cpp, illustrates the use of htool for compression of user defined matrix operator)

  • transform a surface meshS in 2d mesh (warning with overloaping, no test) with movemesh:

    meshS Ths = square3(10,10,[x,y,square(2x-1)+square(2y-1)]); real[int] gzz; mesh Th2 = movemesh(Ths,transfo=[x,y,z],getZ=gzz);// get flat 2d mesh form meshS

  • New 1d finite element P3 hermite (C1) finite element in plugin Element_P3 meshL Th=segment(1,[x*L,0,0]); fespace Vh(Th,P3HL); see example end of example plugin/testFE-P3

  • missing new 1d finite element P4 in plugin Elemnt_P4

  • plugin plugin/seq/MatrixMarket.cpp to read and save matrix in MatrixMarket and add also a binary form (see examples/plugin/MatrixMarket.edp test)

  • add ILU on complex matrix in plugin IncompleteCholesky remark:the IncompleteCholesky is writen but not tested

  • add test of functionnal interface of complex eigen value problem in examples/eigen/LapEigenValueFuncComplex.edp

Changed

  • correct some old code with old version of K.facePermutation() function in plugin/seq/Element_Mixte3d.cpp and plugin/seq/Element_P2bulle3.cpp (not tested)

Deprecated

Removed

Fixed

  • fixe in A.RemoveHalf (alway return a new matrix)

4.11

Added

  • add computation scalar product of R3 example : ( N'*Tl)

  • add tools to do compution with R3 vector see tutorial/calculus.edp

  • add an example tutorial/tgv-test.edp see see what tgv do on matrix build.

  • add R3 Th.be(k).N to get the normal of boundary element (in all mesh type)

  • add R3 Th.be(k)[i].P to get the point (R3) of boundary vertices

  • add R3 Th.be(k).measure to get the measure of the boundary elment

  • add projection function to a mesh , meshL, MeshS or mesh3 with return a R3 point

  • see new example dist-projection.edp example in examples

  • add dxx, dyy, dzz, dxy, .. on P2L finite element

  • add tools to compute solid angle : let R3 O; a given point, Th3 a mesh3 and ThS a meshS. solidangle(O,Th3.be(ke)) // triangular face is the boundary face solidangle(O,Th3[k],nuface) // triangular face is face nuface of tet Th3[k] solidangle(O,ThS[k]) // triangular face is ThS[k] solidangle(O,A,B,C) // triangular face i (A,B,C) Volume(O,Th3.be(ke)) // O, triangular face is the boundary face Volume(O,Th3[k],nuface) // O, triangular face is face nuface of tet Th3[k] Volume(O,ThS[k]) // O, triangular face is ThS[k] Volume(O,A,B,C) // (O,A,B,C) tet ..

  • in bem pluging add array of HMatrix

  • examples/3d/Connectivite-3d.edp or /3dSurf/Connectivite-S.edp of test.

  • 3 function mapk, mapkk, mapkk to set a function in fourier space with k parametre R3 K; // le fourier variable allway 3d (sorry) int n1=16,n2=8, n3=4; real[int] tab1(nx,tab2(nxny),tab3(nxny*nz); mapk(tab1,K,sqr(K.x)); mapkk(tab2,ny,K,K.norm2); mapkkk(tab3,ny,nz,K,K.norm2); // Remark you can change K by P (current point)

  • in SurfaceMesh.ipd fonction to build a Isocaedron and a Sphere from this Isocaedron

  • new finite element on MeshS this finite element is the ortogonal of RT0 on surface, or Nelelec Finite Element on triangle with one DoF per mesh edge and where the DoF is the current on Edge in orientate edge by number of vertices.

  • plugin Element_P3pnc for new 2d finite element P3pnc (P3 + 2 bulles) nonconforming (continuite of P2 mod)
    and add 2 examples with this new finite element examples/plugin/cavityNewtowP3pnc.edp examples/plugin/testFE-P3pnc.edp

  • plugin Element_P3nc for new 2d finite element P3pnc (P3 ) nonconforming (continuite of P2 mod)
    and add examples with this new finite element examples/plugin/testFE-P3pnc.edp

  • function to set dirichlet Boundary conditon on matrix A (real ou compex) trought an real[int] (if none zero => set BC ) setBC(A,au1[],-2); and the example examples/3d/Elasticity-simple-support-BC.edp

Changed

  • the beaviour of linear solver UMFPACK, CHOLMOD in case of error , now FreeFEm exit on ExecError like in MUMPS
  • PETSc 3.17.0

Deprecated

Removed

  • map function in plugin dfft

Fixed

  • pow(int,int) now call int version not complex version..
  • correct the normal the N implicite variable on meshL case
  • correct version dump in banner FreeFem++ - version 4.10 (V ...
  • correct in CPU time on big mesh due to do bad HCode in HashTable.hpp
  • bug in array of finite element on meshhS, meshL (ie. fespace Vh(ThS,[P1,P1]); )

4.10

Added

  • ridgeangle named parameter in ExtractMeshL in msh3 plugin
  • DG formulation in 1d : add integral of all border of element : intallBE(ThL) and unified the notation by adding intallBE(ThS) , intallBE(Th2), intallBE(Th3) nuVertex of now the vertex number of element in intallBE0d integral BoundaryBE, InternalBE to know if border element (BE) is on true boundary of not. update nElementonB in case on no manifold data (value greater > 2) in meshL, MeshS case .. add code to use jump, mean, otherside of test function on MeshL case. ( not in mesh3 ) to compute RHS.
  • add getcwd() function in shell plugin to get the current working dir
  • add nuVertex to get the vextex on element in some int?

Changed

  • PETSc 3.16.1

Deprecated

  • SLEPc and SLEPc-complex have been part of PETSc and PETSc-complex for multiple releases and are now deprecated

Removed

Fixed

  • examples/potential.edp correct problem in times loops and BC
  • tutorial/mortar-DN-4.edp correct problem of region number in meshL
  • fix problem in Curve mesh and intallBE , vertex number is wrong
  • portability issue on arm64-apple with make petsc-slepc
  • fix assertion failure with transfer and transferMat with some finite elements

4.9

Added

  • add P3 lagrange finite element on meshS and meshS
  • add new plugin meshtoolto add tool to compute the number of connected components of a all kind of mesh (mesh,mesh3,meshS,meshL) with 2 kind of connected components ones on interior part of the mesh (default) ans secondly on the closure of the mesh (see examples/hpddm/bConnectedComponents.edp ) add functions int[int] In=iminP1K(Th,u) or int[int] Ix=imaxP1K(Th,u) get the array min/max of value u[i]
    where i is vertex number on each element k, so we have u[Im[k]] = min u[i]/ i in k;
  • add in plugin bfstream to to read binary int (4 bytes) to read fortran file and try to pull tools to share the endiannes in progress
  • add gluemesh of array of MeshL and MeshS type
  • interface to PC_MG_GALERKIN_BOTH
  • Kronecker product of two sparse matrices matrix C = kron(A, B)
  • add lot of finite element on Mesh3, MeshS, MeshL of Discontinous Galerling Element in 3d : P1dc3d, P2dc3d, P3dc3d, P4dc3d , P0edge3d ,P0edgedc3d , P0face3d ,P0facedc3d , P0VF3d ,P0VFdc3d , on Surface : P1dcS, P2dcS, P3dcS, P4dcS , P0edgeS ,P0edgedcS , P0VFS ,P0VFdcS, on Curve : P1dcL, P2dcL, P3dcL, P4dcL , P0VFL ,P0VFdcL remark; the associated generic name existe of P1dc, P2dc, P0edge, P0VF and all dc finite element corresponding to no continuity across element.
  • add code of intallfaces to do Discontinous Galerkin formulation in 3d (in test FH.)
  • add dist function to a mesh , meshL, MeshS or mesh3
  • signeddistfunction to a meshL or meshS
  • add buildmesh functon to build a 2d mesh from a meshL (same as buildmesh see examples/3dCurve/border.edp)

Changed

  • Now the order to find MPI in configure is first if you have PETSC then take MPI from PETSc otherwise use previous method
  • on MeshL defined with buildmeshL now the default label are 2k-1 (resp. 2k) for the begin (resp. end) of curve where k is the order of curve use in buildmeshL. So if you have one curve the labels are 1 and 2. And new the element label are te region number not the label. This element are not really test so be carfull.
  • PETSc 3.15.0

Deprecated

Removed

Fixed

  • bug in Find triangle contening point in 2d (border case), int Mesh::DataFindBoundary::Find(R2 PP,R *l,int & outside) const the parameter l not correclty return due to local variable.
  • set CFLAGS=-Wno-implicit-function-declaration to complie with Apple clang version 12.0.0 (clang-1200.0.32.29) to remove following error: implicit declaration of function correct 3dCurve/basicGlue.edpand add missing test
  • bugs in SLEPc SVDSolve() with a rectangular Mat
  • bugs in nElementonB for DG 3d formulation.

4.8

Added

  • Bilaplacian example using Morley FE with PETSc, see examples/hpddm/bilaplacian-2d-PETSc.edp
  • Oseen problem preconditioned by PCD, see examples/hpddm/oseen-2d-PETSc.edp
  • SLEPc polynomial eigenvalue solver PEPSolve()
  • add trivial example to check periodic boundary condition on meshS , meshL , mesh3 examples/3d/periodic3.edp examples/3dSurf/periodicS.edp examples/3dCurve/periodicL.edp

Changed

  • PETSc version 3.14.2
  • Mmg version 5.5.2
  • link of ffglut so change in configure.ac and Makefile.am LIBS -> FF_LIBS and LIBS become empty to remove default libs
  • change number of save plot in ffglut from 10 to 20 for O. Pironneau

Fixed

  • some memory leaks
  • the periodic boundary condition have wrong before first a sementic level of MeshS and MeshL case. the new syntexe is for example: meshL Tl=segment(10); fespace Vl(Tl,P1,periodic=[[1],[2]]); meshS Th=square3(10,10,[x2pi,y2pi]); fespace Vh2(Th,P1,periodic=[[1,x],[3,x],[2,y],[4,y]]);
  • fixed '*' keyboard trick, to keep the viewpoint in ffglut or not.

4.7-1

Changed

  • change the language definition to use type as a construction function with named arguments for bem plugin
  • PETSc version 3.14.0
  • ARPACK compiled by SLEPc
  • Mmg version 5.5.0
  • -std=c++14 instead of -std=c++11 when possible

Removed

  • plugins thresholdings, symmetrizeCSR, and fflapack and associed example

Fixed

  • problem compilation with gfortran-10 of arpack and mumps (add -fallow-argument-mismatch flags)

4.7

Added

  • new way to build matrix between 2d Finite element 2d and Curve finite element to do mortar (Thank to Axel ) , see first example examples/tutorial/mortar-DN-4-v4.5.edp
  • add Ns normal vector in R^3 on meshS (normal of the surface) of current point (to day Ns of [x,y,0] plan is [0,0,-1]) no be compatible to exterior normal.
  • add Tl tangent vector in R^3 on meshL (tangent vector of the line/curve) of current point
  • compile ffmaster / ffslave example under windows (thanks to [email protected])
  • Boolean parameter spiltpbedge in buildmesh to split in to edge with two boundary vertices
  • interface to PETSc DMPlex, see examples/hpddm/DMPlex-PETSc.edp
  • function MatDestroy
  • function MatPtAP and transferMat for parallel interpolation between non-matching grids, see examples/hpddm/PtAP-2d-PETSc.edp or examples/hpddm/diffusion-mg-2d-PETSc.edp
  • preliminary interface to SVDSolve from SLEPc to compute singular value decompositions, see examples/hpddm/mf-2d-SLEPc.edp or examples/hpddm/helmholtz-2d-SLEPc-complex.edp
  • preliminary interface to NEPSolve from SLEPc to solve nonlinear eigenvalue problems, see examples/hpddm/nonlinear-2d-SLEPc-complex.edp
  • transpose parameter when constructing a Mat for defining a matrix-free transposed operation
  • interface to PetscMemoryGetCurrentUsage
  • add P2b, RT0, RT1 surface FE (P2bS, RT0S, RT1S))
  • add operator interpolate (2d->3d surface)
  • add operator x = A'*b; where x, b are array and A 2 dim array (full matrix) and generate an error in case of b'*A or b'*A expression
  • function MatLoad to load a PETSc Mat from disk, see examples/hpddm/MatLoad-PETSc.edp
  • possibility to assemble a symmetric HMatrix<complex> and to densify a HMatrix<complex> into a Mat<complex>

Changed

  • moved Htool to its new GitHub location
  • ScaLAPACK and MUMPS are not compiled by PETSc anymore if there is no Fortran compiler
  • MPICH is compiled by PETSc if no MPI is detected during configure, see https://community.freefem.org/t/feature-request-use-download-mpich-on-ubuntu/407
  • PETSc version 3.13.5
  • force --with-cudac=0 in make petsc-slepc, see #141
  • change DSL keyword P1dc3dL->P1dcL and P1dc3dS->P1dcS
  • rename view, hasType, changeSchur to respectively ObjectView, HasType, and ChangeSchur

Deprecated

  • rename changeNumbering, globalNumbering, originalNumbering, changeOperator, destroyRecycling, and attachCoarseOperator to respectively ChangeNumbering, GlobalNumbering, OriginalNumbering, ChangeOperator, DestroyRecycling, and AttachCoarseOperator
  • Nt the normal vector of the current (wrong on meshL) use Ns or Tl

Removed

  • augmentation routine from the PETSc plugin
  • MPIF77 variable

Fixed

  • lot of mistake in MeshL element add a example o check lot of thing tutomesh1d.edp
  • fixed problem of change of mesh when rebuild 2d mesh with buildmesh, .... (Thank to P. Jovilet to points this problem)
  • missing METIS library when using SuiteSparse compiled by PETSc
  • missing -fno-stack-protector when building PETSc on Windows, see https://community.freefem.org/t/error-loading-complex-petsc-slepc-library/370
  • fixed ffglut for the plotting of FE array solution
  • fixed ffglut bug on MacOS Catalina , draw inn only half windows screen (Apple Bug ???)
  • correct P0VF finite element
  • abs function of array

4.6

Added

  • new search algorithm for the element containing a point (more safe) in mesh of type mesh3, meshS, or meshL.
  • new function hasType to know if a PETSc component has been installed, e.g., hasType("PC", "hypre")
  • eigenvalue problems on linear elements, cf. examples/eigen/LapEigen1DBeltrami.edp or examples/hpddm/laplace-beltrami-3d-line-SLEPc.edp
  • --download-cmake in PETSc configure if there is no CMake available
  • flags --with-[slepc|slepccomplex]-include and --with-[slepc|slepccomplex]-ldflags for when SLEPc has been built outside of FreeFEM or PETSc
  • interface to KSPSetResidualHistory and KSPGetIterationNumber
  • interface to mpiWaitAll
  • new function extract, allows to build a curve mesh from a 2d mesh (can extract a labeled boundary, apply a geometric transformation)
  • ffglut can plot a vectorial FE function in surface 3d
  • distributed ParMmg interface, cf. examples/hpddm/distributed-parmmg.edp or examples/hpddm/laplace-adapt-dist-3d-PETSc.edp
  • new parallel interpolator on non-matching meshes, cf. examples/hpddm/transfer.edp
  • ability to solve problems in single precision or with 64 bit integers
  • tool to read data form vtk file only in 3d (cf. plugin iovtk a first example examples/plugin/iovtk.edp)
  • tool to read/wrile ply file of meshL, mesh3, MeshS : Polygon File Format / Stanford Triangle Format do load "ioply" see examples//3dSurf/operatorsOnMeshS.edp

Changed

  • new tgv values: -10 => zero row, -20 => zero row/column
  • Windows binary now shipped with PETSc/SLEPc
  • BEM examples are now in examples/mpi
  • plot border type is now in 3d (border 2d and 3d)
  • PETSc version 3.13.0

Fixed

  • --enable-download_package may now be used to download a single package, e.g., --enable-download_metis
  • compilation of PETSc under Windows
  • compilation of plugins when using static libraries
  • correct detection problem in FE type when use a vectorial FE
  • macro concatenation with spaces in arguments
  • correct bug in plugin/seq/Schur-Complement.cpp
  • correct ambiguity bug in plugin/seq/bfstream.cpp (reading real or integer)
  • compilation of plugin libff-mmap-semaphore.c under windows

4.5

Added

  • for windows version: rename under mpi MUMPS in MUMPS_mpi and in sequentiel in MUMPS_seq due to conflict between seq. and mpi version so all MUMPS load become MUMPS_seq or MUMPS_mpiin all examples
  • correct link edition with fortran mpi under windows juste use the msmpi (just use libmsmpi.dll )
  • new mmg and parmmg (parallel mmg) plugins interfacing mmg5 and parmmg libraries, to replace mmg3d-v4.0 and freeyams (Thanks to P-H Tournier)
  • a true 3d anisotropic mesh adaptation examples/3d/Laplace-Adapt-aniso-3d.edp
  • an example to extract surface mesh from isovalue in examples/3dSurf/Pinochio.edp
  • function f.eatspace to reach eof on istream file which return false in case of EOF.
  • function f.length to get the istream file length
  • Interface to PetscLogStagePush()/PetscLogStagePop()
  • Ability to directly assemble a Mat using a varf
  • New bem plugin for the Boundary Element Method (using htool and BemTool libraries)
  • New DSL for BEM (varfbem see examples/bem)
  • add int0d to apply Neumann BC (curve FE), differential operators (dx,dy,...), compute an 1d integral
  • add P1dc FE for Border FEM (possible to define a new FE with plugin)
  • PETSc as a subdomain solver for HPDDM

Changed

  • correct ffglut (bug in case of changing number of nb isovalue)
  • PETSc version 3.12.4
  • Change the point search triangle algorithm to be sure in any case (in test)
  • Sline operator renamed to segment
  • In square3, segment, movemesh functions: geometry transformation can now be [X] or [X,Y] or [X,Y,Z] according to the minimal shape element dim
  • PETSc now download OpenBLAS if there is no BLAS found by FreeFEM configure

Deprecated

  • freeyams plugin
  • mmg3d-v4.0 plugin

Fixed

  • fix plot for curve mesh

4.4-3

Added

  • Preliminary support for symmetric distributed PETSc matrices (MATMPISBAIJ instead of MATMPIAIJ)
  • Interface to AMS, Hiptmair--Xu preconditioner for problems in H(curl), see maxwell-3d-PETSc.edp
  • FEM on curve 3D (in test)
  • P0, P1, P2 curve 3D FE (scalar for the moment)
  • i/o medit and vtk format for curve FE
  • checkMesh() function, allow to remove multiple vertices, elements and border elements (argument: precisvertice(double),removeduplicate(bool))
  • possible to build a curve mesh from a surface, ThS = buildBdMesh(ThS) and define this new mesh by meshL ThL= ThS.Gamma
  • can extract a border part of a meshL (meshL ThL = extract(ThL,label=llabs))
  • Support for optimized boundary conditions with PETSc, see helmholtz-2d-PETSc-complex.edp
  • buildmeshL() function: build meshL from borders
  • mpiCommSelf keyword

Changed

  • function buildSurface(...) renamed by buildBdMesh(...)
  • line3(...) renamed by SLine(...)

Removed

  • FFTW is not compiled by PETSc anymore
  • Spurious outputs in TetGen plugin
  • curve3 type -> border
  • hypre examples since it is not downloaded by FreeFEM for many months (use PETSc instead)
  • dscalprod routine from HPDDM and PETSc plugins, use A(u, v) with A a Mat or a schwarz object
  • export function for macro_ddm.idp, use savevtk as in the sequential iovtk plugin

Fixed

  • plotMPI function for plotting 3D solutions, problem with serialize
  • variable mes in clean_mesh function
  • correct bug verflow in plugin iohdf5
  • correct problem with buffer iostream function (buffer must be out of range )
  • correct i/o vtk and by defaut write at binary format
  • fix an overflow in RT13d FE
  • problem with auto-build of border mesh

4.4-2

Added

  • add matrix and array tools (FH) 
     matrix A=eye(10);
 real[int,int] af = eye(10,10);
 real[int,int] a(10,10);
 int[int] I=[1,3,6];
 real[int] d = a.diag ; // get the diag of full matrix (no copy)
 real[int] dI= d(I); // init a array from renumbering array
 real[int] c= a(:,1)(I); // init a array from renumbering array
 real[int] aa= a.asarray; //  view full the matrice  as an array (no copy)
 a(2:5,3:7).diag= 200;
 a.diag += 100;
  • adding of a global variable lockOrientation to allows the building of mesh without checking the orientation elements (AF)
  • add plugin tool to build matrix edge/P1 with sign mat_edgeP1 (FH)
  • new examples diffusion-2d-mg.edp and helmholtz-2d-mg.edp showing how to use user-defined coarse corrections
  • support for nonzero scalars in PETSc block matrices
  • simpler constructor for sequential HPDDM matrices (no need for the restriction array and the partition of unity)
  • array of Mat and schwarz types
  • add mpi meshS (serialize object)

Changed

  • correct mistake in mpirank in case of broadcast with comm (thank tp PHT)
  • update fftw to v3.3.8 and openblas v0.3.6
  • in movemesh23 correct the argument label -> region to change label
  • new implementation for the moving mesh functions, new arguments: boolean cleanmesh, removemultiple, rebuildborder
  • new PETSc version 3.12
  • templatize movemesh, setMesh functions
  • add conditional tests in make check

Fixed

  • spurious output in PARDISO
  • fix problem in ffglut (AF)
  • detect hdf5 and gsl if no enable-download

4.4

Added

  • interface to TSSolve, DAE/ODE solvers from PETSc
  • interface to TaoSolve, Toolkit for Advance Optimization from PETSc
  • simpler constructor for sequential PETSc matrices (no need for the restriction array and the partition of unity)
  • some unit tests

Changed

  • PETSc version 3.11.3
  • replaced custom implementations (RNM::real, RNM::norm2, and Fem2D::norm) by C++11 functions
  • API of the macro plotMPI
  • switched to inexact coarse operators in HPDDM by default
  • RHS and solution vectors permuted in IterativeMethod and DDM
  • .mesh are now saved using version 2 (which stores floating-point scalars in double precision)

Removed

  • legacy linear solver interfaces using the old matrix type
  • dot products using CBLAS because of errors at link time
  • Newtow function (bad name)

Fixed

  • assertion failure with some 3D meshes when doing trunc(Th, true) (thanks to F. Feppon)
  • compile error when plotting arrays of vectorial functions

4.2.1

Added

  • nested fieldsplit example examples/hpddm/natural-convection-2d-PETSc-fieldsplit.edp
  • int[int][int] array; is now supported (a size was previously needed, i.e., array(0);)
  • check selectivity during make check, depending on available 3rd party librairies
  • new CI/CD tools for develop branch
  • new gestion of mesh3 - meshS coupling
  • square3, buildSurface... operators for meshS

Changed

  • SLEPc is now directly downloaded by PETSc with --download-slepc
  • HPDDM and PETSc API have been simplified, instead of an int[int] and an int[int][int], only a single int[int][int] is now needed
  • build macros for HPDDM and PETSc have been simplified to follow the above API change, two parameters have been permuted as well to match the HPDDM and PETSc constructors
  • PETSc version 3.11.2 and HPDDM with multilevel GenEO

Removed

  • old interfaces that were not maintained anymore (pARMS, PaStiX, hips) and that are available through PETSc
  • spurious outputs when destroying some meshes
  • old surface msh3 type, replaced by meshS

Fixed

  • multiple segmentation faults when using unitialized values (thanks to G. Sadaka)
  • nested fieldsplits in the PETSc interface
  • memory leaks in SNESSolve (nonlinear PETSc solvers)
  • bug fix of Cofactor function
  • various bug fixes on surface mesh

4.1

Fixed

  • missing conj operation is some hermitian operation on complex sparse matrix like A+cB', AB' thanks to P-H Tournier
  • writing CheckAllEdp to be compatible with new tree
  • fix eps in trunc in case of very anisotrope mesh, Thank G. Sadaka

Added

Changed

  • PETSc/SLEPc version 3.11

4.0

  • correct bug in RT1Ortho and RT2Ortho 2d in the computation of derivative (2018-01-30, Thank to [email protected])
  • uniformize 2d/3d in element, method EdgeOrientation(e) now return +1/-1
  • change all the sparse matrix structure
  • remove all map matrix jan 2019

Added

  • surface finite element (in progress)
  • nElementonB (version 2d and 3d of nTonEge)
  • area ( same the lenEdge) in 3d
  • add function labels to get the array of label of a mesh
  • add function regions to get the array of label of a mesh
  • correct big bug in toRarray,toZarray, toCarray transform [ ... ] array to int[int], real[int], complex[int]

3.62 - 2018-08-31

Added

  • add x0=true/false, add veps=eps in solver parameters to initialazed of not the the CG , GMRES algo with 0 or previous value and veps is to get the absolue tolerance
  • A tool of solve adjoint matrix A with only one single LU decomposition with LU, UMFPACK, GMRES
    u[]=A'^-1*b;
  • Add plugin to save matrix in Harwell-Boeing format (see Harwell-Boeing format)

Fixed

  • Fix bug in trunc (2d) in case of very fine mesh (eps too small)

3.61 - 2018-06-20

Added

  • Add name parameter kerneln=, kernelt=, kerneldim= for dissection solver
  • Add option in method toClose function in fquatree to get the nearst point (for intersect meshes)
  • Add missing file curvature.edp
  • Add imax, jmax, imin, jmin to get index of row or column of the min or max coefficient
    We have: A(A.imin,A.jmin) = A.min
  • Add cosmetics in macro (macro name, macro line...)

Changed

  • Pass to PETSc/SLEPc version 3.8.4/3.8.3

Fixed

  • Fix launchff.exe bug under windows 64 to choose a filescrip if no parameter
  • Fix the label definition in case of intalledges in 2d
  • Fix mpi_comm with MUMPS (very rare)

3.60 - 2018-04-13

Changed

  • The main distribution is now on Github