Partial vector substitution

[bqpd]

allows substitutions like m.substitutions[v[0,2:]] = 1 # or np.ones(v[0,2:].shape)

[bqpd]

test models please

[bqpd]

@pgkirsch what do you think of this feature's utility? It seems to me like something that should Just Work, so I'm willing to take the performance hit...but that's easy for me to say.

[pgkirsch]

@bqpd that seems helpful! This seems like it would further simplify the proposed aero changes?

I assume that syntax carries over to the setup method form?

def setup(self):
...
    substitutions = {
        v[0,2:]: 1,
    ...
    }

    return constraints, substitutions

I know I'm sort of pushing the envelope on some of these vector variable subs, but do you have other users with anything similar in use case?

[bqpd]

yup, carries over!

I think this is one of those things where people might have tried to use it, it didn't work, they found an alternative without mentioning...

[pgkirsch]

Awesome!

Have you profiled what the performance hit is for a large model? I can't imagine it's too significant, right?

[bqpd]

it was ~4% of the time in that function...don't remember exactly the timing in my test run, but no more than a couple hundredths of a second.

…

[pgkirsch]

Seems very tolerable!

[bqpd]

I'm just so used to making things quicker when I make changes to core classes :p

[pgkirsch]

Ha sorry!

[bqpd]

test models please