BOOST_FORCEINLINE result_type operator()( A0& a0, A1& s, A2& u,A3 &vt, A4 const jobu , A5 const jobvt) const { result_type that; details::workspace<typename A0::value_type> w; nt2_la_int m = nt2::height(a0); nt2_la_int n = nt2::width(a0); nt2_la_int ld = a0.leading_size(); nt2_la_int ldu = u.leading_size() > 1 ? u.leading_size() : 1 ; nt2_la_int ldvt= vt.leading_size() > 1 ? vt.leading_size() : 1 ; nt2_la_int lwork_query = -1; char c = std::toupper(jobu); magma_vec_t magma_jobu = c == 'A' ? MagmaAllVec : c == 'S' ? MagmaSomeVec : c == 'O' ? MagmaOverwriteVec : c == 'N' ? MagmaNoVec : (magma_vec_t)0; c = std::toupper(jobvt); magma_vec_t magma_jobvt = c == 'A' ? MagmaAllVec : c == 'S' ? MagmaSomeVec : c == 'O' ? MagmaOverwriteVec : c == 'N' ? MagmaNoVec : (magma_vec_t)0; magma_sgesvd(magma_jobu,magma_jobvt,m, n, 0, ld, 0, 0, ldu , 0, ldvt, w.main() , lwork_query, &that ); w.prepare_main(); nt2::gesvd(a0,s,u,vt,jobu,jobvt,w); return that; }
BOOST_FORCEINLINE result_type operator()(A0& a0, A1& s, A2& u,A3 &vt, A4 jobu, A5 jobvt) const { result_type that; details::workspace<typename A0::value_type> w; nt2_la_int m = nt2::height(a0); nt2_la_int n = nt2::width(a0); nt2_la_int ld = a0.leading_size(); nt2_la_int ldu = u.leading_size() > 1 ? u.leading_size() : 1 ; nt2_la_int ldvt= vt.leading_size() > 1 ? vt.leading_size() : 1 ; NT2_F77NAME(sgesvd) (&jobu,&jobvt,&m, &n, 0, &ld, 0, 0, &ldu , 0, &ldvt, w.main() , details::query(), &that ); w.resize_main(5*nt2::max(m, n)); nt2::gesvd(a0,s,u,vt,jobu,jobvt,w); return that; }
void Print() { A1::Print(); A2::Print(); aa3.Print(); cout<<"d="<<d<<endl; }
int CmdAppend::execute (std::string& output) { int rc = 0; int count = 0; // Apply filter. std::vector <Task> filtered; filter (filtered); if (filtered.size () == 0) { context.footnote (STRING_FEEDBACK_NO_TASKS_SP); return 1; } // Apply the command line modifications to the new task. A3 modifications = context.a3.extract_modifications (); if (!modifications.size ()) throw std::string (STRING_CMD_MODIFY_NEED_TEXT); // Accumulated project change notifications. std::map <std::string, std::string> projectChanges; std::vector <Task>::iterator task; for (task = filtered.begin (); task != filtered.end (); ++task) { Task before (*task); // Append to the specified task. std::string question = format (STRING_CMD_APPEND_CONFIRM, task->id, task->get ("description")); modify_task_description_append (*task, modifications); if (permission (*task, taskDifferences (before, *task) + question, filtered.size ())) { context.tdb2.modify (*task); ++count; feedback_affected (STRING_CMD_APPEND_TASK, *task); if (context.verbose ("project")) projectChanges[task->get ("project")] = onProjectChange (*task, false); // Append to siblings. if (task->has ("parent")) { std::vector <Task> siblings = context.tdb2.siblings (*task); if (siblings.size () && confirm (STRING_CMD_APPEND_CONFIRM_R)) { std::vector <Task>::iterator sibling; for (sibling = siblings.begin (); sibling != siblings.end (); ++sibling) { modify_task_description_append (*sibling, modifications); context.tdb2.modify (*sibling); ++count; feedback_affected (STRING_CMD_APPEND_TASK_R, *sibling); } // Append to the parent Task parent; context.tdb2.get (task->get ("parent"), parent); modify_task_description_append (parent, modifications); context.tdb2.modify (parent); } } } else { std::cout << STRING_CMD_APPEND_NO << "\n"; rc = 1; if (_permission_quit) break; } } // Now list the project changes. std::map <std::string, std::string>::iterator i; for (i = projectChanges.begin (); i != projectChanges.end (); ++i) if (i->first != "") context.footnote (i->second); context.tdb2.commit (); feedback_affected (count == 1 ? STRING_CMD_APPEND_1 : STRING_CMD_APPEND_N, count); return rc; }
int CmdModify::execute (std::string& output) { int rc = 0; int count = 0; // Apply filter. std::vector <Task> filtered; filter (filtered); if (filtered.size () == 0) { context.footnote (STRING_FEEDBACK_NO_TASKS_SP); return 1; } // Apply the command line modifications to the new task. A3 modifications = context.a3.extract_modifications (); if (!modifications.size ()) throw std::string (STRING_CMD_MODIFY_NEED_TEXT); // Accumulated project change notifications. std::map <std::string, std::string> projectChanges; std::vector <Task>::iterator task; for (task = filtered.begin (); task != filtered.end (); ++task) { Task before (*task); modify_task_description_replace (*task, modifications); if (taskDiff (before, *task)) { // Perform some logical consistency checks. if (task->has ("recur") && !task->has ("due") && !before.has ("due")) throw std::string (STRING_CMD_MODIFY_NO_DUE); if (before.has ("recur") && before.has ("due") && (!task->has ("due") || task->get ("due") == "")) throw std::string (STRING_CMD_MODIFY_REM_DUE); if (before.has ("recur") && (!task->has ("recur") || task->get ("recur") == "")) throw std::string (STRING_CMD_MODIFY_REC_ALWAYS); // Delete the specified task. std::string question = format (STRING_CMD_MODIFY_CONFIRM, task->id, task->get ("description")); if (permission (*task, taskDifferences (before, *task) + question, filtered.size ())) { updateRecurrenceMask (*task); dependencyChainOnModify (before, *task); ++count; feedback_affected (STRING_CMD_MODIFY_TASK, *task); feedback_unblocked (*task); context.tdb2.modify (*task); if (context.verbose ("project")) projectChanges[task->get ("project")] = onProjectChange (before, *task); // Task potentially has siblings - modify them. if (task->has ("parent")) { std::vector <Task> siblings = context.tdb2.siblings (*task); if (siblings.size () && confirm (STRING_CMD_MODIFY_RECUR)) { std::vector <Task>::iterator sibling; for (sibling = siblings.begin (); sibling != siblings.end (); ++sibling) { Task alternate (*sibling); modify_task_description_replace (*sibling, modifications); updateRecurrenceMask (*sibling); dependencyChainOnModify (alternate, *sibling); ++count; feedback_affected (STRING_CMD_MODIFY_TASK_R, *sibling); feedback_unblocked (*sibling); context.tdb2.modify (*sibling); if (context.verbose ("project")) projectChanges[sibling->get ("project")] = onProjectChange (alternate, *sibling); } } } // Task potentially has child tasks - modify them. else if (task->get ("status") == "recurring") { std::vector <Task> children = context.tdb2.children (*task); if (children.size () && confirm (STRING_CMD_MODIFY_RECUR)) { std::vector <Task>::iterator child; for (child = children.begin (); child != children.end (); ++child) { Task alternate (*child); modify_task_description_replace (*child, modifications); updateRecurrenceMask (*child); context.tdb2.modify (*child); dependencyChainOnModify (alternate, *child); if (context.verbose ("project")) projectChanges[child->get ("project")] = onProjectChange (alternate, *child); ++count; feedback_affected (STRING_CMD_MODIFY_TASK_R, *child); } } } } else { std::cout << STRING_CMD_MODIFY_NO << "\n"; rc = 1; } } } // Now list the project changes. std::map <std::string, std::string>::iterator i; for (i = projectChanges.begin (); i != projectChanges.end (); ++i) if (i->first != "") context.footnote (i->second); context.tdb2.commit (); feedback_affected (count == 1 ? STRING_CMD_MODIFY_1 : STRING_CMD_MODIFY_N, count); return rc; }