kumpeishiraishi

kumpeishiraishi

# include <iostream>
# include <chrono>

int main() {
  auto start = std::chrono::system_clock::now(); // use auto to skip typing std::chrono::system_clock::time_point
  /* do something */
  auto end = std::chrono::system_clock::now();
  auto elapsed = std::chrono::duration_cast<std::chrono::nanoseconds>(end-start).count();
  /*
    Trasnlate time unit to nanoseconds.

kumpeishiraishi

やりたいこと

Eigenでdouble型の行列の固有値を計算し、結果を倍精度で標準出力する。

まずはPythonで正解を計算する

import numpy as np
A = np.array([[4,1],[2,-1]])
print(np.linalg.eig(A)[0][0])
print(np.linalg.eig(A)[0][1])

kumpeishiraishi

import numpy as np
from scipy.special import sph_harm
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# Designate quantum numbers here!
m = 1
l = 3

point_size = 100