chn/ufo
1
0
Fork 0
mirror of https://github.com/CHN-beta/ufo.git synced 2024-10-22 19:58:44 +08:00
Code Issues Packages Projects Releases Wiki Activity

放宽合并模式和去除模式的条件

Browse Source
This commit is contained in:
陈浩南 2023-09-21 15:40:20 +08:00
parent b30fdd69be
commit ace7793fd7
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
main.cpp
View File

@ -219,14 +219,14 @@ int main(int argc, const char** argv)
_.Source = input.QPointData[i_of_qpoint].QPoint;
if (!input.Debug.value_or(false))
{
// 从小到大枚举所有的模式，并将相近的模式（相差小于 0.01 THz）合并
// 从小到大枚举所有的模式，并将相近的模式（相差小于 0.1 THz）合并
std::map<double, double> frequency_to_weight;
for (unsigned i_of_mode = 0; i_of_mode < input.QPointData[i_of_qpoint].ModeData.size(); i_of_mode++)
{
auto frequency = input.QPointData[i_of_qpoint].ModeData[i_of_mode].Frequency;
auto weight = projection_coefficient[i_of_qpoint][i_of_mode][i_of_sub_qpoint];
auto it_lower = frequency_to_weight.lower_bound(frequency - 0.01);
auto it_upper = frequency_to_weight.upper_bound(frequency + 0.01);
auto it_lower = frequency_to_weight.lower_bound(frequency - 0.1);
auto it_upper = frequency_to_weight.upper_bound(frequency + 0.1);
if (it_lower == it_upper)
frequency_to_weight[frequency] = weight;
else
@ -241,9 +241,9 @@ int main(int argc, const char** argv)
frequency_to_weight[frequency_sum / weight_sum] = weight_sum;
}
}
// 仅保留权重大于 0.01 的模式
// 仅保留权重大于 0.1 的模式
for (auto& mode : frequency_to_weight)
if (mode.second > 0.01)
if (mode.second > 0.1)
{
auto& __ = _.ModeData.emplace_back();
__.Frequency = mode.first;