An efficient algorithm is presented for numerically evaluating a self-consistent field theoretic (SCFT) model of block copolymer structure. This algorithm is implemented on a distributed memory parallel cluster in order to solve the SCFT equations on large computational grids. Simulation results are presented for a two-component molten mixture of a symmetric ABA triblock copolymer with an A homopolymer. These results illustrate a case in which simulating a large system is required to resolve features with a wide range of length scales. (C) 2003 Elsevier Ltd. All rights reserved.