mosekcomodel 0.5.0

Library for Conic Optimization Modeling with Mosek
Documentation 
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131


use std::iter::StepBy;
use std::ops::{Range,RangeTo,RangeFrom,RangeFull, RangeToInclusive,RangeInclusive};
use super::variable::*;
use super::ModelItemIndex;

pub enum IndexType {
    Index(usize),
    Range(usize,usize),
    RangeFrom(usize),
    RangeFull(),
}



trait ModelItemIndexElement {
    fn tt(self) -> IndexType;
}

impl ModelItemIndexElement for usize  { 
    fn tt(self) -> IndexType {
        IndexType::Index(self)
    }
}

impl ModelItemIndexElement for Range<usize> { 
    fn tt(self) -> IndexType {
        if self.start <= self.end {
            IndexType::Range(self.start,self.end-self.start)
        }
        else {
            IndexType::RevRange(self.start,0,1)
        }
    }
}

impl ModelItemIndexElement for RangeInclusive<usize> { 
    fn tt(self) -> IndexType {
        if self.start() <= self.end() {
            IndexType::Range(*self.start(), self.end()-self.start()+1)
        }
        else {
            IndexType::Range(self.start())
        }
    }
}

impl ModelItemIndexElement for RangeFrom<usize> { 
    fn tt(self) -> IndexType {
        IndexType::RangeFrom(self.start)
    }
}

impl ModelItemIndexElement for RangeTo<usize> { 
    fn tt(self) -> IndexType {
        IndexType::Range(0,self.end)
    }
}

impl ModelItemIndexElement for RangeToInclusive<usize> { 
    fn tt(self) -> IndexType {
        IndexType::Range(0,self.end+1,1)
    }
}

impl ModelItemIndexElement for RangeFull { 
    fn tt(self) -> IndexType {
        IndexType::RangeFull(1)
    }
}

impl<I> ModelItemIndexElement for I where I : Iterator<Item = usize> {
    fn tt(self) -> IndexType {
        
    }
}









impl<const N : usize> ModelItemIndex<Variable<N>> for [IndexType; N] {
    type Output = Variable<N>;

    fn index(self,obj : &Variable<N>) -> Self::Output {
        unimplemented!("Index not implemented");
    }
}




impl<I1,I2> ModelItemIndex<Variable<2>> for (I1,I2) 
    where
        I1 : ModelItemIndexElement,
        I2 : ModelItemIndexElement
{
    type Output = Variable<2>;
    fn index(self, v : &Variable<1>) -> Variable<1> {
        let n = self.len();
        if let Some(ref sp) = v.sparsity {
            let first = match sp.binary_search(&self.start) {
                Ok(i)  => i,
                Err(i) => i
            };
            let last = match sp.binary_search(&self.start) {
                Ok(i) => i+1,
                Err(i) => i
            };

            Variable{
                idxs     : Rc::new(v.idxs[first..last].to_vec()), 
                sparsity : Some(Rc::new(sp[first..last].iter().map(|&i| i - self.start).collect())),
                shape    : [n]
            }
        }
        else {
            Variable{
                idxs     : Rc::new(v.idxs[self].to_vec()),
                sparsity : None,
                shape    : [n]
            }
        }+1
    }
}