blob: aadf1fd5fb2cba6c13d399711becaec0703721d0 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
package optimization.stopCriteria;
import optimization.gradientBasedMethods.Objective;
import optimization.gradientBasedMethods.ProjectedObjective;
import optimization.util.MathUtils;
public class ProjectedGradientL2Norm implements StopingCriteria{
/**
* Stop if gradientNorm/(originalGradientNorm) smaller
* than gradientConvergenceValue
*/
protected double gradientConvergenceValue;
public ProjectedGradientL2Norm(double gradientConvergenceValue){
this.gradientConvergenceValue = gradientConvergenceValue;
}
public void reset(){
}
double[] projectGradient(ProjectedObjective obj){
if(obj.auxParameters == null){
obj.auxParameters = new double[obj.getNumParameters()];
}
System.arraycopy(obj.getParameters(), 0, obj.auxParameters, 0, obj.getNumParameters());
MathUtils.minusEquals(obj.auxParameters, obj.gradient, 1);
obj.auxParameters = obj.projectPoint(obj.auxParameters);
MathUtils.minusEquals(obj.auxParameters,obj.getParameters(),1);
return obj.auxParameters;
}
public boolean stopOptimization(Objective obj){
if(obj instanceof ProjectedObjective) {
ProjectedObjective o = (ProjectedObjective) obj;
double norm = MathUtils.L2Norm(projectGradient(o));
if(norm < gradientConvergenceValue){
// System.out.println("Gradient norm below treshold: " + norm);
return true;
}else{
// System.out.println("projected gradient norm: " + norm);
return false;
}
}
System.out.println("Not a projected objective");
throw new RuntimeException();
}
}
|