Invited talk at Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences