def test_check_parallel_derivs(self): prob = Problem() root = prob.root = Group() root.add('p1', IndepVarComp('a', 1.0), promotes=['*']) root.add('p2', IndepVarComp('b', 1.0), promotes=['*']) sub1 = root.add('sub1', Group(), promotes=['*']) sub1.ln_solver = LinearGaussSeidel() sub2 = sub1.add('sub2', Group(), promotes=['*']) sub2.add('comp', ExecComp(['x = 2.0*a + 3.0*b', 'y=4.0*a - 1.0*b']), promotes=['*']) sub2.ln_solver = LinearGaussSeidel() root.ln_solver.options['mode'] = 'fwd' sub1.ln_solver.options['mode'] = 'fwd' sub2.ln_solver.options['mode'] = 'fwd' prob.setup(check=False) prob.run() root.ln_solver = LinearGaussSeidel() root.ln_solver.options['single_voi_relevance_reduction'] = True prob.driver.add_desvar('p1.a', 1.0) prob.driver.add_constraint('x', upper=0.0) prob.driver.add_constraint('y', upper=0.0) with warnings.catch_warnings(record=True) as w: if not MPI: # suppress warning about not running under MPI warnings.simplefilter("ignore") prob.driver.parallel_derivs(['x', 'y']) root.ln_solver.options['mode'] = 'rev' sub1.ln_solver.options['mode'] = 'rev' try: mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False) except Exception as err: msg = "Group 'sub2' has mode 'fwd' but the root group has mode 'rev'. Modes must match to use parallel derivative groups." self.assertEqual(text_type(err), msg) else: self.fail('Exception expected') sub1.ln_solver.options['mode'] = 'fwd' sub2.ln_solver.options['mode'] = 'rev' try: mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False) except Exception as err: msg = "Group 'sub1' has mode 'fwd' but the root group has mode 'rev'. Modes must match to use parallel derivative groups." self.assertEqual(text_type(err), msg) else: self.fail('Exception expected') sub1.ln_solver.options['mode'] = 'rev' mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False)
def test_check_parallel_derivs(self): prob = Problem() root = prob.root = Group() root.add('p1', IndepVarComp('a', 1.0), promotes=['*']) root.add('p2', IndepVarComp('b', 1.0), promotes=['*']) sub1 = root.add('sub1', Group(), promotes=['*']) sub1.ln_solver = LinearGaussSeidel() sub2 = sub1.add('sub2', Group(), promotes=['*']) sub2.add('comp', ExecComp(['x = 2.0*a + 3.0*b', 'y=4.0*a - 1.0*b']), promotes=['*']) sub2.ln_solver = LinearGaussSeidel() root.ln_solver.options['mode'] = 'fwd' sub1.ln_solver.options['mode'] = 'fwd' sub2.ln_solver.options['mode'] = 'fwd' prob.setup(check=False) prob.run() root.ln_solver = LinearGaussSeidel() root.ln_solver.options['single_voi_relevance_reduction'] = True prob.driver.add_desvar('p1.a', 1.0) prob.driver.add_constraint('x', upper=0.0) prob.driver.add_constraint('y', upper=0.0) with warnings.catch_warnings(record=True) as w: if not MPI: # suppress warning about not running under MPI warnings.simplefilter("ignore") prob.driver.parallel_derivs(['x','y']) root.ln_solver.options['mode'] = 'rev' sub1.ln_solver.options['mode'] = 'rev' try: mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False) except Exception as err: msg = "Group 'sub2' has mode 'fwd' but the root group has mode 'rev'. Modes must match to use parallel derivative groups." self.assertEqual(text_type(err), msg) else: self.fail('Exception expected') sub1.ln_solver.options['mode'] = 'fwd' sub2.ln_solver.options['mode'] = 'rev' try: mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False) except Exception as err: msg = "Group 'sub1' has mode 'fwd' but the root group has mode 'rev'. Modes must match to use parallel derivative groups." self.assertEqual(text_type(err), msg) else: self.fail('Exception expected') sub1.ln_solver.options['mode'] = 'rev' mode = prob._check_for_parallel_derivs(['a'], ['x'], True, False)
def test_mode_auto(self): # Make sure mode=auto chooses correctly for all prob sizes as well # as for abs/rel/etc paths prob = Problem() root = prob.root = Group() root.add('p1', IndepVarComp('a', 1.0), promotes=['*']) root.add('p2', IndepVarComp('b', 1.0), promotes=['*']) root.add('comp', ExecComp(['x = 2.0*a + 3.0*b', 'y=4.0*a - 1.0*b']), promotes=['*']) root.ln_solver.options['mode'] = 'auto' prob.setup(check=False) prob.run() mode = prob._mode('auto', ['a'], ['x']) self.assertEqual(mode, 'fwd') mode = prob._mode('auto', ['a', 'b'], ['x']) self.assertEqual(mode, 'rev') # make sure _check function does it too #try: #mode = prob._check_for_parallel_derivs(['a'], ['x'], False, False) #except Exception as err: #msg = "Group '' must have the same mode as root to use Matrix Matrix." #self.assertEqual(text_type(err), msg) #else: #self.fail('Exception expected') root.ln_solver.options['mode'] = 'fwd' mode = prob._check_for_parallel_derivs(['a', 'b'], ['x'], False, False) self.assertEqual(mode, 'fwd')
def test_mode_auto(self): # Make sure mode=auto chooses correctly for all prob sizes as well # as for abs/rel/etc paths prob = Problem() root = prob.root = Group() root.add('p1', IndepVarComp('a', 1.0), promotes=['*']) root.add('p2', IndepVarComp('b', 1.0), promotes=['*']) root.add('comp', ExecComp(['x = 2.0*a + 3.0*b', 'y=4.0*a - 1.0*b']), promotes=['*']) root.ln_solver.options['mode'] = 'auto' prob.setup(check=False) prob.run() mode = prob._mode('auto', ['a'], ['x']) self.assertEqual(mode, 'fwd') mode = prob._mode('auto', ['a', 'b'], ['x']) self.assertEqual(mode, 'rev') # make sure _check function does it too #try: #mode = prob._check_for_parallel_derivs(['a'], ['x'], False, False) #except Exception as err: #msg = "Group '' must have the same mode as root to use Matrix Matrix." #self.assertEqual(text_type(err), msg) #else: #self.fail('Exception expected') # Cheat a bit so I can twiddle mode OptionsDictionary.locked = False root.ln_solver.options['mode'] = 'fwd' mode = prob._check_for_parallel_derivs(['a', 'b'], ['x'], False, False) self.assertEqual(mode, 'fwd')