/** calc corrected descent direction */
private void mapDirByInverseHessian() {
int m = sList.size();
if (m != 0) {
/*
* V_k-1 = (I - y_k-1 s_k-1^T/(s_k-1 * y_k-1) ?
*
*/
for (int i = m - 1; i >= 0; --i) {
alphas[i] = -ArrayUtils.dot(sList.get(i), dir) / roList.get(i);
addMult(dir, yList.get(i), alphas[i]);
}
/** diagonal elements ? */
double[] lastY = yList.get(m - 1);
double yDotY = ArrayUtils.dot(lastY, lastY);
double scalar = roList.get(m - 1) / yDotY;
ArrayUtils.times(dir, scalar);
/** V_k-1^top ? */
for (int i = 0; i < m; ++i) {
double beta = ArrayUtils.dot(yList.get(i), dir) / roList.get(i);
addMult(dir, sList.get(i), -alphas[i] - beta);
}
}
}
public LinkedList<DirEdge> constructPath(int[] predecessor) {
Stack<Integer> stack = new Stack<Integer>();
int curNode = sink;
// invariant precondition curNode >= 0
do {
stack.push(curNode);
curNode = predecessor[curNode];
} while (curNode != source);
stack.push(source);
LinkedList<DirEdge> resultPath = new LinkedList<DirEdge>();
int size = stack.size();
// invariant precondition i > 0
for (int i = size - 1; i > 0; i--) {
resultPath.add(getEdge(stack.get(i), stack.get(i - 1)));
}
return resultPath;
}
void beginTask(TraceEntry te) {
BdynCallback cb = bdyn_factory.getCallback(te.getEntryLocation());
if (cb == null) return;
if (trace_writer != null)
trace_writer.println(
"PROCESS " + te.getEntryLocation() + " " + te.isExit() + " " + te.getTime());
if (current_transaction == null) {
OutputTask ot = null;
if (te.isExit()) return;
int i0 = te.getObject1();
int i1 = te.getObject2();
if (i0 != 0) {
ot = object_tasks.get(i0);
if (ot == dummy_task) ot = null;
}
if (ot == null) {
if (i1 != 0) {
ot = object_tasks.get(i1);
if (ot == dummy_task) ot = null;
}
}
if (ot == null) {
ot = new OutputTask(++task_counter, cb);
}
current_transaction = ot;
nest_level = 0;
}
if (!te.isExit()) {
OutputEntry oe = null;
if (nest_level == 0) {
oe = createOutputEntry(te);
} else if (cb.getCallbackType() == CallbackType.KEY
&& BdynFactory.getOptions().useKeyCallback()) {
for (int i = output_stack.size() - 1; i >= 0; --i) {
OutputEntry poe = output_stack.get(i);
if (poe != null) {
poe.finishAt(te.getTime());
if (trace_writer != null)
trace_writer.println(
"FINISHINT "
+ poe.getEntryTask().getId()
+ " "
+ nest_level
+ " "
+ poe.getTotalTime(0)
+ " "
+ poe.hashCode());
break;
}
}
oe = createOutputEntry(te);
} else {
if (trace_writer != null) trace_writer.println("IGNORE " + cb.getId() + " " + nest_level);
// System.err.println("TRACE: Ignore " + cb.getDisplayName() + " " + nest_level);
}
output_stack.push(oe);
active_threads.add(this);
++nest_level;
} else {
if (nest_level <= 0) {
return;
}
--nest_level;
OutputEntry oe = output_stack.pop();
if (oe != null) {
// TODO: max_delta doesn't take into account calls that haven't ended yet
max_delta = Math.max(max_delta, te.getTime() - oe.getStartTime());
oe.finishAt(te.getTime());
if (trace_writer != null)
trace_writer.println(
"FINISH "
+ oe.getEntryTask().getId()
+ " "
+ nest_level
+ " "
+ oe.getTotalTime(0)
+ " "
+ oe.hashCode());
if (!oe.isRelevant()) {
output_set.remove(oe);
if (trace_writer != null)
trace_writer.println("REMOVE KEY " + oe.getEntryTask().getId() + " " + nest_level);
// System.err.println("REMOVE KEY " + oe.getEntryTask().getDisplayName());
oe = null;
for (int i = output_stack.size() - 1; i >= 0; --i) {
OutputEntry poe = output_stack.get(i);
if (poe != null) {
poe.finishAt(0);
if (trace_writer != null)
trace_writer.println(
"NO END "
+ poe.getEntryTask().getId()
+ " "
+ nest_level
+ " "
+ poe.hashCode());
break;
}
}
}
}
if (last_entry != null
&& oe != null
&& last_entry.getEntryTask() == oe.getEntryTask()
&& (te.getTime() - last_entry.getStartTime()) < MERGE_TIME) {
last_entry.mergeWith(oe);
output_set.remove(oe);
if (trace_writer != null)
trace_writer.println("MERGE KEY " + oe.getEntryTask().getId() + " " + nest_level);
// System.err.println("REMOVE " + oe.getEntryTask().getDisplayName());
}
if (oe != null) last_entry = oe;
if (nest_level == 0) endTask();
else if (oe != null) {
for (int i = output_stack.size() - 1; i >= 0; --i) {
OutputEntry poe = output_stack.get(i);
if (poe != null) {
if (!poe.isSignificant()) {
output_set.remove(poe);
}
if (trace_writer != null)
trace_writer.println(
"NEST KEY "
+ poe.getEntryTask().getId()
+ " "
+ nest_level
+ " "
+ poe.isSignificant());
oe = new OutputEntry(te.getTime(), 0, this, current_transaction, poe.getEntryTask());
oe.setDeletable();
output_set.add(oe);
output_stack.set(i, oe);
break;
}
}
}
}
}
