1 #ifndef CPPAD_CG_LANGUAGE_LATEX_INDEX_PATTERNS_INCLUDED 2 #define CPPAD_CG_LANGUAGE_LATEX_INDEX_PATTERNS_INCLUDED 22 inline void LanguageLatex<Base>::generateNames4RandomIndexPatterns(
const std::set<RandomIndexPattern*>& randomPatterns) {
23 std::ostringstream os;
25 std::set<std::string> usedNames;
29 for (RandomIndexPattern* ip : randomPatterns) {
30 if (!ip->getName().empty()) {
31 usedNames.insert(ip->getName());
37 for (RandomIndexPattern* ip : randomPatterns) {
38 if (ip->getName().empty()) {
42 os << _C_STATIC_INDEX_ARRAY << c;
46 }
while (usedNames.find(name) != usedNames.end());
56 const std::string& indentation,
57 const std::set<RandomIndexPattern*>& randomPatterns) {
59 if (ip->getType() == IndexPatternType::Random1D) {
64 const std::map<size_t, size_t>& x2y = ip1->getValues();
66 std::vector<size_t> y(x2y.rbegin()->first + 1);
67 for (
const std::pair<size_t, size_t>& p : x2y)
68 y[p.first] = p.second;
71 printStaticIndexArray(os, ip->getName(), y);
73 CPPADCG_ASSERT_UNKNOWN(ip->getType() == IndexPatternType::Random2D);
79 printStaticIndexMatrix(os, ip->getName(), ip2->getValues());
86 const std::string& name,
87 const std::vector<size_t>& values) {
88 os << name <<
" = \\left\\{";
89 if (!values.empty()) {
91 for (
size_t i = 1; i < values.size(); i++) {
92 os <<
", " << values[i];
95 os <<
"\\right\\}" << _endEq <<
" % size: " << values.size() << _endline;
100 const std::string& name,
101 const std::map<
size_t, std::map<size_t, size_t> >& values) {
105 std::map<size_t, std::map<size_t, size_t> >::const_iterator it;
106 std::map<size_t, size_t>::const_iterator ity2z;
108 if (!values.empty()) {
109 m = values.rbegin()->first + 1;
111 for (it = values.begin(); it != values.end(); ++it) {
112 if (!it->second.empty())
113 n = std::max(n, it->second.rbegin()->first + 1);
117 os << name <<
" = \\left\\{";
119 for (it = values.begin(); it != values.end(); ++it) {
120 if (it->first != x) {
121 while (it->first != x) {
130 for (ity2z = it->second.begin(); ity2z != it->second.end(); ++ity2z) {
131 if (ity2z->first != y) {
132 while (ity2z->first != y) {
139 if (ity2z->first != it->second.rbegin()->first) os <<
",";
145 if (it->first != values.rbegin()->first) os <<
",";
149 os <<
"\\right\\}" << _endEq <<
"% size: " << m <<
" x " << n << _endline;
155 return indexPattern2String(ip,{&index});
161 std::stringstream ss;
162 switch (ip.getType()) {
163 case IndexPatternType::Linear:
165 CPPADCG_ASSERT_KNOWN(indexes.size() == 1,
"Invalid number of indexes");
167 return linearIndexPattern2String(lip, *indexes[0]);
169 case IndexPatternType::Sectioned:
171 CPPADCG_ASSERT_KNOWN(indexes.size() == 1,
"Invalid number of indexes");
173 const std::map<size_t, IndexPattern*>& sections = lip->getLinearSections();
174 size_t sSize = sections.size();
175 CPPADCG_ASSERT_UNKNOWN(sSize > 1);
177 std::map<size_t, IndexPattern*>::const_iterator its = sections.begin();
178 for (
size_t s = 0; s < sSize - 1; s++) {
181 size_t xStart = its->first;
183 ss <<
"\\left(" << (*indexes[0]->getName()) <<
"<" << xStart <<
"\\right)? " 184 << indexPattern2String(*lp, *indexes[0]) <<
": ";
186 ss << indexPattern2String(*its->second, *indexes[0]);
191 case IndexPatternType::Plane2D:
193 CPPADCG_ASSERT_KNOWN(indexes.size() >= 1,
"Invalid number of indexes");
194 std::string indexExpr;
196 bool useParens = pip.getPattern1() !=
nullptr && pip.getPattern2() !=
nullptr;
198 if (useParens) indexExpr +=
"\\left(";
200 if (pip.getPattern1() !=
nullptr)
201 indexExpr += indexPattern2String(*pip.getPattern1(), *indexes[0]);
203 if (useParens) indexExpr +=
"\\right) + \\left(";
205 if (pip.getPattern2() !=
nullptr)
206 indexExpr += indexPattern2String(*pip.getPattern2(), *indexes.back());
208 if (useParens) indexExpr +=
"\\right)";
212 case IndexPatternType::Random1D:
214 CPPADCG_ASSERT_KNOWN(indexes.size() == 1,
"Invalid number of indexes");
216 CPPADCG_ASSERT_KNOWN(!rip.getName().empty(),
"Invalid name for array");
217 return rip.getName() +
"\\left[" + (*indexes[0]->getName()) +
"\\right]";
219 case IndexPatternType::Random2D:
221 CPPADCG_ASSERT_KNOWN(indexes.size() == 2,
"Invalid number of indexes");
223 CPPADCG_ASSERT_KNOWN(!rip.getName().empty(),
"Invalid name for array");
224 return rip.getName() +
"\\left[" + (*indexes[0]->getName()) +
"\\right]\\left[" + (*indexes[1]->getName()) +
"\\right]";
227 CPPADCG_ASSERT_UNKNOWN(
false);
235 long dy = lip.getLinearSlopeDy();
236 long dx = lip.getLinearSlopeDx();
237 long b = lip.getLinearConstantTerm();
238 long xOffset = lip.getXOffset();
240 std::stringstream ss;
247 ss <<
" - " << xOffset <<
"\\right)";
254 ss <<
" \\cdot " << dy;
const std::string * getName() const
void printRandomIndexPatternDeclaration(std::ostringstream &os, const std::string &identation, const std::set< RandomIndexPattern *> &randomPatterns)