- The Lis solver has been updated to the latest version 2.1.6.
- The MUMPS solver has been updated to the latest version 5.7.3.
- Bug fixes.
- A new maximum shear stress fracture propagation criterion was added (*CFCP=4, see PhiPsi keywords manual for details).
- The compiler has been updated to GCC version 14.1.0.

- Bug fixes and enhancements of 3D XFEM.
- Keywords file now supports parameter definition and four arithmetic operations (+,-,*, and /).
- Now you can start a PhiPsi simulation in Terminal using the following instruction, for example:
`PhiPsi_Win64.exe -i XFEM.kpp -n 4 where -n specifies the number of threads for the OpenMP parallel simulation.`

- Enhancements of 3D XFEM.
- Support output vtk files for 3D cracks (Set *Key_Save_vtk to 1).
- Compiler updated to GCC version 13.1.0.
- New method to calculate 3D crack aperture (Set *Key_Crack_Aperture_Method to 1 to active).
- Support degenerated hexahedral elements.
- Some bugs on data conflicting in OpenMP have been fixed.

- Now supports intersection of arbitrary 3D cracks and the simulation of 3D complex fracture network.
- Add a new linear system solver: PCG-EBE (Set *Key_SLOE to 11).
- Supports simulation of contact between3D crack surfaces.
- Dynamic memory management for 3D XFEM simulation.
- Supports generation of random 3D natural fractures.
- Supports the import of FracMan natural cracks file (*.fab file).
- Bug fixes.
- Optimization.

- Kill elements to simulate excavation.
- Bug fixes.
- Optimization.

- Propagation of 3D cracks.

- A new linear solver SuperLU added.
- README.txt added.

- Support field problems with cracks and holes.

- Fixes a critical bug on stiffness matrix assembly of different kinds of enriched elements.
- Support random hole generation.
- Support cracks emerged form holes.

- Support triangular elements degenerated from quadrilateral elements.
- Improve the interactive display of error and warning information.
- Support nonlinear analysis such as plastic analysis.

- Add a linear system solver called Lis which supports OpenMP.
- Computational efficiency optimization (up to 50% as least).
- Add the support for the interactions between cracks and polygonal inclusions.

- Add the peridymanic analysis for simple problems.
- Add the support for the random generation of the irregular polygon inclusions.
- Add a new random generation algorithm which is controlled by a seed number.
- Improvements on the save of the binary data.
- Adds a new integration rule, i.e., subquads integration.
- Change the suffix of the keywords file to kpp.
- Add the support for curved cracks.
- Improvement: expand dimensions of related matrices when calculating the B matrix.

- Fix the bug in subroutine Get_Node_Stress_XFEM.f.
- Improvements of subroutines Vector_belongs_Matrix_Is_Dou.f, Vector_belongs_Matrix_Is_Int.f, and Matrix_Count_Row_Int.f.
- Rewrite the subroutine Tool_Intersection_Line_and_Circle.f with a totally new scheme.
- Improvement: clear all the results files in the current work directory by setting the keyword *Key_Clear_All to be 1.
- Add the support for the simulation of the intersection of cracks and holes.
- Add the support for the crossing crack.
- Improvement: identify the input parameters of the executable program by using "getarg" function.

- Fix a bug for the analysis of the static field problems.
- Add the analysis type for the transient field problems.
- Add the explicit dynamic analysis.
- Enhance OpenMP supports in the following subroutines: a) Get_Node_Stress_FEM.f b) Get_Gauss_Stress_FEM.f c) Assemble_Mass_Matrix_FEM.f d) Assemble_Stiffness_Matrix_FEM.f e) Get_Gauss_Disps.f
- Add the molecular dynamics analysis.
- Add the support for the plasticity material and damaged material.
- Add the support for the cohesive crack.
- Improvement: propagation increment is now allowed to be less than the length of an element.

- Add 3D XFEM analysis for the simple 3D planar crack.
- Add the support for the coupling of degrees of freedom (DOFs).
- Add the implicit dynamic analysis.
- Add the support the thermal stress.
- Add the support for the Biot consolidation.
- Add the analysis type for the static field problems.