How to set number of threads that numpy/Scipy uses
OPENBLAS_NUM_THREADS=5 MKL_NUM_THREADS=5 NUMEXPR_MUM_THREAD=5 python contact.py
Tech
kevin-tofu Unfitted FEM: Nitsche Method, CutFEM and Ghost Penalty
Introduction CutFEM (Cut Finite Element Method) is a framework that allows the finite element method to be applied even when…
Benchmarking: Solve Linear Equations with Julia
Introduction Since Julia has great package Gridap, which is able to deal with FEM computation efficiently, Julia lang is one…
Labeling Boundary with Gridap and export as VTU. What is the feature of Julia?
Data Structure in FaceLabeling I struggled with ChatGPT to resolve bugs in my code, but finally gave up. It just…
Solve Linear Elasticity with Gridap, Julia
Introduction Gridap is a highly useful finite element solver built with Julia.However, since the documentation still feels somewhat incomplete, I…
Create simple mesh with Gridap
Lx, Ly, Lz = (5.0, 5.0, 1.0) nx, ny, nz = (10, 10, 2) domain = (0.0,Lx, 0.0,Ly, 0.0,Lz) partition…
Level set method: Update GridapTopOpt Tutorial
Introduction GridapTopOpt is a library that is able to handle level set method based topology optimization. This library is built…
Compile and Install PETSc with external libraries
Download/Install superlu_dist/mumps/MPI/cmake with metis / parmetis without mpi Output Install Activate PETSc Add dir to PATH Check Activation thru GridapPETSc…
Measure elapsed time of process with bash
start=$(date +%s) # processes you wanna measure sleep 3 end=$(date +%s) echo "elapsed time: $((end - start)) sec" start=$(date +%s%N)…