Skip to content
体验新版
项目
组织
正在加载...
登录
切换导航
打开侧边栏
Greenplum
Opencv
提交
235a6784
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,发现更多精彩内容 >>
提交
235a6784
编写于
4月 04, 2013
作者:
A
Andrey Kamaev
浏览文件
操作
浏览文件
下载
电子邮件补丁
差异文件
SVD: always update W vector for better algorithm convergency
上级
1e9ed142
变更
1
隐藏空白更改
内联
并排
Showing
1 changed file
with
9 addition
and
27 deletion
+9
-27
modules/core/src/lapack.cpp
modules/core/src/lapack.cpp
+9
-27
未找到文件。
modules/core/src/lapack.cpp
浏览文件 @
235a6784
...
...
@@ -577,10 +577,10 @@ JacobiSVDImpl_(_Tp* At, size_t astep, _Tp* _W, _Tp* Vt, size_t vstep,
continue
;
p
*=
2
;
double
beta
=
a
-
b
,
gamma
=
hypot
((
double
)
p
,
beta
)
,
delta
;
double
beta
=
a
-
b
,
gamma
=
hypot
((
double
)
p
,
beta
);
if
(
beta
<
0
)
{
delta
=
(
gamma
-
beta
)
*
0.5
;
d
ouble
d
elta
=
(
gamma
-
beta
)
*
0.5
;
s
=
(
_Tp
)
std
::
sqrt
(
delta
/
gamma
);
c
=
(
_Tp
)(
p
/
(
gamma
*
s
*
2
));
}
...
...
@@ -588,36 +588,18 @@ JacobiSVDImpl_(_Tp* At, size_t astep, _Tp* _W, _Tp* Vt, size_t vstep,
{
c
=
(
_Tp
)
std
::
sqrt
((
gamma
+
beta
)
/
(
gamma
*
2
));
s
=
(
_Tp
)(
p
/
(
gamma
*
c
*
2
));
delta
=
p
*
p
*
0.5
/
(
gamma
+
beta
);
}
W
[
i
]
+=
delta
;
W
[
j
]
-=
delta
;
if
(
iter
%
2
!=
0
&&
W
[
i
]
>
0
&&
W
[
j
]
>
0
)
{
k
=
vblas
.
givens
(
Ai
,
Aj
,
m
,
c
,
s
);
for
(
;
k
<
m
;
k
++
)
{
_Tp
t0
=
c
*
Ai
[
k
]
+
s
*
Aj
[
k
];
_Tp
t1
=
-
s
*
Ai
[
k
]
+
c
*
Aj
[
k
];
Ai
[
k
]
=
t0
;
Aj
[
k
]
=
t1
;
}
}
else
a
=
b
=
0
;
for
(
k
=
0
;
k
<
m
;
k
++
)
{
a
=
b
=
0
;
for
(
k
=
0
;
k
<
m
;
k
++
)
{
_Tp
t0
=
c
*
Ai
[
k
]
+
s
*
Aj
[
k
];
_Tp
t1
=
-
s
*
Ai
[
k
]
+
c
*
Aj
[
k
];
Ai
[
k
]
=
t0
;
Aj
[
k
]
=
t1
;
_Tp
t0
=
c
*
Ai
[
k
]
+
s
*
Aj
[
k
];
_Tp
t1
=
-
s
*
Ai
[
k
]
+
c
*
Aj
[
k
];
Ai
[
k
]
=
t0
;
Aj
[
k
]
=
t1
;
a
+=
(
double
)
t0
*
t0
;
b
+=
(
double
)
t1
*
t1
;
}
W
[
i
]
=
a
;
W
[
j
]
=
b
;
a
+=
(
double
)
t0
*
t0
;
b
+=
(
double
)
t1
*
t1
;
}
W
[
i
]
=
a
;
W
[
j
]
=
b
;
changed
=
true
;
...
...
编辑
预览
Markdown
is supported
0%
请重试
或
添加新附件
.
添加附件
取消
You are about to add
0
people
to the discussion. Proceed with caution.
先完成此消息的编辑!
取消
想要评论请
注册
或
登录