Skip to content
体验新版
项目
组织
正在加载...
登录
切换导航
打开侧边栏
Greenplum
Opencv
提交
5e92a777
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,发现更多精彩内容 >>
提交
5e92a777
编写于
1月 19, 2015
作者:
V
Vadim Pisarevsky
浏览文件
操作
浏览文件
下载
差异文件
Merge pull request #3569 from ilya-lavrenov:sse_mul
上级
17b1152f
68962adc
变更
1
隐藏空白更改
内联
并排
Showing
1 changed file
with
159 addition
and
0 deletion
+159
-0
modules/core/src/arithm.cpp
modules/core/src/arithm.cpp
+159
-0
未找到文件。
modules/core/src/arithm.cpp
浏览文件 @
5e92a777
...
...
@@ -2355,6 +2355,165 @@ struct Mul_SIMD<float, float>
}
};
#elif CV_SSE2
#if CV_SSE4_1
template
<
>
struct
Mul_SIMD
<
ushort
,
float
>
{
Mul_SIMD
()
{
haveSSE
=
checkHardwareSupport
(
CV_CPU_SSE4_1
);
}
int
operator
()
(
const
ushort
*
src1
,
const
ushort
*
src2
,
ushort
*
dst
,
int
width
,
float
scale
)
const
{
int
x
=
0
;
if
(
!
haveSSE
)
return
x
;
__m128i
v_zero
=
_mm_setzero_si128
();
if
(
scale
!=
1.0
f
)
{
__m128
v_scale
=
_mm_set1_ps
(
scale
);
for
(
;
x
<=
width
-
8
;
x
+=
8
)
{
__m128i
v_src1
=
_mm_loadu_si128
((
__m128i
const
*
)(
src1
+
x
));
__m128i
v_src2
=
_mm_loadu_si128
((
__m128i
const
*
)(
src2
+
x
));
__m128
v_dst1
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_unpacklo_epi16
(
v_src1
,
v_zero
)),
_mm_cvtepi32_ps
(
_mm_unpacklo_epi16
(
v_src2
,
v_zero
)));
v_dst1
=
_mm_mul_ps
(
v_dst1
,
v_scale
);
__m128
v_dst2
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_unpackhi_epi16
(
v_src1
,
v_zero
)),
_mm_cvtepi32_ps
(
_mm_unpackhi_epi16
(
v_src2
,
v_zero
)));
v_dst2
=
_mm_mul_ps
(
v_dst2
,
v_scale
);
__m128i
v_dsti
=
_mm_packus_epi32
(
_mm_cvtps_epi32
(
v_dst1
),
_mm_cvtps_epi32
(
v_dst2
));
_mm_storeu_si128
((
__m128i
*
)(
dst
+
x
),
v_dsti
);
}
}
return
x
;
}
bool
haveSSE
;
};
#endif
template
<
>
struct
Mul_SIMD
<
schar
,
float
>
{
Mul_SIMD
()
{
haveSSE
=
checkHardwareSupport
(
CV_CPU_SSE2
);
}
int
operator
()
(
const
schar
*
src1
,
const
schar
*
src2
,
schar
*
dst
,
int
width
,
float
scale
)
const
{
int
x
=
0
;
if
(
!
haveSSE
)
return
x
;
__m128i
v_zero
=
_mm_setzero_si128
();
if
(
scale
==
1.0
f
)
for
(
;
x
<=
width
-
8
;
x
+=
8
)
{
__m128i
v_src1
=
_mm_loadl_epi64
((
__m128i
const
*
)(
src1
+
x
));
__m128i
v_src2
=
_mm_loadl_epi64
((
__m128i
const
*
)(
src2
+
x
));
v_src1
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src1
),
8
);
v_src2
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src2
),
8
);
__m128
v_dst1
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src2
),
16
)));
__m128
v_dst2
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src2
),
16
)));
__m128i
v_dsti
=
_mm_packs_epi32
(
_mm_cvtps_epi32
(
v_dst1
),
_mm_cvtps_epi32
(
v_dst2
));
_mm_storel_epi64
((
__m128i
*
)(
dst
+
x
),
_mm_packs_epi16
(
v_dsti
,
v_zero
));
}
else
{
__m128
v_scale
=
_mm_set1_ps
(
scale
);
for
(
;
x
<=
width
-
8
;
x
+=
8
)
{
__m128i
v_src1
=
_mm_loadl_epi64
((
__m128i
const
*
)(
src1
+
x
));
__m128i
v_src2
=
_mm_loadl_epi64
((
__m128i
const
*
)(
src2
+
x
));
v_src1
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src1
),
8
);
v_src2
=
_mm_srai_epi16
(
_mm_unpacklo_epi8
(
v_zero
,
v_src2
),
8
);
__m128
v_dst1
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src2
),
16
)));
v_dst1
=
_mm_mul_ps
(
v_dst1
,
v_scale
);
__m128
v_dst2
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src2
),
16
)));
v_dst2
=
_mm_mul_ps
(
v_dst2
,
v_scale
);
__m128i
v_dsti
=
_mm_packs_epi32
(
_mm_cvtps_epi32
(
v_dst1
),
_mm_cvtps_epi32
(
v_dst2
));
_mm_storel_epi64
((
__m128i
*
)(
dst
+
x
),
_mm_packs_epi16
(
v_dsti
,
v_zero
));
}
}
return
x
;
}
bool
haveSSE
;
};
template
<
>
struct
Mul_SIMD
<
short
,
float
>
{
Mul_SIMD
()
{
haveSSE
=
checkHardwareSupport
(
CV_CPU_SSE2
);
}
int
operator
()
(
const
short
*
src1
,
const
short
*
src2
,
short
*
dst
,
int
width
,
float
scale
)
const
{
int
x
=
0
;
if
(
!
haveSSE
)
return
x
;
__m128i
v_zero
=
_mm_setzero_si128
();
if
(
scale
!=
1.0
f
)
{
__m128
v_scale
=
_mm_set1_ps
(
scale
);
for
(
;
x
<=
width
-
8
;
x
+=
8
)
{
__m128i
v_src1
=
_mm_loadu_si128
((
__m128i
const
*
)(
src1
+
x
));
__m128i
v_src2
=
_mm_loadu_si128
((
__m128i
const
*
)(
src2
+
x
));
__m128
v_dst1
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpacklo_epi16
(
v_zero
,
v_src2
),
16
)));
v_dst1
=
_mm_mul_ps
(
v_dst1
,
v_scale
);
__m128
v_dst2
=
_mm_mul_ps
(
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src1
),
16
)),
_mm_cvtepi32_ps
(
_mm_srai_epi32
(
_mm_unpackhi_epi16
(
v_zero
,
v_src2
),
16
)));
v_dst2
=
_mm_mul_ps
(
v_dst2
,
v_scale
);
__m128i
v_dsti
=
_mm_packs_epi32
(
_mm_cvtps_epi32
(
v_dst1
),
_mm_cvtps_epi32
(
v_dst2
));
_mm_storeu_si128
((
__m128i
*
)(
dst
+
x
),
v_dsti
);
}
}
return
x
;
}
bool
haveSSE
;
};
#endif
template
<
typename
T
,
typename
WT
>
static
void
...
...
编辑
预览
Markdown
is supported
0%
请重试
或
添加新附件
.
添加附件
取消
You are about to add
0
people
to the discussion. Proceed with caution.
先完成此消息的编辑!
取消
想要评论请
注册
或
登录