# lap

**Repository Path**: luyi07/lap

## Basic Information

- **Project Name**: lap
- **Description**: LAP: linear algebraic projection.
With given kernel values, the code solves the Hill-Wheeler equation and end up with eigen-values and eigen-functions.
- **Primary Language**: C++
- **License**: MIT
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2023-09-19
- **Last Updated**: 2024-06-15

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

# lap

#### Introduction
LAP is short for linear algebraic projection, used in nuclear structure codes. 
Given parameters Jmax, nPC, ..., the code generates a set of Euler angles, and ask the user to provide kernel values with those angles, in a well defined format. 
Then the code calculates matrix elements out of those kernel values, and solves the eigen-values (energies) and eigen-vectors (wavefunctions).

#### structure

dsygvx.cpp	pack dsygvx from lapack, and use it to eliminate linear denpendence among bases
LAP.cpp 	linear algebraic projection for good angular momentum, given kernel values
UI.py		User Interface, written with tkinter in python
Wigner\_d.py, lap.py	not really used in this version, written in an unsuccessful try with python codes

UI		excutable

#### User's Guide
1. run the code in linux desktop
   * 1.1 Under a linux desktop, in the directory lap/, double click the excutable "UI", or in command line type "./UI".
   * 1.2 press the button "EulerAngles.txt"
       A set of Euler angles will be generated and shown in the text window, and also written into EulerAngles.txt.
   * 1.3 Provide kernel values
       An example is shown in kernels.txt. Calculate your own kernels and format them as in the example into a file. Then modify the entry aside from the button "Show Kernels". Press that button if you'd like to see the contents.
   * 1.4 Run LAP
       Press the button "Run LAP". It may take a while to compute, and show the eigen values in the text window. Eigen values are written into "eig.txt+", while eigen waves are written into "wave.txt+".

2. run the code by command line
   * In directory lap/, find a script "gEuler.sh", to generate Euler angles in EulerAngles.txt 
   * In "gEuler.sh", modify the parameters "Jmax, nPC, ..."
   * Type "bash gEuler.sh", to run.

   * In directory lap/, find a script "projection.sh", to do projection, given kernel values. 
   * In "projection.sh", modify the parameters "Jmax, nPC, ..."
   * Type "bash projection.sh", to run.

#### Join us

欢迎交流讨论！如果有批评建议，欢迎发送到 luyi@qfnu.edu.cn