From 76a3acc06c772c9c9c2efdf2d6ecedb17586083c Mon Sep 17 00:00:00 2001 From: Yi Zhang Date: Thu, 28 Nov 2024 13:24:07 +0800 Subject: [PATCH] tmp --- lib/potential/mkernel_tricone.cpp | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/lib/potential/mkernel_tricone.cpp b/lib/potential/mkernel_tricone.cpp index 1e00e48..556429b 100644 --- a/lib/potential/mkernel_tricone.cpp +++ b/lib/potential/mkernel_tricone.cpp @@ -212,8 +212,8 @@ void gctl::magkernel(matrix &kernel, const array &top_ele, for (j = 0; j < e_size; j++) { 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 == By) kernel[i][j] = mag_b.y; + if (comp_type == Bx) 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; } } @@ -227,8 +227,8 @@ void gctl::magkernel(matrix &kernel, const array &top_ele, for (j = 0; j < e_size; j++) { 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 == By) kernel[i][j] -= mag_b.y; + if (comp_type == Bx) 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; } } @@ -300,8 +300,8 @@ void gctl::magkernel(spmat &kernel, const array &top_ele, c 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)); - if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val = mag_b.z; - if (comp_type == By) triplts[tri_idx[j + i*e_size]].val = mag_b.y; + 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.z; if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val = mag_b.x; } } @@ -319,8 +319,8 @@ void gctl::magkernel(spmat &kernel, const array &top_ele, c 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)); - if (comp_type == Bx) triplts[tri_idx[j + i*e_size]].val -= mag_b.z; - if (comp_type == By) triplts[tri_idx[j + i*e_size]].val -= mag_b.y; + 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.z; if (comp_type == Bz) triplts[tri_idx[j + i*e_size]].val -= mag_b.x; } }