DGW: при сложении или вычитании матриц не квадратного порядка вылетает сообщение с ошибкой

private void button1_Click(object sender, EventArgs e)
        {
            double[,] a, b, c;
            int cola = int.Parse(textBox1.Text);
            int rowa = int.Parse(textBox2.Text);
            int colb = int.Parse(textBox3.Text);
            int rowb = int.Parse(textBox4.Text);
            a = new double[cola, rowa];
            b = new double[colb, rowb];
            c = new double[cola, rowa];
            string elem, elem1;
 
            for (int i = 0; i < cola; i++)
                for (int j = 0; j < rowa; j++)
                {
                    elem = dGV1.Rows[i].Cells[j].Value.ToString();
                    a[i, j] = Convert.ToDouble(elem);
                }
 
            for (int i = 0; i < colb; i++)
                for (int j = 0; j < rowb; j++)
                    {
                        elem1 = dGV2.Rows[i].Cells[j].Value.ToString();
                        b[i, j] = Convert.ToDouble(elem1);
 
                    }
 
            switch (comboBox1.SelectedIndex)
            {
                case 0:
                    dGV3.RowCount = rowa;
                    dGV3.ColumnCount = cola;
                    for (int i = 0; i < cola; i++)
                        for (int j = 0; j < rowa; j++)
                        {
                            c[i, j] = a[i, j] + b[i, j];
                            dGV3.Rows[i].Cells[j].Value = c[i, j];
                        }
                    break;
                case 1:
                    dGV3.RowCount = rowa;
                    dGV3.ColumnCount = cola;
                    for (int i = 0; i < cola; i++)
                        for (int j = 0; j < rowa; j++)
                        {
                            c[i, j] = a[i, j] - b[i, j];
                            dGV3.Rows[i].Cells[j].Value = c[i, j];
                        }
                    break;
            }
        }