Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
CoMISo
CoMISo
Commits
48ab9d26
Commit
48ab9d26
authored
Nov 23, 2016
by
Max Lyon
Browse files
Merge branch 'VCI/master' into linear_dependent_constraint_removal
parents
b0a8180c
bee43685
Changes
1
Hide whitespace changes
Inline
Side-by-side
NSolver/ConstraintTools.cc
View file @
48ab9d26
...
...
@@ -132,12 +132,11 @@ remove_dependent_linear_constraints_only_linear_equality( std::vector<NConstrain
SVectorGMM
col
=
Ac
.
col
(
j
);
// copy row
SVectorGMM
row
;
gmm
::
copy
(
A
.
row
(
i
),
row
);
SVectorGMM
row
=
A
.
row
(
i
);
// iterate over column
typename
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_it
=
gmm
::
vect_const_begin
(
col
);
typename
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_end
=
gmm
::
vect_const_end
(
col
);
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_it
=
gmm
::
vect_const_begin
(
col
);
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_end
=
gmm
::
vect_const_end
(
col
);
for
(;
c_it
!=
c_end
;
++
c_it
)
if
(
row_status
[
c_it
.
index
()]
==
-
1
)
// only process unvisited rows
...
...
@@ -182,8 +181,8 @@ find_max_abs_coeff(SVectorGMM& _v)
unsigned
int
imax
(
0
);
double
vmax
(
0.0
);
typename
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_it
=
gmm
::
vect_const_begin
(
_v
);
typename
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_end
=
gmm
::
vect_const_end
(
_v
);
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_it
=
gmm
::
vect_const_begin
(
_v
);
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
c_end
=
gmm
::
vect_const_end
(
_v
);
for
(;
c_it
!=
c_end
;
++
c_it
)
if
(
c_it
.
index
()
!=
n
-
1
)
...
...
@@ -209,7 +208,7 @@ add_row_simultaneously( int _row_i,
CMatrixGMM
&
_cmat
,
const
double
_eps
)
{
typedef
typename
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
RIter
;
typedef
gmm
::
linalg_traits
<
SVectorGMM
>::
const_iterator
RIter
;
RIter
r_it
=
gmm
::
vect_const_begin
(
_row
);
RIter
r_end
=
gmm
::
vect_const_end
(
_row
);
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment