2017-09-08 Nicolas Roche <roche@adacore.com>

	* gcc-interface/Make-lang.in, gcc-interface/Makefile.in: Find runtime
	source in libgnat/
	* a-lfztio.ads, g-timsta.ads, g-sercom-linux.adb, s-osprim-solaris.adb,
	a-inteio.ads, s-stchop-rtems.adb, s-casuti.adb, s-pack39.adb,
	i-vxwork-x86.ads, a-strbou.adb, a-stzmap.adb, s-assert.adb,
	a-sfecin.ads, a-cohama.adb, s-casuti.ads, a-suenco.adb, s-pack39.ads,
	a-stzmap.ads, a-strbou.ads, s-stalib.adb, s-trasym.adb, g-comver.adb,
	s-assert.ads, s-vector.ads, g-cgi.adb, a-cohama.ads, s-wchcnv.adb,
	a-titest.adb, s-pack48.adb, a-suenco.ads, a-strunb.adb, s-stalib.ads,
	s-trasym.ads, a-nudira.adb, g-comver.ads, a-nuflra.adb, g-cgi.ads,
	a-chacon.adb, s-wchcnv.ads, a-excach.adb, s-pack48.ads, a-titest.ads,
	a-strunb.ads, s-dwalin.adb, a-nudira.ads, a-chtgbo.adb, s-resfil.adb,
	a-scteio.ads, a-nuflra.ads, g-soliop-mingw.ads, s-pack57.adb,
	a-chacon.ads, s-bytswa.ads, s-pooloc.adb, g-os_lib.adb, s-dwalin.ads,
	a-szuzha.adb, s-resfil.ads, a-chtgbo.ads, s-spsufi.adb, s-pack57.ads,
	s-pooloc.ads, g-os_lib.ads, a-stfiha.ads, a-lcteio.ads, a-wtcoau.adb,
	a-szuzha.ads, s-mmosin-unix.adb, a-stmaco.ads, s-spsufi.ads,
	s-stchop-limit.ads, a-wtcoau.ads, a-exctra.adb, s-mmosin-unix.ads,
	s-sequio.adb, s-conca2.adb, g-table.adb, s-imglli.adb,
	a-numaux-x86.adb, a-strsea.adb, s-wchstw.adb, a-clrefi.adb,
	a-wwboio.adb, a-exctra.ads, s-sequio.ads, s-conca2.ads, a-wwunio.ads,
	system-linux-hppa.ads, g-table.ads, s-dimkio.ads, s-imglli.ads,
	a-cofove.adb, a-numaux-x86.ads, s-wchstw.ads, a-strsea.ads,
	a-clrefi.ads, a-wwboio.ads, s-stratt-xdr.adb, s-crc32.adb,
	s-excmac-arm.adb, g-busora.adb, a-cofove.ads, s-osprim-unix.adb,
	g-io.adb, s-pack49.adb, s-crc32.ads, s-excmac-arm.ads, a-fzteio.ads,
	g-busora.ads, s-stausa.adb, system-linux-mips.ads, sequenio.ads,
	g-exctra.adb, g-rewdat.adb, a-cgaaso.adb, g-io.ads, s-pack49.ads,
	a-wtflau.adb, a-undesu.adb, s-stausa.ads, a-ztenau.adb, g-enutst.ads,
	calendar.ads, s-pack58.adb, g-rewdat.ads, g-exctra.ads, s-ststop.adb,
	a-cgaaso.ads, a-strfix.adb, a-comlin.adb, a-strunb-shared.adb,
	a-wtflau.ads, a-undesu.ads, a-cbhase.adb, a-ztenau.ads, s-os_lib.adb,
	a-coorse.adb, a-chlat1.ads, s-pack58.ads, s-ststop.ads, a-strfix.ads,
	a-comlin.ads, a-strunb-shared.ads, a-nscefu.ads, s-valboo.adb,
	directio.ads, a-chtgke.adb, a-cbhase.ads, a-wtinau.adb,
	system-linux-alpha.ads, s-os_lib.ads, a-coorse.ads,
	system-linux-s390.ads, s-imgwiu.adb, a-chtgop.adb, s-valboo.ads,
	a-chtgke.ads, a-tienio.adb, s-conca3.adb, a-wtinau.ads,
	system-darwin-ppc.ads, i-c.adb, s-expllu.adb, g-expect.adb,
	g-sha256.ads, s-vallld.adb, s-imgwiu.ads, a-chtgop.ads, a-strmap.adb,
	a-tienio.ads, s-conca3.ads, s-imgint.adb, i-c.ads, s-expllu.ads,
	s-osprim-darwin.adb, a-cogeso.adb, g-expect.ads, a-iwteio.ads,
	s-vallld.ads, a-coinho-shared.adb, g-shsh64.adb, a-strmap.ads,
	g-comlin.adb, a-excpol.adb, s-imgint.ads, a-ztdeau.adb, a-cogeso.ads,
	a-coinho-shared.ads, g-shsh64.ads, g-comlin.ads, a-stzsup.adb,
	a-rbtgbk.adb, a-wtmoau.adb, a-ztdeau.ads, s-exnlli.adb, g-tty.adb,
	g-heasor.adb, g-socthi-dummy.adb, s-llflex.ads, a-zchara.ads,
	a-stzsup.ads, a-ztcstr.adb, a-rbtgbk.ads, a-sfwtio.ads, a-wtmoau.ads,
	a-sulcin.adb, s-exnlli.ads, system-freebsd.ads, a-stunha.adb,
	a-charac.ads, g-tty.ads, g-heasor.ads, s-exctra.adb,
	g-socthi-dummy.ads, a-coboho.adb, a-ztcstr.ads, a-tideio.adb,
	a-sulcin.ads, a-wrstfi.adb, g-alleve.adb, s-pack59.adb, a-ngrear.adb,
	a-stboha.adb, a-stunau-shared.adb, a-stunha.ads, a-lfwtio.ads,
	s-fileio.adb, s-exctra.ads, a-coboho.ads, a-ioexce.ads, a-tideio.ads,
	a-ngrear.ads, a-wrstfi.ads, s-pack59.ads, g-alleve.ads, a-stboha.ads,
	s-poosiz.adb, g-traceb.adb, g-rannum.adb, machcode.ads, s-purexc.ads,
	s-fileio.ads, a-cfinve.adb, a-crbtgk.adb, system-solaris-x86.ads,
	s-poosiz.ads, g-rannum.ads, g-traceb.ads, a-except.adb, s-conca4.adb,
	a-stream.adb, a-cfinve.ads, a-crbtgk.ads, s-wchwts.adb,
	system-mingw.ads, a-except.ads, s-conca4.ads, a-chzla9.ads,
	s-valenu.adb, s-soflin.adb, a-stream.ads, a-cgarso.adb, s-valllu.adb,
	g-crc32.adb, s-wchwts.ads, s-fatflt.ads, s-imguns.adb, s-strcom.adb,
	g-decstr.adb, s-valenu.ads, s-soflin.ads, a-cgarso.ads, a-cwila1.ads,
	s-valllu.ads, g-crc32.ads, s-imguns.ads, g-spipat.adb, s-valwch.adb,
	s-strcom.ads, g-decstr.ads, text_io.ads, g-debuti.adb, s-stchop.adb,
	g-spipat.ads, s-valwch.ads, a-string.ads, s-exnint.adb, g-awk.adb,
	g-tasloc.adb, s-wwdenu.adb, s-boustr.adb, a-zchuni.adb, s-stchop.ads,
	g-debuti.ads, s-stopoo.adb, system-dragonfly-x86_64.ads,
	system-linux-x86.ads, s-exnint.ads, g-awk.ads, a-stzhas.adb,
	g-tasloc.ads, s-wwdenu.ads, g-debpoo.adb, g-except.ads,
	g-sse.ads, s-boustr.ads, a-zchuni.ads, s-bitops.adb, s-wwdwch.adb,
	s-stopoo.ads, a-catizo.adb, a-stzhas.ads, a-nlcefu.ads, g-debpoo.ads,
	i-vxwoio.adb, s-bitops.ads, g-io-put-vxworks.adb, s-wwdwch.ads,
	g-sehamd.adb, a-ssicst.adb, a-catizo.ads, s-mmap.adb, g-string.adb,
	s-traceb.adb, a-swunau.adb, s-rannum.adb, a-ticoau.adb, i-vxwoio.ads,
	g-sehamd.ads, a-stwiun.adb, a-ssicst.ads, s-conca5.adb, a-ssitio.ads,
	s-mmap.ads, a-zttest.adb, g-string.ads, g-sercom.adb, a-cdlili.adb,
	a-swunau.ads, s-traceb.ads, s-rannum.ads, a-ticoau.ads, system-aix.ads,
	a-cforma.adb, a-stwiun.ads, s-conca5.ads, s-carsi8.adb, a-zttest.ads,
	g-sercom.ads, a-cdlili.ads, a-cihama.adb, g-sptain.ads, a-cforma.ads,
	s-maccod.ads, s-carsi8.ads, a-strsup.adb, g-sha1.adb, a-cihama.ads,
	g-stseme.adb, s-traent.adb, s-valcha.adb, g-curexc.ads, a-strsup.ads,
	g-sha1.ads, a-sflcin.ads, s-traent.ads, s-pack10.adb, s-valcha.ads,
	a-coteio.ads, s-tasloc.adb, g-utf_32.adb, a-suteio.adb, s-except.adb,
	a-direct.adb, g-stsifd-sockets.adb, a-numaux-vxworks.ads, s-winext.ads,
	s-pack10.ads, a-ztexio.adb, a-tiflau.adb, system-vxworks-arm.ads,
	s-tasloc.ads, a-suteio.ads, g-utf_32.ads, s-except.ads,
	a-direct.ads, a-swbwha.adb, g-hesorg.adb, s-wwdcha.adb, a-wtedit.adb,
	a-ztexio.ads, a-wtcoio.adb, a-tiflau.ads, a-ssizti.ads, s-casi32.adb,
	a-swbwha.ads, s-veboop.adb, g-hesorg.ads, s-parame-rtems.adb,
	s-wwdcha.ads, a-wtedit.ads, a-stuten.adb, a-coinve.adb, a-wtcoio.ads,
	s-casi32.ads, s-string.adb, a-tiinau.adb, a-cusyqu.adb, s-conca6.adb,
	s-veboop.ads, a-cgcaso.adb, a-numaux-darwin.adb, a-envvar.adb,
	a-stuten.ads, s-secsta.adb, a-coinve.ads, s-string.ads, a-cusyqu.ads,
	a-tiinau.ads, s-osprim-vxworks.adb, s-conca6.ads, g-spchge.adb,
	s-parint.adb, a-cuprqu.adb, a-cgcaso.ads, a-numaux-darwin.ads,
	a-envvar.ads, s-secsta.ads, g-spchge.ads, s-parint.ads, a-cuprqu.ads,
	a-swuwti.adb, a-flteio.ads, a-sbhcin.adb, a-coprnu.adb, g-u3spch.adb,
	s-atocou.adb, g-ctrl_c.adb, a-swuwti.ads, a-calend.adb, a-sbhcin.ads,
	a-coprnu.ads, g-dirope.adb, g-sha512.ads, g-u3spch.ads, s-atocou.ads,
	g-ctrl_c.ads, a-timoau.adb, a-witeio.adb, s-pack11.adb, a-strhas.adb,
	a-wtflio.adb, g-spitbo.adb, a-calend.ads, a-ztenio.adb, g-dirope.ads,
	a-slcain.adb, g-sechas.adb, a-timoau.ads, a-witeio.ads, s-pack11.ads,
	s-shasto.adb, s-traceb-mastop.adb, a-ciorse.adb, s-utf_32.adb,
	a-strhas.ads, a-wtflio.ads, g-spitbo.ads, a-ztenio.ads, a-slcain.ads,
	g-sechas.ads, s-gearop.adb, a-siztio.ads, s-pack20.adb, s-shasto.ads,
	a-ciorse.ads, s-utf_32.ads, s-crtl.ads, a-wtinio.adb, s-elaall.adb,
	s-explli.adb, s-chepoo.ads, s-gearop.ads, a-einuoc.adb, s-pack20.ads,
	system-linux-ia64.ads, a-swunau-shared.adb, a-wtinio.ads, g-alvety.ads,
	a-liztio.ads, g-calend.adb, s-conca7.adb, s-elaall.ads, s-explli.ads,
	a-einuoc.ads, s-widboo.adb, s-imgdec.adb, a-cbhama.adb, g-calend.ads,
	s-conca7.ads, a-llitio.ads, i-cexten.ads, a-coorma.adb, s-widboo.ads,
	s-diflio.adb, g-souinf.ads, s-imgdec.ads, g-strhas.ads, a-cbhama.ads,
	g-shshco.adb, a-ztdeio.adb, s-gloloc.adb, a-coorma.ads, g-wispch.adb,
	s-pack03.adb, g-eacodu.adb, s-casi16.adb, s-diflio.ads, a-colien.adb,
	g-shshco.ads, a-wtmoio.adb, a-rbtgbo.adb, a-ztdeio.ads,
	system-rtems.ads, s-gloloc.ads, a-csquin.ads, a-cofuse.adb,
	g-wispch.ads, s-pack03.ads, s-casi16.ads, s-io.adb, a-colien.ads,
	g-alveop.adb, gnat.ads, s-diinio.adb, a-cfdlli.adb, g-pehage.adb,
	a-wtmoio.ads, a-stwiha.adb, a-locale.adb, a-tirsfi.adb, a-nscoty.ads,
	a-rbtgbo.ads, s-pack12.adb, a-cofuse.ads, a-sfteio.ads, s-io.ads,
	g-alveop.ads, a-cfdlli.ads, s-diinio.ads, a-stwiha.ads, g-pehage.ads,
	a-locale.ads, a-tirsfi.ads, s-pack12.ads, s-valuti.adb, g-cppexc.adb,
	system-vxworks-ppc.ads, g-memdum.adb, a-lfteio.ads, s-pack21.adb,
	s-unstyp.ads, s-valuti.ads, g-cppexc.ads, system-hpux-ia64.ads,
	g-memdum.ads, g-soccon.ads, g-altive.ads, a-crbtgo.adb, s-pack21.ads,
	a-llizti.ads, a-numaux-libc-x86.ads, s-expint.adb, s-conca8.adb,
	a-crbtgo.ads, s-pack30.adb, s-vallli.adb, s-geveop.adb, s-expint.ads,
	a-direio.adb, s-conca8.ads, a-widcha.ads, s-pack30.ads, s-vallli.ads,
	s-strhas.adb, s-geveop.ads, g-md5.adb, a-direio.ads, a-numaux.ads,
	s-ransee.adb, a-szbzha.adb, i-cobol.adb, g-busorg.adb, s-strhas.ads,
	g-md5.ads, s-widenu.adb, s-ransee.ads, s-widllu.adb, a-szbzha.ads,
	a-ststio.adb, i-cobol.ads, g-busorg.ads, g-regpat.adb, s-widenu.ads,
	a-secain.adb, s-widllu.ads, s-pack13.adb, g-encstr.adb, a-ztcoau.adb,
	a-ststio.ads, s-widwch.adb, g-regpat.ads, s-atacco.adb, a-cborse.adb,
	a-secain.ads, s-pack13.ads, g-encstr.ads, a-ztcoau.ads, s-widwch.ads,
	g-io_aux.adb, s-atacco.ads, a-ncelfu.ads, interfac.ads, a-cborse.ads,
	g-regexp.adb, s-pack22.adb, a-szuzti.adb, g-io_aux.ads, s-caun32.adb,
	a-nselfu.ads, g-regexp.ads, s-pack22.ads, a-ticoio.adb, a-szuzti.ads,
	g-diopit.adb, s-caun32.ads, s-conca9.adb, a-tags.adb, a-swmwco.ads,
	a-sbecin.adb, s-pack31.adb, s-expuns.adb, a-ticoio.ads, s-valint.adb,
	s-conca9.ads, g-diopit.ads, a-tags.ads, a-nllcef.ads, a-izteio.ads,
	a-sbecin.ads, s-expuns.ads, s-pack31.ads, g-dyntab.adb, s-powtab.ads,
	s-flocon-none.adb, s-valint.ads, a-ssiwti.ads, s-mmosin-mingw.adb,
	s-pack40.adb, s-pack05.adb, a-ztflau.adb, g-dyntab.ads,
	a-szuzti-shared.adb, g-alvevi.ads, a-stwise.adb, s-mmosin-mingw.ads,
	s-pack40.ads, a-diocst.adb, a-ztflau.ads, s-pack05.ads, a-nlcoty.ads,
	a-contai.ads, a-stwisu.adb, g-byorma.adb, a-siwtio.ads, a-stwise.ads,
	s-regpat.adb, g-mbdira.adb, s-pack14.adb, a-diocst.ads, g-flocon.ads,
	g-mbflra.adb, a-ztinau.adb, s-dim.ads, s-mantis.adb, a-stwisu.ads,
	g-byorma.ads, s-atopri.adb, g-wistsp.ads, a-uncdea.ads, s-widcha.adb,
	a-caldel.adb, s-regpat.ads, g-mbdira.ads, a-tiflio.adb, s-pack14.ads,
	s-parame.adb, a-liwtio.ads, s-memory.adb, g-mbflra.ads, a-ztinau.ads,
	a-wtgeau.adb, s-direio.adb, s-mantis.ads, s-atopri.ads, s-widcha.ads,
	a-caldel.ads, s-pack23.adb, a-unccon.ads, a-tiflio.ads, s-parame.ads,
	a-llftio.ads, s-memory.ads, s-regexp.adb, a-wtgeau.ads, a-exexda.adb,
	s-direio.ads, s-pack23.ads, g-stheme.adb, a-tiinio.adb, g-sestin.ads,
	s-regexp.ads, a-wtfiio.adb, a-comutr.adb, a-exexpr.adb, a-tiinio.ads,
	a-ztmoau.adb, a-cohata.ads, a-wtfiio.ads, s-imgrea.adb, ada.ads,
	a-szunau-shared.adb, a-comutr.ads, s-valuns.adb, a-ztmoau.ads,
	system-linux-arm.ads, s-osprim-x32.adb, s-pack41.adb, s-pack06.adb,
	s-imgrea.ads, s-valuns.ads, s-finroo.adb, s-caun16.adb, s-pooglo.adb,
	a-zrstfi.adb, a-suenst.adb, s-pack41.ads, g-binenv.adb, s-pack06.ads,
	a-calari.adb, a-nlcoar.ads, s-finroo.ads, a-timoio.adb, s-caun16.ads,
	s-flocon.adb, a-suenst.ads, a-zrstfi.ads, s-pooglo.ads, s-wchcon.adb,
	s-traceb-hpux.adb, s-pack50.adb, i-fortra.adb, s-pack15.adb,
	a-ngcefu.adb, g-sptavs.ads, g-binenv.ads, s-wchjis.adb, a-calari.ads,
	a-timoio.ads, a-decima.adb, s-flocon.ads, s-wchcon.ads, a-llfzti.ads,
	i-fortra.ads, s-pack50.ads, s-pack15.ads, a-ngcefu.ads, a-cfhase.adb,
	s-wchjis.ads, g-soliop.ads, a-decima.ads, a-chlat9.ads, s-pack24.adb,
	a-nlelfu.ads, a-cfhase.ads, g-locfil.adb, s-atocou-builtin.adb,
	s-memcop.ads, a-szunau.adb, s-pack24.ads, s-imgllb.adb, s-auxdec.adb,
	g-locfil.ads, s-pack33.adb, a-szunau.ads, s-parame-vxworks.adb,
	s-imgllb.ads, a-ciorma.adb, s-auxdec.ads, a-cobove.adb, s-dsaser.ads,
	a-elchha.adb, s-pack33.ads, a-cofuve.adb, s-parame-vxworks.ads,
	a-ciorma.ads, system-darwin-x86.ads, s-multip.adb, a-stwiun-shared.adb,
	a-wichun.adb, a-cobove.ads, s-imgbiu.adb, s-tsmona-mingw.adb,
	a-coormu.adb, a-siocst.adb, s-win32.ads, a-elchha.ads, s-pack42.adb,
	s-pack07.adb, a-cofuve.ads, system-hpux.ads, a-teioed.adb,
	a-convec.adb, g-speche.adb, s-multip.ads, a-stwiun-shared.ads,
	a-wichun.ads, s-imgbiu.ads, a-numeri.ads, a-siocst.ads, a-coormu.ads,
	a-lliwti.ads, s-pack42.ads, s-pack07.ads, a-teioed.ads, a-convec.ads,
	g-speche.ads, g-socthi.adb, a-nucoty.ads, a-szmzco.ads, s-pack51.adb,
	s-osprim-mingw.adb, s-casi64.adb, g-strspl.ads, g-socthi.ads,
	g-socket-dummy.adb, s-pack51.ads, s-dimmks.ads, s-casi64.ads,
	a-wtenau.adb, s-stchop-vxworks.adb, s-pack60.adb,
	system-solaris-sparc.ads, s-pack25.adb, g-socket-dummy.ads,
	a-exstat.adb, a-cofuma.adb, s-tsmona-linux.adb, a-wtenau.ads,
	s-pack60.ads, s-pack25.ads, i-cstrea.adb, a-cofuma.ads, g-exptty.adb,
	a-chzla1.ads, s-pack34.adb, i-cstrea.ads, s-excdeb.adb, a-iteint.ads,
	g-exptty.ads, i-pacdec.adb, s-pack34.ads, s-rident.ads, s-sopco3.adb,
	i-vxwork.ads, s-excdeb.ads, system-linux-ppc.ads, a-swuwti-shared.adb,
	s-widlli.adb, s-pack43.adb, i-pacdec.ads, a-cwila9.ads, s-sopco3.ads,
	a-fwteio.ads, s-widlli.ads, s-pack43.ads, a-suhcin.adb, a-wtdeau.adb,
	g-allein.ads, a-suezst.adb, a-dirval-mingw.adb, g-zspche.adb,
	s-bignum.adb, a-ztedit.adb, g-regist.adb, a-nllefu.ads, a-ztcoio.adb,
	s-pack52.adb, a-llctio.ads, a-nucoar.ads, s-pack17.adb, a-suhcin.ads,
	a-wtdeau.ads, a-suezst.ads, a-dirval.adb, g-zspche.ads, g-regist.ads,
	a-ztedit.ads, s-bignum.ads, a-wtcstr.adb, system.ads, s-pack52.ads,
	a-ztcoio.ads, s-pack17.ads, s-imgboo.adb, a-rbtgso.adb, a-dirval.ads,
	a-cohase.adb, s-pack61.adb, a-wtcstr.ads, s-pack26.adb, s-osprim.ads,
	a-tigeau.adb, s-imgboo.ads, a-nuelfu.ads, a-swfwha.ads, s-commun.adb,
	g-socthi-vxworks.adb, a-rbtgso.ads, a-cohase.ads, g-zstspl.ads,
	s-pack61.ads, s-pack26.ads, a-intnam-dragonfly.ads, s-imglld.adb,
	a-tigeau.ads, s-commun.ads, g-socthi-vxworks.ads, a-cborma.adb,
	a-stwifi.adb, g-moreex.adb, s-pack35.adb, s-imglld.ads, s-valdec.adb,
	a-tifiio.adb, a-cborma.ads, g-moreex.ads, a-stwifi.ads, s-pack35.ads,
	s-sopco4.adb, g-sha224.ads, g-socket.adb, a-intnam-rtems.ads,
	s-finmas.adb, s-valdec.ads, s-addima.adb, a-finali.adb, a-tifiio.ads,
	s-rpc.adb, a-ztflio.adb, s-pack44.adb, s-pack09.adb, a-sblcin.adb,
	s-sopco4.ads, a-textio.adb, g-socket.ads, g-sptabo.ads, s-finmas.ads,
	g-shsh32.adb, s-addima.ads, a-finali.ads, s-mmauni-long.ads, s-rpc.ads,
	a-ztflio.ads, system-djgpp.ads, s-stache.adb, s-pack44.ads,
	s-pack09.ads, a-sblcin.ads, a-textio.ads, a-cidlli.adb, g-shsh32.ads,
	a-chtgbk.adb, a-tiocst.adb, s-pack53.adb, s-pack18.adb, s-stache.ads,
	a-zchhan.adb, s-fatlfl.ads, a-ztinio.adb, s-strops.adb, a-siteio.ads,
	a-cidlli.ads, a-chtgbk.ads, g-ssvety.ads, a-tiocst.ads, s-pack53.ads,
	s-parame-hpux.ads, s-pack18.ads, a-zchhan.ads, s-strops.ads,
	a-ztinio.ads, a-wichha.adb, a-stwima.adb, a-nlrear.ads, a-liteio.ads,
	s-pack62.adb, s-pack27.adb, s-fore.adb, s-vercon.adb, a-wichha.ads,
	a-stwima.ads, s-pack62.ads, system-linux-sparc.ads, s-pack27.ads,
	g-dynhta.adb, s-fore.ads, s-vercon.ads, a-cofuba.adb, a-cimutr.adb,
	i-cpoint.adb, s-imgenu.adb, a-stwibo.adb, s-pack36.adb, i-cstrin.adb,
	s-imgllu.adb, a-suteio-shared.adb, g-excact.adb, s-stoele.adb,
	s-addope.adb, g-dynhta.ads, a-cofuba.ads, a-ztmoio.adb, a-llfwti.ads,
	a-cimutr.ads, i-cpoint.ads, s-imgenu.ads, a-stwibo.ads, a-wttest.adb,
	s-pack36.ads, a-tgdico.ads, s-sopco5.adb, s-scaval.adb, i-cstrin.ads,
	s-imgllu.ads, g-excact.ads, s-stoele.ads, g-deutst.ads, s-addope.ads,
	s-imgwch.adb, g-sha384.ads, a-ztmoio.ads, s-pack45.adb, a-wttest.ads,
	s-sopco5.ads, s-excmac-gcc.adb, s-scaval.ads, a-storio.adb,
	a-coinho.adb, a-btgbso.adb, s-imgwch.ads, s-carun8.adb, memtrack.adb,
	s-pack45.ads, a-sfhcin.ads, s-excmac-gcc.ads, a-storio.ads,
	a-coinho.ads, a-btgbso.ads, s-stratt.adb, s-carun8.ads, a-shcain.adb,
	s-pack54.adb, s-pack19.adb, a-colire.adb, a-tigeli.adb, s-caun64.adb,
	s-stratt.ads, s-fatgen.adb, a-shcain.ads, a-stzunb-shared.adb,
	s-pack54.ads, s-pack19.ads, a-colire.ads, a-calcon.adb, s-caun64.ads,
	s-fatgen.ads, s-pack63.adb, g-arrspl.adb, a-stzunb-shared.ads,
	s-pack28.adb, a-nllrar.ads, a-zzboio.adb, a-zzunio.ads, a-stunau.adb,
	a-calcon.ads, g-cgideb.adb, s-objrea.adb, s-mastop.adb, a-tienau.adb,
	g-altcon.adb, g-arrspl.ads, s-pack63.ads, s-restri.adb, s-pack28.ads,
	a-zzboio.ads, a-stunau.ads, g-cgideb.ads, g-htable.adb, g-sothco.adb,
	s-objrea.ads, g-soliop-solaris.ads, s-mastop.ads, a-tienau.ads,
	system-linux-m68k.ads, g-altcon.ads, s-dmotpr.ads, s-memory-mingw.adb,
	g-cgicoo.adb, s-pack37.adb, s-restri.ads, s-fatllf.ads, s-expmod.adb,
	a-swuwha.adb, a-exextr.adb, a-cfhama.adb, s-gloloc-mingw.adb,
	a-tiboio.adb, g-forstr.adb, g-sothco.ads, a-stzbou.adb, a-nllcty.ads,
	a-suecin.adb, g-htable.ads, s-exctab.adb, a-tiunio.ads, g-cgicoo.ads,
	s-osprim-posix.adb, s-pack37.ads, a-ciormu.adb, s-atocou-x86.adb,
	a-swuwha.ads, s-expmod.ads, a-cfhama.ads, s-ficobl.ads, a-ngcoty.adb,
	g-forstr.ads, a-tiboio.ads, a-calfor.adb, a-stzbou.ads, a-suecin.ads,
	a-conhel.adb, a-crbltr.ads, s-exctab.ads, a-dhfina.ads, s-imgcha.adb,
	s-pack46.adb, a-ciormu.ads, system-linux-sh4.ads, a-chahan.adb,
	a-ngcoty.ads, a-stzunb.adb, a-szfzha.ads, a-calfor.ads, a-cbdlli.adb,
	a-conhel.ads, s-imgcha.ads, s-pack46.ads, a-assert.adb, a-chahan.ads,
	a-stzunb.ads, a-crdlli.adb, s-pack55.adb, a-cbdlli.ads, a-tideau.adb,
	a-assert.ads, ioexcept.ads, s-boarop.ads, g-hesora.adb, a-crdlli.ads,
	s-pack55.ads, a-tideau.ads, g-bubsor.adb, a-wtenio.adb, a-cbsyqu.adb,
	g-hesora.ads, s-pack29.adb, a-nurear.ads, g-catiio.adb, s-stposu.adb,
	g-bubsor.ads, a-wtenio.ads, a-cbsyqu.ads, a-suewst.adb,
	system-vxworks-x86.ads, s-pack29.ads, a-cbmutr.adb, a-cbprqu.adb,
	s-imenne.adb, g-sothco-dummy.adb, g-casuti.adb, g-catiio.ads,
	s-stposu.ads, a-stzsea.adb, s-pack38.adb, a-suewst.ads, s-imgllw.adb,
	a-cbprqu.ads, a-cbmutr.ads, s-imenne.ads, g-sothco-dummy.ads,
	g-casuti.ads, s-htable.adb, s-fatsfl.ads, g-trasym.adb, unchconv.ads,
	a-stzsea.ads, s-arit64.adb, s-pack38.ads, a-nllcar.ads, s-valrea.adb,
	s-imgllw.ads, s-htable.ads, a-sequio.adb, g-trasym.ads, a-ngcoar.adb,
	s-exnllf.adb, s-pack47.adb, s-arit64.ads, g-sercom-mingw.adb,
	s-valrea.ads, g-socthi-mingw.adb, g-bytswa.adb, g-sehash.adb,
	unchdeal.ads, a-sequio.ads, a-ngcoar.ads, s-exnllf.ads, a-wtdeio.adb,
	s-pack47.ads, g-socthi-mingw.ads, a-excpol-abort.adb, a-ztgeau.adb,
	g-bytswa.ads, g-sehash.ads, s-pack56.adb, a-wtdeio.ads, a-ngelfu.adb,
	a-ztgeau.ads, a-cforse.adb, s-filatt.ads, a-stzfix.adb, a-cihase.adb,
	s-pack56.ads, a-sfztio.ads, a-ngelfu.ads, s-trasym-dwarf.adb,
	a-cforse.ads, a-ztfiio.adb, g-timsta.adb, a-stzfix.ads, a-cihase.ads,
	a-ztfiio.ads, system-darwin-arm.ads: Move non-tasking runtime sources
	to libgnat subdirectory.



git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@251902 138bc75d-0d04-0410-961f-82ee72b054a4

1 files changed, 0 insertions, 934 deletions

diff --git a/gcc/ada/s-gearop.adb b/gcc/ada/s-gearop.adb
deleted file mode 100644
index b6d6f22d51b..00000000000
--- a/gcc/ada/s-gearop.adb
+++ /dev/null
@@ -1,934 +0,0 @@
-------------------------------------------------------------------------------
---                                                                          --
---                         GNAT RUN-TIME COMPONENTS                         --
---                                                                          --
---       S Y S T E M . G E N E R I C _ A R R A Y _ O P E R A T I O N S      --
---                                                                          --
---                                 B o d y                                  --
---                                                                          --
---         Copyright (C) 2006-2016, Free Software Foundation, Inc.          --
---                                                                          --
--- GNAT is free software;  you can  redistribute it  and/or modify it under --
--- terms of the  GNU General Public License as published  by the Free Soft- --
--- ware  Foundation;  either version 3,  or (at your option) any later ver- --
--- sion.  GNAT is distributed in the hope that it will be useful, but WITH- --
--- OUT ANY WARRANTY;  without even the  implied warranty of MERCHANTABILITY --
--- or FITNESS FOR A PARTICULAR PURPOSE.                                     --
---                                                                          --
--- As a special exception under Section 7 of GPL version 3, you are granted --
--- additional permissions described in the GCC Runtime Library Exception,   --
--- version 3.1, as published by the Free Software Foundation.               --
---                                                                          --
--- You should have received a copy of the GNU General Public License and    --
--- a copy of the GCC Runtime Library Exception along with this program;     --
--- see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see    --
--- <http://www.gnu.org/licenses/>.                                          --
---                                                                          --
--- GNAT was originally developed  by the GNAT team at  New York University. --
--- Extensive contributions were provided by Ada Core Technologies Inc.      --
---                                                                          --
-------------------------------------------------------------------------------
-
-with Ada.Numerics; use Ada.Numerics;
-package body System.Generic_Array_Operations is
-   function Check_Unit_Last
-     (Index : Integer;
-      Order : Positive;
-      First : Integer) return Integer;
-   pragma Inline_Always (Check_Unit_Last);
-   --  Compute index of last element returned by Unit_Vector or Unit_Matrix.
-   --  A separate function is needed to allow raising Constraint_Error before
-   --  declaring the function result variable. The result variable needs to be
-   --  declared first, to allow front-end inlining.
-
-   --------------
-   -- Diagonal --
-   --------------
-
-   function Diagonal (A : Matrix) return Vector is
-      N : constant Natural := Natural'Min (A'Length (1), A'Length (2));
-   begin
-      return R : Vector (A'First (1) .. A'First (1) + N - 1) do
-         for J in 0 .. N - 1 loop
-            R (R'First + J) := A (A'First (1) + J, A'First (2) + J);
-         end loop;
-      end return;
-   end Diagonal;
-
-   --------------------------
-   -- Square_Matrix_Length --
-   --------------------------
-
-   function Square_Matrix_Length (A : Matrix) return Natural is
-   begin
-      if A'Length (1) /= A'Length (2) then
-         raise Constraint_Error with "matrix is not square";
-      else
-         return A'Length (1);
-      end if;
-   end Square_Matrix_Length;
-
-   ---------------------
-   -- Check_Unit_Last --
-   ---------------------
-
-   function Check_Unit_Last
-      (Index : Integer;
-       Order : Positive;
-       First : Integer) return Integer
-   is
-   begin
-      --  Order the tests carefully to avoid overflow
-
-      if Index < First
-        or else First > Integer'Last - Order + 1
-        or else Index > First + (Order - 1)
-      then
-         raise Constraint_Error;
-      end if;
-
-      return First + (Order - 1);
-   end Check_Unit_Last;
-
-   ---------------------
-   -- Back_Substitute --
-   ---------------------
-
-   procedure Back_Substitute (M, N : in out Matrix) is
-      pragma Assert (M'First (1) = N'First (1)
-                       and then
-                     M'Last  (1) = N'Last (1));
-
-      procedure Sub_Row
-        (M      : in out Matrix;
-         Target : Integer;
-         Source : Integer;
-         Factor : Scalar);
-      --  Elementary row operation that subtracts Factor * M (Source, <>) from
-      --  M (Target, <>)
-
-      -------------
-      -- Sub_Row --
-      -------------
-
-      procedure Sub_Row
-        (M      : in out Matrix;
-         Target : Integer;
-         Source : Integer;
-         Factor : Scalar)
-      is
-      begin
-         for J in M'Range (2) loop
-            M (Target, J) := M (Target, J) - Factor * M (Source, J);
-         end loop;
-      end Sub_Row;
-
-      --  Local declarations
-
-      Max_Col : Integer := M'Last (2);
-
-   --  Start of processing for Back_Substitute
-
-   begin
-      Do_Rows : for Row in reverse M'Range (1) loop
-         Find_Non_Zero : for Col in reverse M'First (2) .. Max_Col loop
-            if Is_Non_Zero (M (Row, Col)) then
-
-               --  Found first non-zero element, so subtract a multiple of this
-               --  element  from all higher rows, to reduce all other elements
-               --  in this column to zero.
-
-               declare
-                  --  We can't use a for loop, as we'd need to iterate to
-                  --  Row - 1, but that expression will overflow if M'First
-                  --  equals Integer'First, which is true for aggregates
-                  --  without explicit bounds..
-
-                  J : Integer := M'First (1);
-
-               begin
-                  while J < Row loop
-                     Sub_Row (N, J, Row, (M (J, Col) / M (Row, Col)));
-                     Sub_Row (M, J, Row, (M (J, Col) / M (Row, Col)));
-                     J := J + 1;
-                  end loop;
-               end;
-
-               --  Avoid potential overflow in the subtraction below
-
-               exit Do_Rows when Col = M'First (2);
-
-               Max_Col := Col - 1;
-
-               exit Find_Non_Zero;
-            end if;
-         end loop Find_Non_Zero;
-      end loop Do_Rows;
-   end Back_Substitute;
-
-   -----------------------
-   -- Forward_Eliminate --
-   -----------------------
-
-   procedure Forward_Eliminate
-     (M   : in out Matrix;
-      N   : in out Matrix;
-      Det : out Scalar)
-   is
-      pragma Assert (M'First (1) = N'First (1)
-                       and then
-                     M'Last  (1) = N'Last (1));
-
-      --  The following are variations of the elementary matrix row operations:
-      --  row switching, row multiplication and row addition. Because in this
-      --  algorithm the addition factor is always a negated value, we chose to
-      --  use  row subtraction instead. Similarly, instead of multiplying by
-      --  a reciprocal, we divide.
-
-      procedure Sub_Row
-        (M      : in out Matrix;
-         Target : Integer;
-         Source : Integer;
-         Factor : Scalar);
-      --  Subtrace Factor * M (Source, <>) from M (Target, <>)
-
-      procedure Divide_Row
-        (M, N  : in out Matrix;
-         Row   : Integer;
-         Scale : Scalar);
-      --  Divide M (Row) and N (Row) by Scale, and update Det
-
-      procedure Switch_Row
-        (M, N  : in out Matrix;
-         Row_1 : Integer;
-         Row_2 : Integer);
-      --  Exchange M (Row_1) and N (Row_1) with M (Row_2) and N (Row_2),
-      --  negating Det in the process.
-
-      -------------
-      -- Sub_Row --
-      -------------
-
-      procedure Sub_Row
-        (M      : in out Matrix;
-         Target : Integer;
-         Source : Integer;
-         Factor : Scalar)
-      is
-      begin
-         for J in M'Range (2) loop
-            M (Target, J) := M (Target, J) - Factor * M (Source, J);
-         end loop;
-      end Sub_Row;
-
-      ----------------
-      -- Divide_Row --
-      ----------------
-
-      procedure Divide_Row
-        (M, N  : in out Matrix;
-         Row   : Integer;
-         Scale : Scalar)
-      is
-      begin
-         Det := Det * Scale;
-
-         for J in M'Range (2) loop
-            M (Row, J) := M (Row, J) / Scale;
-         end loop;
-
-         for J in N'Range (2) loop
-            N (Row - M'First (1) + N'First (1), J) :=
-              N (Row - M'First (1) + N'First (1), J) / Scale;
-         end loop;
-      end Divide_Row;
-
-      ----------------
-      -- Switch_Row --
-      ----------------
-
-      procedure Switch_Row
-        (M, N  : in out Matrix;
-         Row_1 : Integer;
-         Row_2 : Integer)
-      is
-         procedure Swap (X, Y : in out Scalar);
-         --  Exchange the values of X and Y
-
-         ----------
-         -- Swap --
-         ----------
-
-         procedure Swap (X, Y : in out Scalar) is
-            T : constant Scalar := X;
-         begin
-            X := Y;
-            Y := T;
-         end Swap;
-
-      --  Start of processing for Switch_Row
-
-      begin
-         if Row_1 /= Row_2 then
-            Det := Zero - Det;
-
-            for J in M'Range (2) loop
-               Swap (M (Row_1, J), M (Row_2, J));
-            end loop;
-
-            for J in N'Range (2) loop
-               Swap (N (Row_1 - M'First (1) + N'First (1), J),
-                     N (Row_2 - M'First (1) + N'First (1), J));
-            end loop;
-         end if;
-      end Switch_Row;
-
-      --  Local declarations
-
-      Row : Integer := M'First (1);
-
-   --  Start of processing for Forward_Eliminate
-
-   begin
-      Det := One;
-
-      for J in M'Range (2) loop
-         declare
-            Max_Row : Integer := Row;
-            Max_Abs : Real'Base := 0.0;
-
-         begin
-            --  Find best pivot in column J, starting in row Row
-
-            for K in Row .. M'Last (1) loop
-               declare
-                  New_Abs : constant Real'Base := abs M (K, J);
-               begin
-                  if Max_Abs < New_Abs then
-                     Max_Abs := New_Abs;
-                     Max_Row := K;
-                  end if;
-               end;
-            end loop;
-
-            if Max_Abs > 0.0 then
-               Switch_Row (M, N, Row, Max_Row);
-
-               --  The temporaries below are necessary to force a copy of the
-               --  value and avoid improper aliasing.
-
-               declare
-                  Scale : constant Scalar := M (Row, J);
-               begin
-                  Divide_Row (M, N, Row, Scale);
-               end;
-
-               for U in Row + 1 .. M'Last (1) loop
-                  declare
-                     Factor : constant Scalar := M (U, J);
-                  begin
-                     Sub_Row (N, U, Row, Factor);
-                     Sub_Row (M, U, Row, Factor);
-                  end;
-               end loop;
-
-               exit when Row >= M'Last (1);
-
-               Row := Row + 1;
-
-            else
-               --  Set zero (note that we do not have literals)
-
-               Det := Zero;
-            end if;
-         end;
-      end loop;
-   end Forward_Eliminate;
-
-   -------------------
-   -- Inner_Product --
-   -------------------
-
-   function Inner_Product
-     (Left  : Left_Vector;
-      Right : Right_Vector) return  Result_Scalar
-   is
-      R : Result_Scalar := Zero;
-
-   begin
-      if Left'Length /= Right'Length then
-         raise Constraint_Error with
-            "vectors are of different length in inner product";
-      end if;
-
-      for J in Left'Range loop
-         R := R + Left (J) * Right (J - Left'First + Right'First);
-      end loop;
-
-      return R;
-   end Inner_Product;
-
-   -------------
-   -- L2_Norm --
-   -------------
-
-   function L2_Norm (X : X_Vector) return Result_Real'Base is
-      Sum : Result_Real'Base := 0.0;
-
-   begin
-      for J in X'Range loop
-         Sum := Sum + Result_Real'Base (abs X (J))**2;
-      end loop;
-
-      return Sqrt (Sum);
-   end L2_Norm;
-
-   ----------------------------------
-   -- Matrix_Elementwise_Operation --
-   ----------------------------------
-
-   function Matrix_Elementwise_Operation (X : X_Matrix) return Result_Matrix is
-   begin
-      return R : Result_Matrix (X'Range (1), X'Range (2)) do
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) := Operation (X (J, K));
-            end loop;
-         end loop;
-      end return;
-   end Matrix_Elementwise_Operation;
-
-   ----------------------------------
-   -- Vector_Elementwise_Operation --
-   ----------------------------------
-
-   function Vector_Elementwise_Operation (X : X_Vector) return Result_Vector is
-   begin
-      return R : Result_Vector (X'Range) do
-         for J in R'Range loop
-            R (J) := Operation (X (J));
-         end loop;
-      end return;
-   end Vector_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Matrix_Matrix_Elementwise_Operation --
-   -----------------------------------------
-
-   function Matrix_Matrix_Elementwise_Operation
-     (Left  : Left_Matrix;
-      Right : Right_Matrix) return Result_Matrix
-   is
-   begin
-      return R : Result_Matrix (Left'Range (1), Left'Range (2)) do
-         if Left'Length (1) /= Right'Length (1)
-              or else
-            Left'Length (2) /= Right'Length (2)
-         then
-            raise Constraint_Error with
-              "matrices are of different dimension in elementwise operation";
-         end if;
-
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) :=
-                 Operation
-                   (Left (J, K),
-                    Right
-                      (J - R'First (1) + Right'First (1),
-                       K - R'First (2) + Right'First (2)));
-            end loop;
-         end loop;
-      end return;
-   end Matrix_Matrix_Elementwise_Operation;
-
-   ------------------------------------------------
-   -- Matrix_Matrix_Scalar_Elementwise_Operation --
-   ------------------------------------------------
-
-   function Matrix_Matrix_Scalar_Elementwise_Operation
-     (X : X_Matrix;
-      Y : Y_Matrix;
-      Z : Z_Scalar) return Result_Matrix
-   is
-   begin
-      return R : Result_Matrix (X'Range (1), X'Range (2)) do
-         if X'Length (1) /= Y'Length (1)
-              or else
-            X'Length (2) /= Y'Length (2)
-         then
-            raise Constraint_Error with
-              "matrices are of different dimension in elementwise operation";
-         end if;
-
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) :=
-                 Operation
-                   (X (J, K),
-                    Y (J - R'First (1) + Y'First (1),
-                       K - R'First (2) + Y'First (2)),
-                    Z);
-            end loop;
-         end loop;
-      end return;
-   end Matrix_Matrix_Scalar_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Vector_Vector_Elementwise_Operation --
-   -----------------------------------------
-
-   function Vector_Vector_Elementwise_Operation
-     (Left  : Left_Vector;
-      Right : Right_Vector) return Result_Vector
-   is
-   begin
-      return R : Result_Vector (Left'Range) do
-         if Left'Length /= Right'Length then
-            raise Constraint_Error with
-              "vectors are of different length in elementwise operation";
-         end if;
-
-         for J in R'Range loop
-            R (J) := Operation (Left (J), Right (J - R'First + Right'First));
-         end loop;
-      end return;
-   end Vector_Vector_Elementwise_Operation;
-
-   ------------------------------------------------
-   -- Vector_Vector_Scalar_Elementwise_Operation --
-   ------------------------------------------------
-
-   function Vector_Vector_Scalar_Elementwise_Operation
-     (X : X_Vector;
-      Y : Y_Vector;
-      Z : Z_Scalar) return Result_Vector is
-   begin
-      return R : Result_Vector (X'Range) do
-         if X'Length /= Y'Length then
-            raise Constraint_Error with
-              "vectors are of different length in elementwise operation";
-         end if;
-
-         for J in R'Range loop
-            R (J) := Operation (X (J), Y (J - X'First + Y'First), Z);
-         end loop;
-      end return;
-   end Vector_Vector_Scalar_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Matrix_Scalar_Elementwise_Operation --
-   -----------------------------------------
-
-   function Matrix_Scalar_Elementwise_Operation
-     (Left  : Left_Matrix;
-      Right : Right_Scalar) return Result_Matrix
-   is
-   begin
-      return R : Result_Matrix (Left'Range (1), Left'Range (2)) do
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) := Operation (Left (J, K), Right);
-            end loop;
-         end loop;
-      end return;
-   end Matrix_Scalar_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Vector_Scalar_Elementwise_Operation --
-   -----------------------------------------
-
-   function Vector_Scalar_Elementwise_Operation
-     (Left  : Left_Vector;
-      Right : Right_Scalar) return Result_Vector
-   is
-   begin
-      return R : Result_Vector (Left'Range) do
-         for J in R'Range loop
-            R (J) := Operation (Left (J), Right);
-         end loop;
-      end return;
-   end Vector_Scalar_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Scalar_Matrix_Elementwise_Operation --
-   -----------------------------------------
-
-   function Scalar_Matrix_Elementwise_Operation
-     (Left  : Left_Scalar;
-      Right : Right_Matrix) return Result_Matrix
-   is
-   begin
-      return R : Result_Matrix (Right'Range (1), Right'Range (2)) do
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) := Operation (Left, Right (J, K));
-            end loop;
-         end loop;
-      end return;
-   end Scalar_Matrix_Elementwise_Operation;
-
-   -----------------------------------------
-   -- Scalar_Vector_Elementwise_Operation --
-   -----------------------------------------
-
-   function Scalar_Vector_Elementwise_Operation
-     (Left  : Left_Scalar;
-      Right : Right_Vector) return Result_Vector
-   is
-   begin
-      return R : Result_Vector (Right'Range) do
-         for J in R'Range loop
-            R (J) := Operation (Left, Right (J));
-         end loop;
-      end return;
-   end Scalar_Vector_Elementwise_Operation;
-
-   ----------
-   -- Sqrt --
-   ----------
-
-   function Sqrt (X : Real'Base) return Real'Base is
-      Root, Next : Real'Base;
-
-   begin
-      --  Be defensive: any comparisons with NaN values will yield False.
-
-      if not (X > 0.0) then
-         if X = 0.0 then
-            return X;
-         else
-            raise Argument_Error;
-         end if;
-
-      elsif X > Real'Base'Last then
-
-         --  X is infinity, which is its own square root
-
-         return X;
-      end if;
-
-      --  Compute an initial estimate based on:
-
-      --     X = M * R**E and Sqrt (X) = Sqrt (M) * R**(E / 2.0),
-
-      --  where M is the mantissa, R is the radix and E the exponent.
-
-      --  By ignoring the mantissa and ignoring the case of an odd
-      --  exponent, we get a final error that is at most R. In other words,
-      --  the result has about a single bit precision.
-
-      Root := Real'Base (Real'Machine_Radix) ** (Real'Exponent (X) / 2);
-
-      --  Because of the poor initial estimate, use the Babylonian method of
-      --  computing the square root, as it is stable for all inputs. Every step
-      --  will roughly double the precision of the result. Just a few steps
-      --  suffice in most cases. Eight iterations should give about 2**8 bits
-      --  of precision.
-
-      for J in 1 .. 8 loop
-         Next := (Root + X / Root) / 2.0;
-         exit when Root = Next;
-         Root := Next;
-      end loop;
-
-      return Root;
-   end Sqrt;
-
-   ---------------------------
-   -- Matrix_Matrix_Product --
-   ---------------------------
-
-   function Matrix_Matrix_Product
-     (Left  : Left_Matrix;
-      Right : Right_Matrix) return Result_Matrix
-   is
-   begin
-      return R : Result_Matrix (Left'Range (1), Right'Range (2)) do
-         if Left'Length (2) /= Right'Length (1) then
-            raise Constraint_Error with
-              "incompatible dimensions in matrix multiplication";
-         end if;
-
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               declare
-                  S : Result_Scalar := Zero;
-
-               begin
-                  for M in Left'Range (2) loop
-                     S := S + Left (J, M) *
-                                Right
-                                  (M - Left'First (2) + Right'First (1), K);
-                  end loop;
-
-                  R (J, K) := S;
-               end;
-            end loop;
-         end loop;
-      end return;
-   end  Matrix_Matrix_Product;
-
-   ----------------------------
-   -- Matrix_Vector_Solution --
-   ----------------------------
-
-   function Matrix_Vector_Solution (A : Matrix; X : Vector) return Vector is
-      N   : constant Natural := A'Length (1);
-      MA  : Matrix := A;
-      MX  : Matrix (A'Range (1), 1 .. 1);
-      R   : Vector (A'Range (2));
-      Det : Scalar;
-
-   begin
-      if A'Length (2) /= N then
-         raise Constraint_Error with "matrix is not square";
-      end if;
-
-      if X'Length /= N then
-         raise Constraint_Error with "incompatible vector length";
-      end if;
-
-      for J in 0 .. MX'Length (1) - 1 loop
-         MX (MX'First (1) + J, 1) := X (X'First + J);
-      end loop;
-
-      Forward_Eliminate (MA, MX, Det);
-
-      if Det = Zero then
-         raise Constraint_Error with "matrix is singular";
-      end if;
-
-      Back_Substitute (MA, MX);
-
-      for J in 0 .. R'Length - 1 loop
-         R (R'First + J) := MX (MX'First (1) + J, 1);
-      end loop;
-
-      return R;
-   end Matrix_Vector_Solution;
-
-   ----------------------------
-   -- Matrix_Matrix_Solution --
-   ----------------------------
-
-   function Matrix_Matrix_Solution (A, X : Matrix) return Matrix is
-      N   : constant Natural := A'Length (1);
-      MA  : Matrix (A'Range (2), A'Range (2));
-      MB  : Matrix (A'Range (2), X'Range (2));
-      Det : Scalar;
-
-   begin
-      if A'Length (2) /= N then
-         raise Constraint_Error with "matrix is not square";
-      end if;
-
-      if X'Length (1) /= N then
-         raise Constraint_Error with "matrices have unequal number of rows";
-      end if;
-
-      for J in 0 .. A'Length (1) - 1 loop
-         for K in MA'Range (2) loop
-            MA (MA'First (1) + J, K) := A (A'First (1) + J, K);
-         end loop;
-
-         for K in MB'Range (2) loop
-            MB (MB'First (1) + J, K) := X (X'First (1) + J, K);
-         end loop;
-      end loop;
-
-      Forward_Eliminate (MA, MB, Det);
-
-      if Det = Zero then
-         raise Constraint_Error with "matrix is singular";
-      end if;
-
-      Back_Substitute (MA, MB);
-
-      return MB;
-   end Matrix_Matrix_Solution;
-
-   ---------------------------
-   -- Matrix_Vector_Product --
-   ---------------------------
-
-   function Matrix_Vector_Product
-     (Left  : Matrix;
-      Right : Right_Vector) return Result_Vector
-   is
-   begin
-      return R : Result_Vector (Left'Range (1)) do
-         if Left'Length (2) /= Right'Length then
-            raise Constraint_Error with
-              "incompatible dimensions in matrix-vector multiplication";
-         end if;
-
-         for J in Left'Range (1) loop
-            declare
-               S : Result_Scalar := Zero;
-
-            begin
-               for K in Left'Range (2) loop
-                  S := S + Left (J, K)
-                         * Right (K - Left'First (2) + Right'First);
-               end loop;
-
-               R (J) := S;
-            end;
-         end loop;
-      end return;
-   end Matrix_Vector_Product;
-
-   -------------------
-   -- Outer_Product --
-   -------------------
-
-   function Outer_Product
-     (Left  : Left_Vector;
-      Right : Right_Vector) return Matrix
-   is
-   begin
-      return R : Matrix (Left'Range, Right'Range) do
-         for J in R'Range (1) loop
-            for K in R'Range (2) loop
-               R (J, K) := Left (J) * Right (K);
-            end loop;
-         end loop;
-      end return;
-   end Outer_Product;
-
-   -----------------
-   -- Swap_Column --
-   -----------------
-
-   procedure Swap_Column (A : in out Matrix; Left, Right : Integer) is
-      Temp : Scalar;
-   begin
-      for J in A'Range (1) loop
-         Temp := A (J, Left);
-         A (J, Left) := A (J, Right);
-         A (J, Right) := Temp;
-      end loop;
-   end Swap_Column;
-
-   ---------------
-   -- Transpose --
-   ---------------
-
-   procedure Transpose (A : Matrix; R : out Matrix) is
-   begin
-      for J in R'Range (1) loop
-         for K in R'Range (2) loop
-            R (J, K) := A (K - R'First (2) + A'First (1),
-                           J - R'First (1) + A'First (2));
-         end loop;
-      end loop;
-   end Transpose;
-
-   -------------------------------
-   -- Update_Matrix_With_Matrix --
-   -------------------------------
-
-   procedure Update_Matrix_With_Matrix (X : in out X_Matrix; Y : Y_Matrix) is
-   begin
-      if X'Length (1) /= Y'Length (1)
-           or else
-         X'Length (2) /= Y'Length (2)
-      then
-         raise Constraint_Error with
-           "matrices are of different dimension in update operation";
-      end if;
-
-      for J in X'Range (1) loop
-         for K in X'Range (2) loop
-            Update (X (J, K), Y (J - X'First (1) + Y'First (1),
-                                 K - X'First (2) + Y'First (2)));
-         end loop;
-      end loop;
-   end Update_Matrix_With_Matrix;
-
-   -------------------------------
-   -- Update_Vector_With_Vector --
-   -------------------------------
-
-   procedure Update_Vector_With_Vector (X : in out X_Vector; Y : Y_Vector) is
-   begin
-      if X'Length /= Y'Length then
-         raise Constraint_Error with
-           "vectors are of different length in update operation";
-      end if;
-
-      for J in X'Range loop
-         Update (X (J), Y (J - X'First + Y'First));
-      end loop;
-   end Update_Vector_With_Vector;
-
-   -----------------
-   -- Unit_Matrix --
-   -----------------
-
-   function Unit_Matrix
-     (Order   : Positive;
-      First_1 : Integer := 1;
-      First_2 : Integer := 1) return Matrix
-   is
-   begin
-      return R : Matrix (First_1 .. Check_Unit_Last (First_1, Order, First_1),
-                         First_2 .. Check_Unit_Last (First_2, Order, First_2))
-      do
-         R := (others => (others => Zero));
-
-         for J in 0 .. Order - 1 loop
-            R (First_1 + J, First_2 + J) := One;
-         end loop;
-      end return;
-   end Unit_Matrix;
-
-   -----------------
-   -- Unit_Vector --
-   -----------------
-
-   function Unit_Vector
-     (Index : Integer;
-      Order : Positive;
-      First : Integer := 1) return Vector
-   is
-   begin
-      return R : Vector (First .. Check_Unit_Last (Index, Order, First)) do
-         R := (others => Zero);
-         R (Index) := One;
-      end return;
-   end Unit_Vector;
-
-   ---------------------------
-   -- Vector_Matrix_Product --
-   ---------------------------
-
-   function Vector_Matrix_Product
-     (Left  : Left_Vector;
-      Right : Matrix) return Result_Vector
-   is
-   begin
-      return R : Result_Vector (Right'Range (2)) do
-         if Left'Length /= Right'Length (1) then
-            raise Constraint_Error with
-              "incompatible dimensions in vector-matrix multiplication";
-         end if;
-
-         for J in Right'Range (2) loop
-            declare
-               S : Result_Scalar := Zero;
-
-            begin
-               for K in Right'Range (1) loop
-                  S := S + Left (K - Right'First (1)
-                                   + Left'First) * Right (K, J);
-               end loop;
-
-               R (J) := S;
-            end;
-         end loop;
-      end return;
-   end Vector_Matrix_Product;
-
-end System.Generic_Array_Operations;