/** 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; } } } } }