Commit b2a76014 by Zoltan Karsa

fix p

parent aaf31c54
...@@ -24,7 +24,7 @@ def gen_angels_to_pick(n, plot = False): ...@@ -24,7 +24,7 @@ def gen_angels_to_pick(n, plot = False):
parosit = cp.RawKernel(r''' parosit = cp.RawKernel(r'''
extern "C" extern "C"
__global__ __global__
void parosit(const double* x1, const double* x2, double* a, double* b, const int m, const float PI) { void parosit(const double* x1, const double* x2, double* a, double* b, const int m, const double PI) {
int tid = blockDim.x * blockIdx.x + threadIdx.x; int tid = blockDim.x * blockIdx.x + threadIdx.x;
float alpha = x1[tid]; float alpha = x1[tid];
for (int i = 0; i < m; i++) { for (int i = 0; i < m; i++) {
...@@ -43,7 +43,7 @@ void parosit(const double* x1, const double* x2, double* a, double* b, const int ...@@ -43,7 +43,7 @@ void parosit(const double* x1, const double* x2, double* a, double* b, const int
parosit2 = cp.RawKernel(r''' parosit2 = cp.RawKernel(r'''
extern "C" extern "C"
__global__ __global__
void parosit2(const double* x1, const double* x2, double* a, double* b, const int m, const float PI) { void parosit2(const double* x1, const double* x2, double* a, double* b, const int m, const double PI) {
int tid = blockDim.x * blockIdx.x + threadIdx.x; int tid = blockDim.x * blockIdx.x + threadIdx.x;
float alpha = x1[tid]; float alpha = x1[tid];
for (int i = 0; i < m; i++) { for (int i = 0; i < m; i++) {
...@@ -78,9 +78,9 @@ def angles_alap(anglestopick, plot = False): ...@@ -78,9 +78,9 @@ def angles_alap(anglestopick, plot = False):
tgtA = cp.tan(tompa_alpha_arr) tgtA = cp.tan(tompa_alpha_arr)
# hegyes # hegyes
sztgh = tghB*tghA sztgh = tghB * tghA
otgh = tghA + tghB otgh = tghA + tghB
hCy = sztgh/otgh hCy = sztgh / otgh
hCx = hCy / tghA hCx = hCy / tghA
# tompa # tompa
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or sign in to comment