This work was supported by the Applied Mathematical Sciences Program and the HPCC Grand Challenge Program of the U.S. DOE Office of Mathematics;
Information;
and Computational Sciences under Contract DE-AC03-76SF00098. Additional support was provided by the Applied Mathematical Sciences Program of the U.S. DOE Office of Mathematics;
and Computational Sciences under Subcontract 3643100;