Hi @Chenxi_Fu
Thanks for trying the software. It is possible that there is a memory leak somewhere in the software, including the solver.
If you want to try increasing the voltage with larger steps, and you are using some of the existing scripts, you may try switching from logarithmic to default damping using:
opts = get_equation_command(name="PotentialEquation", device=device, region=region)
opts['variable_update']="default"
equation(**opts)
After the simulation has been setup.
The reset_devsim command can be used to clear the internal storage. However, this would require reloading your mesh and physics afterward. If that doesn’t help, then it is possible there is a leak within devsim, or even the solver.
If you are using the Intel MKL, it may be possible to disable their memory manager, by setting the appropriate environment variable before running your script:
https://www.intel.com/content/www/us/en/docs/onemkl/developer-reference-c/2023-1/mkl-disable-fast-mm.html
It is also possible to use UMFPACK as an alternative solver, to help diagnose the issue. By running the script with:
python -mdevsim.umfpack.umfshim ssac_cap.py
where the last argument is the name of your script.
If you are able to provide a minimum reproducible example, please share it, preferably at:
https://github.com/devsim/devsim/issues
Regards,
Juan