Skip to content
体验新版
项目
组织
正在加载...
登录
切换导航
打开侧边栏
Greenplum
Opencv
提交
1a18a336
O
Opencv
项目概览
Greenplum
/
Opencv
11 个月 前同步成功
通知
7
Star
0
Fork
0
代码
文件
提交
分支
Tags
贡献者
分支图
Diff
Issue
0
列表
看板
标记
里程碑
合并请求
0
DevOps
流水线
流水线任务
计划
Wiki
0
Wiki
分析
仓库
DevOps
项目成员
Pages
O
Opencv
项目概览
项目概览
详情
发布
仓库
仓库
文件
提交
分支
标签
贡献者
分支图
比较
Issue
0
Issue
0
列表
看板
标记
里程碑
合并请求
0
合并请求
0
Pages
DevOps
DevOps
流水线
流水线任务
计划
分析
分析
仓库分析
DevOps
Wiki
0
Wiki
成员
成员
收起侧边栏
关闭侧边栏
动态
分支图
创建新Issue
流水线任务
提交
Issue看板
体验新版 GitCode,发现更多精彩内容 >>
提交
1a18a336
编写于
7月 29, 2016
作者:
K
k-shinotsuka
浏览文件
操作
浏览文件
下载
电子邮件补丁
差异文件
improve SumSqr_SIMD<schar, int, int>()
上级
b34272f8
变更
1
隐藏空白更改
内联
并排
Showing
1 changed file
with
22 addition
and
24 deletion
+22
-24
modules/core/src/stat.cpp
modules/core/src/stat.cpp
+22
-24
未找到文件。
modules/core/src/stat.cpp
浏览文件 @
1a18a336
...
...
@@ -816,38 +816,36 @@ struct SumSqr_SIMD<schar, int, int>
int
x
=
0
;
__m128i
v_zero
=
_mm_setzero_si128
(),
v_sum
=
v_zero
,
v_sqsum
=
v_zero
;
const
int
len_16
=
len
&
~
15
;
for
(
;
x
<=
len
-
16
;
x
+=
16
)
for
(
;
x
<=
len
_16
-
16
;
)
{
__m128i
v_src
=
_mm_loadu_si128
((
const
__m128i
*
)(
src0
+
x
));
__m128i
v_half
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src
),
8
);
__m128i
v_mullo
=
_mm_mullo_epi16
(
v_half
,
v_half
);
__m128i
v_mulhi
=
_mm_mulhi_epi16
(
v_half
,
v_half
);
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_half
),
16
));
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_half
),
16
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpacklo_epi16
(
v_mullo
,
v_mulhi
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpackhi_epi16
(
v_mullo
,
v_mulhi
));
v_half
=
_mm_srai_epi16
(
_mm_unpackhi_epi8
(
v_zero
,
v_src
),
8
);
v_mullo
=
_mm_mullo_epi16
(
v_half
,
v_half
);
v_mulhi
=
_mm_mulhi_epi16
(
v_half
,
v_half
);
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_half
),
16
));
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_half
),
16
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpacklo_epi16
(
v_mullo
,
v_mulhi
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpackhi_epi16
(
v_mullo
,
v_mulhi
));
const
int
len_tmp
=
min
(
x
+
2048
,
len_16
);
__m128i
v_sum_tmp
=
v_zero
;
for
(
;
x
<=
len_tmp
-
16
;
x
+=
16
)
{
__m128i
v_src
=
_mm_loadu_si128
((
const
__m128i
*
)(
src0
+
x
));
__m128i
v_half_0
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src
),
8
);
__m128i
v_half_1
=
_mm_srai_epi16
(
_mm_unpackhi_epi8
(
v_zero
,
v_src
),
8
);
v_sum_tmp
=
_mm_add_epi16
(
v_sum_tmp
,
_mm_add_epi16
(
v_half_0
,
v_half_1
));
__m128i
v_half_2
=
_mm_unpacklo_epi16
(
v_half_0
,
v_half_1
);
__m128i
v_half_3
=
_mm_unpackhi_epi16
(
v_half_0
,
v_half_1
);
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_madd_epi16
(
v_half_2
,
v_half_2
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_madd_epi16
(
v_half_3
,
v_half_3
));
}
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_sum_tmp
),
16
));
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_sum_tmp
),
16
));
}
for
(
;
x
<=
len
-
8
;
x
+=
8
)
{
__m128i
v_src
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
_mm_loadl_epi64
((
__m128i
const
*
)(
src0
+
x
))),
8
);
__m128i
v_half_0
=
_mm_unpackhi_epi64
(
v_src
,
v_src
);
__m128i
v_sum_tmp
=
_mm_add_epi16
(
v_src
,
v_half_0
);
__m128i
v_half_1
=
_mm_unpacklo_epi16
(
v_src
,
v_half_0
);
__m128i
v_mullo
=
_mm_mullo_epi16
(
v_src
,
v_src
);
__m128i
v_mulhi
=
_mm_mulhi_epi16
(
v_src
,
v_src
);
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src
),
16
));
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src
),
16
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpacklo_epi16
(
v_mullo
,
v_mulhi
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_unpackhi_epi16
(
v_mullo
,
v_mulhi
));
v_sum
=
_mm_add_epi32
(
v_sum
,
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_sum_tmp
),
16
));
v_sqsum
=
_mm_add_epi32
(
v_sqsum
,
_mm_madd_epi16
(
v_half_1
,
v_half_1
));
}
int
CV_DECL_ALIGNED
(
16
)
ar
[
8
];
...
...
编辑
预览
Markdown
is supported
0%
请重试
或
添加新附件
.
添加附件
取消
You are about to add
0
people
to the discussion. Proceed with caution.
先完成此消息的编辑!
取消
想要评论请
注册
或
登录