zhangyiss/gctl_potential
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

tmp

Browse Source
This commit is contained in:
张壹 2024-11-28 13:24:07 +08:00
parent 28f6e37527
commit 76a3acc06c
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
lib/potential/mkernel_tricone.cpp
View File

@ -212,8 +212,8 @@ void gctl::magkernel(matrix<double> &kernel, const array<mag_tricone> &top_ele,
for (j = 0; j < e_size; j++) for (j = 0; j < e_size; j++)
{ {
mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i)); mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i));
if (comp_type == Bx) kernel[i][j] = mag_b.z; if (comp_type == Bx) kernel[i][j] = mag_b.y;
if (comp_type == By) kernel[i][j] = mag_b.y; if (comp_type == By) kernel[i][j] = mag_b.z;
if (comp_type == Bz) kernel[i][j] = mag_b.x; if (comp_type == Bz) kernel[i][j] = mag_b.x;
} }
} }
@ -227,8 +227,8 @@ void gctl::magkernel(matrix<double> &kernel, const array<mag_tricone> &top_ele,
for (j = 0; j < e_size; j++) for (j = 0; j < e_size; j++)
{ {
mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i)); mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i));
if (comp_type == Bx) kernel[i][j] -= mag_b.z; if (comp_type == Bx) kernel[i][j] -= mag_b.y;
if (comp_type == By) kernel[i][j] -= mag_b.y; if (comp_type == By) kernel[i][j] -= mag_b.z;
if (comp_type == Bz) kernel[i][j] -= mag_b.x; if (comp_type == Bz) kernel[i][j] -= mag_b.x;
} }
} }
@ -300,8 +300,8 @@ void gctl::magkernel(spmat<double> &kernel, const array<mag_tricone> &top_ele, c
if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0) if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0)
{ {
mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i)); mag_b = magkernel_single(top_ele[j], obsp[i], Rs.get(i));
if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val = mag_b.z; if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val = mag_b.y;
if (comp_type == By) triplts[tri_idx[j + i*e_size]].val = mag_b.y; if (comp_type == By) triplts[tri_idx[j + i*e_size]].val = mag_b.z;
if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val = mag_b.x; if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val = mag_b.x;
} }
} }
@ -319,8 +319,8 @@ void gctl::magkernel(spmat<double> &kernel, const array<mag_tricone> &top_ele, c
if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0) if (geometry3d::angle(obsp[i].s2c(), cen) < cut_angle*GCTL_Pi/180.0)
{ {
mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i)); mag_b = magkernel_single(btm_ele[j], obsp[i], Rs.get(i));
if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val -= mag_b.z; if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val -= mag_b.y;
if (comp_type == By) triplts[tri_idx[j + i*e_size]].val -= mag_b.y; if (comp_type == By) triplts[tri_idx[j + i*e_size]].val -= mag_b.z;
if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val -= mag_b.x; if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val -= mag_b.x;
} }
} }