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

update src

Browse Source
This commit is contained in:
张壹
2021-07-28 15:47:18 +08:00
parent 5b65eab987
commit 755282373a
4 changed files with 19608 additions and 19608 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13122
data/prism/prism_mag_grad.txt
View File

File diff suppressed because it is too large Load Diff

12970
data/prism/prism_mag_pot.txt
View File

File diff suppressed because it is too large Load Diff

13122
data/prism/prism_mag_tensor.txt
View File

File diff suppressed because it is too large Load Diff

2
src/magtet.cpp
View File

@@ -144,7 +144,7 @@ void magtet::cal_tensors()
v3 = *ele_[e].fget(i, (j+1)%3) - *ele_[e].fget(i, j); v3 = *ele_[e].fget(i, (j+1)%3) - *ele_[e].fget(i, j);
ne = gctl::cross(v3, nf).normal(); ne = gctl::cross(v3, nf).normal();
ele_para_[e].enorm[j+i*3] = ne; ele_para_[e].enorm[j+i*3] = ne;
ele_para_[e].etang[j+i*3] = gctl::cross(nf, ne); ele_para_[e].etang[j+i*3] = gctl::cross(nf, ne).normal();
} }
ele_para_[e].mag_amp[i] = gctl::dot(magz_[e], nf); ele_para_[e].mag_amp[i] = gctl::dot(magz_[e], nf);