zhangyiss/liblbfgs
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Added owlqn_project() function.

Browse Source 
git-svn-id: file:///home/svnrepos/software/liblbfgs/trunk@19 ecf4c44f-38d1-4fa4-9757-a0b4dd0349fc
This commit is contained in:
naoaki 2008-07-10 06:17:35 +00:00
parent 85e940b2ca
commit 36b9ae4a62
1 changed files with 25 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
lib/lbfgs.c
View File

@ -192,6 +192,13 @@ static void owlqn_direction(
const int n const int n
); );
static void owlqn_project(
lbfgsfloatval_t* d,
const lbfgsfloatval_t* sign,
const int start,
const int n
);
#if defined(USE_SSE) && (defined(__SSE__) || defined(__SSE2__)) #if defined(USE_SSE) && (defined(__SSE__) || defined(__SSE2__))
static int round_out_variables(int n) static int round_out_variables(int n)
@ -600,11 +607,7 @@ static int line_search_backtracking(
if (param->orthantwise_c != 0.) { if (param->orthantwise_c != 0.) {
/* The current point is projected onto the orthant of the initial one. */ /* The current point is projected onto the orthant of the initial one. */
for (i = param->orthantwise_start;i < n;++i) { owlqn_project(x, xp, param->orthantwise_start, n);
if (x[i] * xp[i] < 0.) {
x[i] = 0.;
}
}
} }
/* Evaluate the function and gradient values. */ /* Evaluate the function and gradient values. */
@ -767,11 +770,7 @@ static int line_search_morethuente(
if (param->orthantwise_c != 0.) { if (param->orthantwise_c != 0.) {
/* The current point is projected onto the orthant of the previous one. */ /* The current point is projected onto the orthant of the previous one. */
for (i = param->orthantwise_start;i < n;++i) { owlqn_project(x, wa, param->orthantwise_start, n);
if (x[i] * wa[i] < 0.) {
x[i] = 0.;
}
}
} }
/* Evaluate the function and gradient values. */ /* Evaluate the function and gradient values. */
@ -1280,3 +1279,19 @@ static void owlqn_direction(
} }
} }
} }
static void owlqn_project(
lbfgsfloatval_t* d,
const lbfgsfloatval_t* sign,
const int start,
const int n
)
{
int i;
for (i = start;i < n;++i) {
if (d[i] * sign[i] < 0) {
d[i] = 0;
}
}
}