280 do 100 mrs = 1, nrs
281 C
282 do 10 mq = 1, nq
283 do 10 mi = 1, morb
284 10 xrsiq(mi, mq) = zero
285 do 40 mp = 1, np
286 do 30 mq = 1, mp
287 mrspq = mrspq+1
288 val = xrspq(mrspq)
289 if (val .eq. zero)
289 goto 30
290 do 20 mi = 1, morb
291 xrsiq(mi, mq) = xrsiq(mi, mq)+val*v(mp, mi)
292 xrsiq(mi, mp) = xrsiq(mi, mp)+val*v(mq, mi)
293 20 continue
294 30 continue
295 40 continue
296 C
297 mrsij = mrsij0
298 do 90 mi = 1, morb
299 do 50 mj = 1, mi
300 50 xij(mj) = zero
301 do 70 mq = 1, nq
302 val = xrsiq(mi, mq)
303 if (val .eq. zero)
303 goto 70
304 do 60 mj = 1, mi
305 60 xij(mj) = xij(mj)+val*v(mq, mj)
306 70 continue
307 do 80 mj = 1, mi
308 mrsij = mrsij+1
309 80 xrsij(mrsij) = xij(mj)
310 90 continue
311 C
312 mrsij0 = mrsij0+nrs
313 100 continue
![]()