--- /dev/null
+#! /bin/sh
+
+# Guess values for system-dependent variables and create Makefiles.
+# Generated automatically using autoconf version 2.13
+# Copyright (C) 1992, 93, 94, 95, 96 Free Software Foundation, Inc.
+#
+# This configure script is free software; the Free Software Foundation
+# gives unlimited permission to copy, distribute and modify it.
+
+# Defaults:
+ac_help=
+ac_default_prefix=/usr/local
+# Any additions from configure.in:
+ac_help="$ac_help
+ --with-deal=DIR where DIR is the root of the deal.II library"
+
+# Initialize some variables set by options.
+# The variables have the same names as the options, with
+# dashes changed to underlines.
+build=NONE
+cache_file=./config.cache
+exec_prefix=NONE
+host=NONE
+no_create=
+nonopt=NONE
+no_recursion=
+prefix=NONE
+program_prefix=NONE
+program_suffix=NONE
+program_transform_name=s,x,x,
+silent=
+site=
+srcdir=
+target=NONE
+verbose=
+x_includes=NONE
+x_libraries=NONE
+bindir='${exec_prefix}/bin'
+sbindir='${exec_prefix}/sbin'
+libexecdir='${exec_prefix}/libexec'
+datadir='${prefix}/share'
+sysconfdir='${prefix}/etc'
+sharedstatedir='${prefix}/com'
+localstatedir='${prefix}/var'
+libdir='${exec_prefix}/lib'
+includedir='${prefix}/include'
+oldincludedir='/usr/include'
+infodir='${prefix}/info'
+mandir='${prefix}/man'
+
+# Initialize some other variables.
+subdirs=
+MFLAGS= MAKEFLAGS=
+SHELL=${CONFIG_SHELL-/bin/sh}
+# Maximum number of lines to put in a shell here document.
+ac_max_here_lines=12
+
+ac_prev=
+for ac_option
+do
+
+ # If the previous option needs an argument, assign it.
+ if test -n "$ac_prev"; then
+ eval "$ac_prev=\$ac_option"
+ ac_prev=
+ continue
+ fi
+
+ case "$ac_option" in
+ -*=*) ac_optarg=`echo "$ac_option" | sed 's/[-_a-zA-Z0-9]*=//'` ;;
+ *) ac_optarg= ;;
+ esac
+
+ # Accept the important Cygnus configure options, so we can diagnose typos.
+
+ case "$ac_option" in
+
+ -bindir | --bindir | --bindi | --bind | --bin | --bi)
+ ac_prev=bindir ;;
+ -bindir=* | --bindir=* | --bindi=* | --bind=* | --bin=* | --bi=*)
+ bindir="$ac_optarg" ;;
+
+ -build | --build | --buil | --bui | --bu)
+ ac_prev=build ;;
+ -build=* | --build=* | --buil=* | --bui=* | --bu=*)
+ build="$ac_optarg" ;;
+
+ -cache-file | --cache-file | --cache-fil | --cache-fi \
+ | --cache-f | --cache- | --cache | --cach | --cac | --ca | --c)
+ ac_prev=cache_file ;;
+ -cache-file=* | --cache-file=* | --cache-fil=* | --cache-fi=* \
+ | --cache-f=* | --cache-=* | --cache=* | --cach=* | --cac=* | --ca=* | --c=*)
+ cache_file="$ac_optarg" ;;
+
+ -datadir | --datadir | --datadi | --datad | --data | --dat | --da)
+ ac_prev=datadir ;;
+ -datadir=* | --datadir=* | --datadi=* | --datad=* | --data=* | --dat=* \
+ | --da=*)
+ datadir="$ac_optarg" ;;
+
+ -disable-* | --disable-*)
+ ac_feature=`echo $ac_option|sed -e 's/-*disable-//'`
+ # Reject names that are not valid shell variable names.
+ if test -n "`echo $ac_feature| sed 's/[-a-zA-Z0-9_]//g'`"; then
+ { echo "configure: error: $ac_feature: invalid feature name" 1>&2; exit 1; }
+ fi
+ ac_feature=`echo $ac_feature| sed 's/-/_/g'`
+ eval "enable_${ac_feature}=no" ;;
+
+ -enable-* | --enable-*)
+ ac_feature=`echo $ac_option|sed -e 's/-*enable-//' -e 's/=.*//'`
+ # Reject names that are not valid shell variable names.
+ if test -n "`echo $ac_feature| sed 's/[-_a-zA-Z0-9]//g'`"; then
+ { echo "configure: error: $ac_feature: invalid feature name" 1>&2; exit 1; }
+ fi
+ ac_feature=`echo $ac_feature| sed 's/-/_/g'`
+ case "$ac_option" in
+ *=*) ;;
+ *) ac_optarg=yes ;;
+ esac
+ eval "enable_${ac_feature}='$ac_optarg'" ;;
+
+ -exec-prefix | --exec_prefix | --exec-prefix | --exec-prefi \
+ | --exec-pref | --exec-pre | --exec-pr | --exec-p | --exec- \
+ | --exec | --exe | --ex)
+ ac_prev=exec_prefix ;;
+ -exec-prefix=* | --exec_prefix=* | --exec-prefix=* | --exec-prefi=* \
+ | --exec-pref=* | --exec-pre=* | --exec-pr=* | --exec-p=* | --exec-=* \
+ | --exec=* | --exe=* | --ex=*)
+ exec_prefix="$ac_optarg" ;;
+
+ -gas | --gas | --ga | --g)
+ # Obsolete; use --with-gas.
+ with_gas=yes ;;
+
+ -help | --help | --hel | --he)
+ # Omit some internal or obsolete options to make the list less imposing.
+ # This message is too long to be a string in the A/UX 3.1 sh.
+ cat << EOF
+Usage: configure [options] [host]
+Options: [defaults in brackets after descriptions]
+Configuration:
+ --cache-file=FILE cache test results in FILE
+ --help print this message
+ --no-create do not create output files
+ --quiet, --silent do not print \`checking...' messages
+ --version print the version of autoconf that created configure
+Directory and file names:
+ --prefix=PREFIX install architecture-independent files in PREFIX
+ [$ac_default_prefix]
+ --exec-prefix=EPREFIX install architecture-dependent files in EPREFIX
+ [same as prefix]
+ --bindir=DIR user executables in DIR [EPREFIX/bin]
+ --sbindir=DIR system admin executables in DIR [EPREFIX/sbin]
+ --libexecdir=DIR program executables in DIR [EPREFIX/libexec]
+ --datadir=DIR read-only architecture-independent data in DIR
+ [PREFIX/share]
+ --sysconfdir=DIR read-only single-machine data in DIR [PREFIX/etc]
+ --sharedstatedir=DIR modifiable architecture-independent data in DIR
+ [PREFIX/com]
+ --localstatedir=DIR modifiable single-machine data in DIR [PREFIX/var]
+ --libdir=DIR object code libraries in DIR [EPREFIX/lib]
+ --includedir=DIR C header files in DIR [PREFIX/include]
+ --oldincludedir=DIR C header files for non-gcc in DIR [/usr/include]
+ --infodir=DIR info documentation in DIR [PREFIX/info]
+ --mandir=DIR man documentation in DIR [PREFIX/man]
+ --srcdir=DIR find the sources in DIR [configure dir or ..]
+ --program-prefix=PREFIX prepend PREFIX to installed program names
+ --program-suffix=SUFFIX append SUFFIX to installed program names
+ --program-transform-name=PROGRAM
+ run sed PROGRAM on installed program names
+EOF
+ cat << EOF
+Host type:
+ --build=BUILD configure for building on BUILD [BUILD=HOST]
+ --host=HOST configure for HOST [guessed]
+ --target=TARGET configure for TARGET [TARGET=HOST]
+Features and packages:
+ --disable-FEATURE do not include FEATURE (same as --enable-FEATURE=no)
+ --enable-FEATURE[=ARG] include FEATURE [ARG=yes]
+ --with-PACKAGE[=ARG] use PACKAGE [ARG=yes]
+ --without-PACKAGE do not use PACKAGE (same as --with-PACKAGE=no)
+ --x-includes=DIR X include files are in DIR
+ --x-libraries=DIR X library files are in DIR
+EOF
+ if test -n "$ac_help"; then
+ echo "--enable and --with options recognized:$ac_help"
+ fi
+ exit 0 ;;
+
+ -host | --host | --hos | --ho)
+ ac_prev=host ;;
+ -host=* | --host=* | --hos=* | --ho=*)
+ host="$ac_optarg" ;;
+
+ -includedir | --includedir | --includedi | --included | --include \
+ | --includ | --inclu | --incl | --inc)
+ ac_prev=includedir ;;
+ -includedir=* | --includedir=* | --includedi=* | --included=* | --include=* \
+ | --includ=* | --inclu=* | --incl=* | --inc=*)
+ includedir="$ac_optarg" ;;
+
+ -infodir | --infodir | --infodi | --infod | --info | --inf)
+ ac_prev=infodir ;;
+ -infodir=* | --infodir=* | --infodi=* | --infod=* | --info=* | --inf=*)
+ infodir="$ac_optarg" ;;
+
+ -libdir | --libdir | --libdi | --libd)
+ ac_prev=libdir ;;
+ -libdir=* | --libdir=* | --libdi=* | --libd=*)
+ libdir="$ac_optarg" ;;
+
+ -libexecdir | --libexecdir | --libexecdi | --libexecd | --libexec \
+ | --libexe | --libex | --libe)
+ ac_prev=libexecdir ;;
+ -libexecdir=* | --libexecdir=* | --libexecdi=* | --libexecd=* | --libexec=* \
+ | --libexe=* | --libex=* | --libe=*)
+ libexecdir="$ac_optarg" ;;
+
+ -localstatedir | --localstatedir | --localstatedi | --localstated \
+ | --localstate | --localstat | --localsta | --localst \
+ | --locals | --local | --loca | --loc | --lo)
+ ac_prev=localstatedir ;;
+ -localstatedir=* | --localstatedir=* | --localstatedi=* | --localstated=* \
+ | --localstate=* | --localstat=* | --localsta=* | --localst=* \
+ | --locals=* | --local=* | --loca=* | --loc=* | --lo=*)
+ localstatedir="$ac_optarg" ;;
+
+ -mandir | --mandir | --mandi | --mand | --man | --ma | --m)
+ ac_prev=mandir ;;
+ -mandir=* | --mandir=* | --mandi=* | --mand=* | --man=* | --ma=* | --m=*)
+ mandir="$ac_optarg" ;;
+
+ -nfp | --nfp | --nf)
+ # Obsolete; use --without-fp.
+ with_fp=no ;;
+
+ -no-create | --no-create | --no-creat | --no-crea | --no-cre \
+ | --no-cr | --no-c)
+ no_create=yes ;;
+
+ -no-recursion | --no-recursion | --no-recursio | --no-recursi \
+ | --no-recurs | --no-recur | --no-recu | --no-rec | --no-re | --no-r)
+ no_recursion=yes ;;
+
+ -oldincludedir | --oldincludedir | --oldincludedi | --oldincluded \
+ | --oldinclude | --oldinclud | --oldinclu | --oldincl | --oldinc \
+ | --oldin | --oldi | --old | --ol | --o)
+ ac_prev=oldincludedir ;;
+ -oldincludedir=* | --oldincludedir=* | --oldincludedi=* | --oldincluded=* \
+ | --oldinclude=* | --oldinclud=* | --oldinclu=* | --oldincl=* | --oldinc=* \
+ | --oldin=* | --oldi=* | --old=* | --ol=* | --o=*)
+ oldincludedir="$ac_optarg" ;;
+
+ -prefix | --prefix | --prefi | --pref | --pre | --pr | --p)
+ ac_prev=prefix ;;
+ -prefix=* | --prefix=* | --prefi=* | --pref=* | --pre=* | --pr=* | --p=*)
+ prefix="$ac_optarg" ;;
+
+ -program-prefix | --program-prefix | --program-prefi | --program-pref \
+ | --program-pre | --program-pr | --program-p)
+ ac_prev=program_prefix ;;
+ -program-prefix=* | --program-prefix=* | --program-prefi=* \
+ | --program-pref=* | --program-pre=* | --program-pr=* | --program-p=*)
+ program_prefix="$ac_optarg" ;;
+
+ -program-suffix | --program-suffix | --program-suffi | --program-suff \
+ | --program-suf | --program-su | --program-s)
+ ac_prev=program_suffix ;;
+ -program-suffix=* | --program-suffix=* | --program-suffi=* \
+ | --program-suff=* | --program-suf=* | --program-su=* | --program-s=*)
+ program_suffix="$ac_optarg" ;;
+
+ -program-transform-name | --program-transform-name \
+ | --program-transform-nam | --program-transform-na \
+ | --program-transform-n | --program-transform- \
+ | --program-transform | --program-transfor \
+ | --program-transfo | --program-transf \
+ | --program-trans | --program-tran \
+ | --progr-tra | --program-tr | --program-t)
+ ac_prev=program_transform_name ;;
+ -program-transform-name=* | --program-transform-name=* \
+ | --program-transform-nam=* | --program-transform-na=* \
+ | --program-transform-n=* | --program-transform-=* \
+ | --program-transform=* | --program-transfor=* \
+ | --program-transfo=* | --program-transf=* \
+ | --program-trans=* | --program-tran=* \
+ | --progr-tra=* | --program-tr=* | --program-t=*)
+ program_transform_name="$ac_optarg" ;;
+
+ -q | -quiet | --quiet | --quie | --qui | --qu | --q \
+ | -silent | --silent | --silen | --sile | --sil)
+ silent=yes ;;
+
+ -sbindir | --sbindir | --sbindi | --sbind | --sbin | --sbi | --sb)
+ ac_prev=sbindir ;;
+ -sbindir=* | --sbindir=* | --sbindi=* | --sbind=* | --sbin=* \
+ | --sbi=* | --sb=*)
+ sbindir="$ac_optarg" ;;
+
+ -sharedstatedir | --sharedstatedir | --sharedstatedi \
+ | --sharedstated | --sharedstate | --sharedstat | --sharedsta \
+ | --sharedst | --shareds | --shared | --share | --shar \
+ | --sha | --sh)
+ ac_prev=sharedstatedir ;;
+ -sharedstatedir=* | --sharedstatedir=* | --sharedstatedi=* \
+ | --sharedstated=* | --sharedstate=* | --sharedstat=* | --sharedsta=* \
+ | --sharedst=* | --shareds=* | --shared=* | --share=* | --shar=* \
+ | --sha=* | --sh=*)
+ sharedstatedir="$ac_optarg" ;;
+
+ -site | --site | --sit)
+ ac_prev=site ;;
+ -site=* | --site=* | --sit=*)
+ site="$ac_optarg" ;;
+
+ -srcdir | --srcdir | --srcdi | --srcd | --src | --sr)
+ ac_prev=srcdir ;;
+ -srcdir=* | --srcdir=* | --srcdi=* | --srcd=* | --src=* | --sr=*)
+ srcdir="$ac_optarg" ;;
+
+ -sysconfdir | --sysconfdir | --sysconfdi | --sysconfd | --sysconf \
+ | --syscon | --sysco | --sysc | --sys | --sy)
+ ac_prev=sysconfdir ;;
+ -sysconfdir=* | --sysconfdir=* | --sysconfdi=* | --sysconfd=* | --sysconf=* \
+ | --syscon=* | --sysco=* | --sysc=* | --sys=* | --sy=*)
+ sysconfdir="$ac_optarg" ;;
+
+ -target | --target | --targe | --targ | --tar | --ta | --t)
+ ac_prev=target ;;
+ -target=* | --target=* | --targe=* | --targ=* | --tar=* | --ta=* | --t=*)
+ target="$ac_optarg" ;;
+
+ -v | -verbose | --verbose | --verbos | --verbo | --verb)
+ verbose=yes ;;
+
+ -version | --version | --versio | --versi | --vers)
+ echo "configure generated by autoconf version 2.13"
+ exit 0 ;;
+
+ -with-* | --with-*)
+ ac_package=`echo $ac_option|sed -e 's/-*with-//' -e 's/=.*//'`
+ # Reject names that are not valid shell variable names.
+ if test -n "`echo $ac_package| sed 's/[-_a-zA-Z0-9]//g'`"; then
+ { echo "configure: error: $ac_package: invalid package name" 1>&2; exit 1; }
+ fi
+ ac_package=`echo $ac_package| sed 's/-/_/g'`
+ case "$ac_option" in
+ *=*) ;;
+ *) ac_optarg=yes ;;
+ esac
+ eval "with_${ac_package}='$ac_optarg'" ;;
+
+ -without-* | --without-*)
+ ac_package=`echo $ac_option|sed -e 's/-*without-//'`
+ # Reject names that are not valid shell variable names.
+ if test -n "`echo $ac_package| sed 's/[-a-zA-Z0-9_]//g'`"; then
+ { echo "configure: error: $ac_package: invalid package name" 1>&2; exit 1; }
+ fi
+ ac_package=`echo $ac_package| sed 's/-/_/g'`
+ eval "with_${ac_package}=no" ;;
+
+ --x)
+ # Obsolete; use --with-x.
+ with_x=yes ;;
+
+ -x-includes | --x-includes | --x-include | --x-includ | --x-inclu \
+ | --x-incl | --x-inc | --x-in | --x-i)
+ ac_prev=x_includes ;;
+ -x-includes=* | --x-includes=* | --x-include=* | --x-includ=* | --x-inclu=* \
+ | --x-incl=* | --x-inc=* | --x-in=* | --x-i=*)
+ x_includes="$ac_optarg" ;;
+
+ -x-libraries | --x-libraries | --x-librarie | --x-librari \
+ | --x-librar | --x-libra | --x-libr | --x-lib | --x-li | --x-l)
+ ac_prev=x_libraries ;;
+ -x-libraries=* | --x-libraries=* | --x-librarie=* | --x-librari=* \
+ | --x-librar=* | --x-libra=* | --x-libr=* | --x-lib=* | --x-li=* | --x-l=*)
+ x_libraries="$ac_optarg" ;;
+
+ -*) { echo "configure: error: $ac_option: invalid option; use --help to show usage" 1>&2; exit 1; }
+ ;;
+
+ *)
+ if test -n "`echo $ac_option| sed 's/[-a-z0-9.]//g'`"; then
+ echo "configure: warning: $ac_option: invalid host type" 1>&2
+ fi
+ if test "x$nonopt" != xNONE; then
+ { echo "configure: error: can only configure for one host and one target at a time" 1>&2; exit 1; }
+ fi
+ nonopt="$ac_option"
+ ;;
+
+ esac
+done
+
+if test -n "$ac_prev"; then
+ { echo "configure: error: missing argument to --`echo $ac_prev | sed 's/_/-/g'`" 1>&2; exit 1; }
+fi
+
+trap 'rm -fr conftest* confdefs* core core.* *.core $ac_clean_files; exit 1' 1 2 15
+
+# File descriptor usage:
+# 0 standard input
+# 1 file creation
+# 2 errors and warnings
+# 3 some systems may open it to /dev/tty
+# 4 used on the Kubota Titan
+# 6 checking for... messages and results
+# 5 compiler messages saved in config.log
+if test "$silent" = yes; then
+ exec 6>/dev/null
+else
+ exec 6>&1
+fi
+exec 5>./config.log
+
+echo "\
+This file contains any messages produced by compilers while
+running configure, to aid debugging if configure makes a mistake.
+" 1>&5
+
+# Strip out --no-create and --no-recursion so they do not pile up.
+# Also quote any args containing shell metacharacters.
+ac_configure_args=
+for ac_arg
+do
+ case "$ac_arg" in
+ -no-create | --no-create | --no-creat | --no-crea | --no-cre \
+ | --no-cr | --no-c) ;;
+ -no-recursion | --no-recursion | --no-recursio | --no-recursi \
+ | --no-recurs | --no-recur | --no-recu | --no-rec | --no-re | --no-r) ;;
+ *" "*|*" "*|*[\[\]\~\#\$\^\&\*\(\)\{\}\\\|\;\<\>\?]*)
+ ac_configure_args="$ac_configure_args '$ac_arg'" ;;
+ *) ac_configure_args="$ac_configure_args $ac_arg" ;;
+ esac
+done
+
+# NLS nuisances.
+# Only set these to C if already set. These must not be set unconditionally
+# because not all systems understand e.g. LANG=C (notably SCO).
+# Fixing LC_MESSAGES prevents Solaris sh from translating var values in `set'!
+# Non-C LC_CTYPE values break the ctype check.
+if test "${LANG+set}" = set; then LANG=C; export LANG; fi
+if test "${LC_ALL+set}" = set; then LC_ALL=C; export LC_ALL; fi
+if test "${LC_MESSAGES+set}" = set; then LC_MESSAGES=C; export LC_MESSAGES; fi
+if test "${LC_CTYPE+set}" = set; then LC_CTYPE=C; export LC_CTYPE; fi
+
+# confdefs.h avoids OS command line length limits that DEFS can exceed.
+rm -rf conftest* confdefs.h
+# AIX cpp loses on an empty file, so make sure it contains at least a newline.
+echo > confdefs.h
+
+# A filename unique to this package, relative to the directory that
+# configure is in, which we can look for to find out if srcdir is correct.
+ac_unique_file=embedding.cc
+
+# Find the source files, if location was not specified.
+if test -z "$srcdir"; then
+ ac_srcdir_defaulted=yes
+ # Try the directory containing this script, then its parent.
+ ac_prog=$0
+ ac_confdir=`echo $ac_prog|sed 's%/[^/][^/]*$%%'`
+ test "x$ac_confdir" = "x$ac_prog" && ac_confdir=.
+ srcdir=$ac_confdir
+ if test ! -r $srcdir/$ac_unique_file; then
+ srcdir=..
+ fi
+else
+ ac_srcdir_defaulted=no
+fi
+if test ! -r $srcdir/$ac_unique_file; then
+ if test "$ac_srcdir_defaulted" = yes; then
+ { echo "configure: error: can not find sources in $ac_confdir or .." 1>&2; exit 1; }
+ else
+ { echo "configure: error: can not find sources in $srcdir" 1>&2; exit 1; }
+ fi
+fi
+srcdir=`echo "${srcdir}" | sed 's%\([^/]\)/*$%\1%'`
+
+# Prefer explicitly selected file to automatically selected ones.
+if test -z "$CONFIG_SITE"; then
+ if test "x$prefix" != xNONE; then
+ CONFIG_SITE="$prefix/share/config.site $prefix/etc/config.site"
+ else
+ CONFIG_SITE="$ac_default_prefix/share/config.site $ac_default_prefix/etc/config.site"
+ fi
+fi
+for ac_site_file in $CONFIG_SITE; do
+ if test -r "$ac_site_file"; then
+ echo "loading site script $ac_site_file"
+ . "$ac_site_file"
+ fi
+done
+
+
+ac_ext=c
+# CFLAGS is not in ac_cpp because -g, -O, etc. are not valid cpp options.
+ac_cpp='$CPP $CPPFLAGS'
+ac_compile='${CC-cc} -c $CFLAGS $CPPFLAGS conftest.$ac_ext 1>&5'
+ac_link='${CC-cc} -o conftest${ac_exeext} $CFLAGS $CPPFLAGS $LDFLAGS conftest.$ac_ext $LIBS 1>&5'
+cross_compiling=$ac_cv_prog_cc_cross
+
+ac_exeext=
+ac_objext=o
+if (echo "testing\c"; echo 1,2,3) | grep c >/dev/null; then
+ # Stardent Vistra SVR4 grep lacks -e, says ghazi@caip.rutgers.edu.
+ if (echo -n testing; echo 1,2,3) | sed s/-n/xn/ | grep xn >/dev/null; then
+ ac_n= ac_c='
+' ac_t=' '
+ else
+ ac_n=-n ac_c= ac_t=
+ fi
+else
+ ac_n= ac_c='\c' ac_t=
+fi
+
+
+
+# Check whether --with-deal or --without-deal was given.
+if test "${with_deal+set}" = set; then
+ withval="$with_deal"
+ DEAL2_DIR=$withval
+else
+ if test -r ../../common/Make.global_options ; then
+ DEAL2_DIR=`pwd`/../.. ;
+ fi
+ if test -r ../../deal.II/common/Make.global_options ; then
+ DEAL2_DIR=`pwd`/../../deal.II ;
+ fi
+ if test -r ../../deal/common/Make.global_options ; then
+ DEAL2_DIR=`pwd`/../../deal ;
+ fi
+fi
+
+
+
+trap '' 1 2 15
+
+trap 'rm -fr conftest* confdefs* core core.* *.core $ac_clean_files; exit 1' 1 2 15
+
+test "x$prefix" = xNONE && prefix=$ac_default_prefix
+# Let make expand exec_prefix.
+test "x$exec_prefix" = xNONE && exec_prefix='${prefix}'
+
+# Any assignment to VPATH causes Sun make to only execute
+# the first set of double-colon rules, so remove it if not needed.
+# If there is a colon in the path, we need to keep it.
+if test "x$srcdir" = x.; then
+ ac_vpsub='/^[ ]*VPATH[ ]*=[^:]*$/d'
+fi
+
+trap 'rm -f $CONFIG_STATUS conftest*; exit 1' 1 2 15
+
+# Transform confdefs.h into DEFS.
+# Protect against shell expansion while executing Makefile rules.
+# Protect against Makefile macro expansion.
+cat > conftest.defs <<\EOF
+s%#define \([A-Za-z_][A-Za-z0-9_]*\) *\(.*\)%-D\1=\2%g
+s%[ `~#$^&*(){}\\|;'"<>?]%\\&%g
+s%\[%\\&%g
+s%\]%\\&%g
+s%\$%$$%g
+EOF
+DEFS=`sed -f conftest.defs confdefs.h | tr '\012' ' '`
+rm -f conftest.defs
+
+
+# Without the "./", some shells look in PATH for config.status.
+: ${CONFIG_STATUS=./config.status}
+
+echo creating $CONFIG_STATUS
+rm -f $CONFIG_STATUS
+cat > $CONFIG_STATUS <<EOF
+#! /bin/sh
+# Generated automatically by configure.
+# Run this file to recreate the current configuration.
+# This directory was configured as follows,
+# on host `(hostname || uname -n) 2>/dev/null | sed 1q`:
+#
+# $0 $ac_configure_args
+#
+# Compiler output produced by configure, useful for debugging
+# configure, is in ./config.log if it exists.
+
+ac_cs_usage="Usage: $CONFIG_STATUS [--recheck] [--version] [--help]"
+for ac_option
+do
+ case "\$ac_option" in
+ -recheck | --recheck | --rechec | --reche | --rech | --rec | --re | --r)
+ echo "running \${CONFIG_SHELL-/bin/sh} $0 $ac_configure_args --no-create --no-recursion"
+ exec \${CONFIG_SHELL-/bin/sh} $0 $ac_configure_args --no-create --no-recursion ;;
+ -version | --version | --versio | --versi | --vers | --ver | --ve | --v)
+ echo "$CONFIG_STATUS generated by autoconf version 2.13"
+ exit 0 ;;
+ -help | --help | --hel | --he | --h)
+ echo "\$ac_cs_usage"; exit 0 ;;
+ *) echo "\$ac_cs_usage"; exit 1 ;;
+ esac
+done
+
+ac_given_srcdir=$srcdir
+
+trap 'rm -fr `echo "Makefile" | sed "s/:[^ ]*//g"` conftest*; exit 1' 1 2 15
+EOF
+cat >> $CONFIG_STATUS <<EOF
+
+# Protect against being on the right side of a sed subst in config.status.
+sed 's/%@/@@/; s/@%/@@/; s/%g\$/@g/; /@g\$/s/[\\\\&%]/\\\\&/g;
+ s/@@/%@/; s/@@/@%/; s/@g\$/%g/' > conftest.subs <<\\CEOF
+$ac_vpsub
+$extrasub
+s%@SHELL@%$SHELL%g
+s%@CFLAGS@%$CFLAGS%g
+s%@CPPFLAGS@%$CPPFLAGS%g
+s%@CXXFLAGS@%$CXXFLAGS%g
+s%@FFLAGS@%$FFLAGS%g
+s%@DEFS@%$DEFS%g
+s%@LDFLAGS@%$LDFLAGS%g
+s%@LIBS@%$LIBS%g
+s%@exec_prefix@%$exec_prefix%g
+s%@prefix@%$prefix%g
+s%@program_transform_name@%$program_transform_name%g
+s%@bindir@%$bindir%g
+s%@sbindir@%$sbindir%g
+s%@libexecdir@%$libexecdir%g
+s%@datadir@%$datadir%g
+s%@sysconfdir@%$sysconfdir%g
+s%@sharedstatedir@%$sharedstatedir%g
+s%@localstatedir@%$localstatedir%g
+s%@libdir@%$libdir%g
+s%@includedir@%$includedir%g
+s%@oldincludedir@%$oldincludedir%g
+s%@infodir@%$infodir%g
+s%@mandir@%$mandir%g
+s%@DEAL2_DIR@%$DEAL2_DIR%g
+
+CEOF
+EOF
+
+cat >> $CONFIG_STATUS <<\EOF
+
+# Split the substitutions into bite-sized pieces for seds with
+# small command number limits, like on Digital OSF/1 and HP-UX.
+ac_max_sed_cmds=90 # Maximum number of lines to put in a sed script.
+ac_file=1 # Number of current file.
+ac_beg=1 # First line for current file.
+ac_end=$ac_max_sed_cmds # Line after last line for current file.
+ac_more_lines=:
+ac_sed_cmds=""
+while $ac_more_lines; do
+ if test $ac_beg -gt 1; then
+ sed "1,${ac_beg}d; ${ac_end}q" conftest.subs > conftest.s$ac_file
+ else
+ sed "${ac_end}q" conftest.subs > conftest.s$ac_file
+ fi
+ if test ! -s conftest.s$ac_file; then
+ ac_more_lines=false
+ rm -f conftest.s$ac_file
+ else
+ if test -z "$ac_sed_cmds"; then
+ ac_sed_cmds="sed -f conftest.s$ac_file"
+ else
+ ac_sed_cmds="$ac_sed_cmds | sed -f conftest.s$ac_file"
+ fi
+ ac_file=`expr $ac_file + 1`
+ ac_beg=$ac_end
+ ac_end=`expr $ac_end + $ac_max_sed_cmds`
+ fi
+done
+if test -z "$ac_sed_cmds"; then
+ ac_sed_cmds=cat
+fi
+EOF
+
+cat >> $CONFIG_STATUS <<EOF
+
+CONFIG_FILES=\${CONFIG_FILES-"Makefile"}
+EOF
+cat >> $CONFIG_STATUS <<\EOF
+for ac_file in .. $CONFIG_FILES; do if test "x$ac_file" != x..; then
+ # Support "outfile[:infile[:infile...]]", defaulting infile="outfile.in".
+ case "$ac_file" in
+ *:*) ac_file_in=`echo "$ac_file"|sed 's%[^:]*:%%'`
+ ac_file=`echo "$ac_file"|sed 's%:.*%%'` ;;
+ *) ac_file_in="${ac_file}.in" ;;
+ esac
+
+ # Adjust a relative srcdir, top_srcdir, and INSTALL for subdirectories.
+
+ # Remove last slash and all that follows it. Not all systems have dirname.
+ ac_dir=`echo $ac_file|sed 's%/[^/][^/]*$%%'`
+ if test "$ac_dir" != "$ac_file" && test "$ac_dir" != .; then
+ # The file is in a subdirectory.
+ test ! -d "$ac_dir" && mkdir "$ac_dir"
+ ac_dir_suffix="/`echo $ac_dir|sed 's%^\./%%'`"
+ # A "../" for each directory in $ac_dir_suffix.
+ ac_dots=`echo $ac_dir_suffix|sed 's%/[^/]*%../%g'`
+ else
+ ac_dir_suffix= ac_dots=
+ fi
+
+ case "$ac_given_srcdir" in
+ .) srcdir=.
+ if test -z "$ac_dots"; then top_srcdir=.
+ else top_srcdir=`echo $ac_dots|sed 's%/$%%'`; fi ;;
+ /*) srcdir="$ac_given_srcdir$ac_dir_suffix"; top_srcdir="$ac_given_srcdir" ;;
+ *) # Relative path.
+ srcdir="$ac_dots$ac_given_srcdir$ac_dir_suffix"
+ top_srcdir="$ac_dots$ac_given_srcdir" ;;
+ esac
+
+
+ echo creating "$ac_file"
+ rm -f "$ac_file"
+ configure_input="Generated automatically from `echo $ac_file_in|sed 's%.*/%%'` by configure."
+ case "$ac_file" in
+ *Makefile*) ac_comsub="1i\\
+# $configure_input" ;;
+ *) ac_comsub= ;;
+ esac
+
+ ac_file_inputs=`echo $ac_file_in|sed -e "s%^%$ac_given_srcdir/%" -e "s%:% $ac_given_srcdir/%g"`
+ sed -e "$ac_comsub
+s%@configure_input@%$configure_input%g
+s%@srcdir@%$srcdir%g
+s%@top_srcdir@%$top_srcdir%g
+" $ac_file_inputs | (eval "$ac_sed_cmds") > $ac_file
+fi; done
+rm -f conftest.s*
+
+EOF
+cat >> $CONFIG_STATUS <<EOF
+
+EOF
+cat >> $CONFIG_STATUS <<\EOF
+
+exit 0
+EOF
+chmod +x $CONFIG_STATUS
+rm -fr confdefs* $ac_clean_files
+test "$no_create" = yes || ${CONFIG_SHELL-/bin/sh} $CONFIG_STATUS || exit 1
+
--- /dev/null
+// $Id$
+// (c) Guido Kanschat
+//
+
+#include <base/quadrature_lib.h>
+#include <lac/vector.h>
+#include <lac/full_matrix.h>
+#include <lac/full_matrix.templates.h>
+#include <lac/solver_richardson.h>
+#include <lac/vector_memory.h>
+#include <grid/tria.h>
+#include <grid/tria_iterator.h>
+#include <dofs/dof_accessor.h>
+#include <grid/grid_generator.h>
+#include <fe/fe_q.h>
+#include <fe/fe_dgq.h>
+#include <fe/fe_system.h>
+#include <fe/mapping_q1.h>
+#include <fe/fe_values.h>
+#include <vector>
+#include <iomanip>
+#include <fstream>
+#include <strstream>
+#include <string>
+
+
+template <int dim>
+void compute_embedding (Triangulation<dim>& tr_coarse,
+ Triangulation<dim>& tr_fine,
+ const FiniteElement<dim>& fe_coarse,
+ const FiniteElement<dim>& fe_fine,
+ const char* name)
+{
+ cerr << name << '<' << dim << '>' << endl;
+
+ DoFHandler<dim> dof_coarse(tr_coarse);
+ dof_coarse.distribute_dofs(fe_coarse);
+ DoFHandler<dim> dof_fine(tr_fine);
+ dof_fine.distribute_dofs(fe_fine);
+
+ FullMatrix<double> A(dof_fine.n_dofs());
+ Vector<long double> f(dof_fine.n_dofs());
+ Vector<long double> u(dof_fine.n_dofs());
+ vector<vector<vector<double> > >
+ result(2<<dim,
+ vector<vector<double> >(fe_coarse.dofs_per_cell,
+ vector<double>(fe_fine.dofs_per_cell, 0.)));
+
+ DoFHandler<dim>::active_cell_iterator coarse_cell
+ = dof_coarse.begin_active();
+ vector<unsigned int> indices(fe_fine.dofs_per_cell);
+
+ MappingQ1<dim> mapping;
+ QGauss<dim> q_fine(12);
+
+ FEValues<dim> fine (mapping, fe_fine, q_fine,
+ update_q_points
+ | update_JxW_values
+ | update_values);
+
+
+ DoFHandler<dim>::active_cell_iterator c;
+
+ for (unsigned int coarse_no=0;coarse_no<fe_coarse.dofs_per_cell;
+ ++coarse_no)
+ {
+ A.clear();
+ f.clear();
+ u.clear();
+
+ for (c = dof_fine.begin_active();
+ c != dof_fine.end();
+ ++c)
+ {
+ c->get_dof_indices(indices);
+ fine.reinit(c);
+ // Build mass matrix and RHS
+
+ Quadrature<dim> q_coarse (fine.get_quadrature_points(),
+ fine.get_JxW_values());
+ FEValues<dim> coarse (mapping, fe_coarse, q_coarse,
+ update_values);
+ coarse.reinit(coarse_cell);
+
+ for (unsigned int k=0;k<fine.n_quadrature_points;++k)
+ {
+ double dx = fine.JxW(k);
+ for (unsigned int i=0;i<fe_fine.dofs_per_cell;++i)
+ {
+ double v = fine.shape_value(i,k);
+ f(indices[i]) += dx *
+ v * coarse.shape_value(coarse_no, k);
+ for (unsigned int j=0;j<fe_fine.dofs_per_cell;++j)
+ A(indices[i],indices[j]) += dx*v*fine.shape_value(j,k);
+ }
+ }
+ }
+// A.print_formatted(cout, 2, false, 4, "~", 9);
+ FullMatrix<double> P(A.n());
+ P.invert(A);
+ SolverControl control (100, 1.e-20, true, false);
+ PrimitiveVectorMemory<Vector<long double> > mem;
+ SolverRichardson<Vector<long double> > solver(control, mem);
+ solver.solve(A,u,f,P);
+
+ unsigned int cell_i=0;
+ for (c = dof_fine.begin_active();
+ c != dof_fine.end();
+ ++c, ++cell_i)
+ {
+ c->get_dof_indices(indices);
+ for (unsigned int i=0;i<fe_fine.dofs_per_cell;++i)
+ result[cell_i][coarse_no][i] = (fabs(u(indices[i])) > 1.e-16)
+ ? (27.*u(indices[i])) : 0.;
+ }
+ }
+
+ for (unsigned int cell=0;cell<tr_fine.n_active_cells();++cell)
+ {
+ cout << "static double[] "
+ << name << "_"
+ << dim << "d_"
+ << cell << " =" << endl << '{' << endl;
+ for (unsigned int i=0;i<fe_fine.dofs_per_cell;++i)
+ {
+ for (unsigned int j=0;j<fe_coarse.dofs_per_cell;++j)
+ {
+ cout << ' ' << setprecision(10) << result[cell][j][i] << "/27.,";
+ }
+ cout << endl;
+ }
+ cout << "};" << endl << endl;
+ }
+}
+
+#define PUSH_Q(k) FE_Q<dim> q ## k (k);\
+ elements.push_back (ElPair(&q ## k, "q" # k))
+
+#define PUSH_DGQ(k) FE_DGQ<dim> dgq ## k(k);\
+ elements.push_back (ElPair(&dgq ## k, "dgq" # k))
+
+template <int dim>
+void loop ()
+{
+ Triangulation<dim> tr_coarse;
+ Triangulation<dim> tr_fine;
+ GridGenerator::hyper_cube (tr_coarse, 0, 1);
+ GridGenerator::hyper_cube (tr_fine, 0, 1);
+ tr_fine.refine_global(1);
+
+ typedef pair<const FiniteElement<dim>*, const char*> ElPair ;
+ vector <ElPair> elements(0);
+
+ PUSH_DGQ(0);
+ PUSH_DGQ(1);
+ PUSH_DGQ(2);
+// PUSH_DGQ(3);
+ PUSH_DGQ(5);
+ PUSH_DGQ(6);
+ PUSH_DGQ(7);
+ // PUSH_Q(0);
+ PUSH_Q(1);
+ PUSH_Q(2);
+ PUSH_Q(3);
+ PUSH_Q(4);
+
+ char* name = new char[100];
+
+ // Embed all lower spaces into higher
+ unsigned int n = elements.size();
+ for (unsigned int i=0;i<n;++i)
+ for (unsigned int j=i;j<=i;++j)
+ {
+ ostrstream os (name, 99);
+ os << elements[i].second << "_into_"
+ << elements[j].second << "_refined" << ends;
+
+ compute_embedding (tr_coarse, tr_fine,
+ *elements[i].first,
+ *elements[j].first,
+ name);
+ }
+
+ delete [] name;
+}
+
+int main ()
+{
+ loop<1> ();
+ loop<2> ();
+ loop<3> ();
+}
--- /dev/null
+// $Id$
+// (c) Guido Kanschat
+//
+
+#include <base/quadrature_lib.h>
+#include <lac/vector.h>
+#include <lac/full_matrix.h>
+#include <lac/full_matrix.templates.h>
+#include <grid/tria.h>
+#include <grid/tria_iterator.h>
+#include <dofs/dof_accessor.h>
+#include <grid/grid_generator.h>
+#include <fe/fe_q.h>
+#include <fe/fe_dgq.h>
+#include <fe/fe_system.h>
+#include <fe/mapping_q1.h>
+#include <fe/fe_values.h>
+#include <vector>
+#include <iomanip>
+#include <fstream>
+#include <strstream>
+#include <string>
+
+
+template <int dim>
+inline bool in_cube (const Point<dim>& p)
+{
+ if ((p(0) < -1e-10) || (p(0) > 1.000000001))
+ return false;
+ if ((dim>1) && ((p(1) < -1e-10) || (p(1) > 1.000000001)))
+ return false;
+ if ((dim>2) && ((p(2) < -1e-10) || (p(2) > 1.000000001)))
+ return false;
+ return true;
+}
+
+template <int dim>
+inline void
+move_points (vector<Point<dim> >& v, unsigned int i)
+{
+ Point<dim> p;
+ switch (dim)
+ {
+ case 1:
+ if (i) p(0) = -1.;
+ break;
+ case 2:
+ if ((i==1) || (i==2))
+ p(0) = -1.;
+ if ((i==2) || (i==3))
+ p(1) = -1.;
+ break;
+ case 3:
+ if ((i==1) || (i==2) || (i==5) || (i==6))
+ p(0) = -1.;
+ if ((i==4) || (i==5) || (i==6) || (i==7))
+ p(1) = -1.;
+ if ((i==2) || (i==3) || (i==6) || (i==7))
+ p(2) = -1.;
+ break;
+ default:
+ Assert(false, ExcNotImplemented());
+ }
+ for (unsigned int j=0;j<v.size();++j)
+ v[j] += p;
+}
+
+
+template <int dim>
+void compute_interpolation (Triangulation<dim>& tr_coarse,
+ Triangulation<dim>& tr_fine,
+ const FiniteElement<dim>& fe_coarse,
+ const FiniteElement<dim>& fe_fine,
+ const char* name)
+{
+ DoFHandler<dim> dof_coarse(tr_coarse);
+ dof_coarse.distribute_dofs(fe_coarse);
+ DoFHandler<dim> dof_fine(tr_fine);
+ dof_fine.distribute_dofs(fe_fine);
+
+ vector<vector<vector<double> > >
+ result(2<<dim,
+ vector<vector<double> >(fe_coarse.dofs_per_cell,
+ vector<double>(fe_fine.dofs_per_cell, 0.)));
+
+ vector<unsigned int> indices(fe_fine.dofs_per_cell);
+
+ MappingQ1<dim> mapping;
+
+ vector<Point<dim> > points;
+ fe_coarse.get_unit_support_points (points);
+ vector<Point<dim> > q_points(points.size());
+ vector<double> dummy_weights(points.size(), 0.);
+
+ // Coarse level cell of fine grid
+ DoFHandler<dim>::cell_iterator father = dof_fine.begin (0);
+
+ for (unsigned int child_no = 0; child_no < GeometryInfo<dim>::children_per_cell;
+ ++child_no)
+ {
+ DoFHandler<dim>::active_cell_iterator c = father->child(child_no);
+
+ c->get_dof_indices(indices);
+ // Evaluate at support points
+ // of father cell.
+ for (unsigned int i=0;i<points.size();++i)
+ q_points[i] = 2.*points[i];
+ move_points (q_points, child_no);
+
+ Quadrature<dim> q_fine (q_points, dummy_weights);
+
+ FEValues<dim> fine (mapping, fe_fine, q_fine,
+ update_values);
+ fine.reinit(c);
+ // Build mass matrix and RHS
+
+ for (unsigned int k=0;k<fine.n_quadrature_points;++k)
+ {
+// cerr << k << '\t' << q_points[k];
+ if (in_cube(q_points[k]))
+ {
+ for (unsigned int i=0;i<fe_fine.dofs_per_cell;++i)
+ {
+ double v = fine.shape_value(i,k);
+// cerr << '[' << i << ' ' << v << ']';
+ result[child_no][k][i] = v;
+ }
+ }
+// cerr << endl;
+ }
+ }
+
+ for (unsigned int cell=0;cell<tr_fine.n_active_cells();++cell)
+ {
+ cout << "static double[] "
+ << name << "_"
+ << dim << "d_"
+ << cell << " =" << endl << '{' << endl;
+ for (unsigned int i=0;i<fe_coarse.dofs_per_cell;++i)
+ {
+ for (unsigned int j=0;j<fe_fine.dofs_per_cell;++j)
+ {
+ double a = result[cell][i][j];
+ if ((a<1.e-14) && (a>-1.e-14))
+ a = 0.;
+// cout << ' ' << setprecision(10) << a << ",";
+ if (a != 0.)
+ cout << "restriction[" << cell << "]("
+ << i << ',' << j << ") = "
+ << a << ';' << endl;
+ }
+ cout << endl;
+ }
+ cout << "};" << endl << endl;
+ }
+}
+
+
+template <int dim>
+void compute_projection (Triangulation<dim>& tr_coarse,
+ Triangulation<dim>& tr_fine,
+ const FiniteElement<dim>& fe_coarse,
+ const FiniteElement<dim>& fe_fine,
+ const char* name)
+{
+ DoFHandler<dim> dof_coarse(tr_coarse);
+ dof_coarse.distribute_dofs(fe_coarse);
+ DoFHandler<dim> dof_fine(tr_fine);
+ dof_fine.distribute_dofs(fe_fine);
+ DoFHandler<dim>::cell_iterator cell_coarse = dof_coarse.begin_active();
+
+ vector<vector<vector<double> > >
+ result(2<<dim,
+ vector<vector<double> >(fe_coarse.dofs_per_cell,
+ vector<double>(fe_fine.dofs_per_cell, 0.)));
+
+ vector<unsigned int> indices(fe_fine.dofs_per_cell);
+
+ MappingQ1<dim> mapping;
+ QGauss<dim> quadrature (9);
+ FullMatrix<double> mass (fe_coarse.dofs_per_cell);
+
+ FEValues<dim> coarse (mapping, fe_coarse, quadrature,
+ update_values
+ | update_JxW_values);
+ coarse.reinit (cell_coarse);
+
+ // Build coarse level mass matrix
+ for (unsigned int k=0;k<coarse.n_quadrature_points;++k)
+ {
+ for (unsigned int i=0;i<fe_coarse.dofs_per_cell;++i)
+ for (unsigned int j=0;j<fe_coarse.dofs_per_cell;++j)
+ mass(i,j) += coarse.JxW(k)
+ * coarse.shape_value(i,k) * coarse.shape_value(j,k);
+ }
+
+ FullMatrix<double> inverse (fe_coarse.dofs_per_cell);
+ inverse.invert(mass);
+ Vector<double> rhs (fe_coarse.dofs_per_cell);
+ Vector<double> u (fe_coarse.dofs_per_cell);
+
+ // Coarse level cell of fine grid
+ DoFHandler<dim>::cell_iterator father = dof_fine.begin (0);
+
+ FEValues<dim> fine (mapping, fe_fine, quadrature,
+ update_values
+ | update_JxW_values
+ | update_q_points);
+
+ // loop over all fine shape functions
+ for (unsigned int n=0;n<fe_fine.dofs_per_cell;++n)
+ {
+ rhs = 0.;
+ DoFHandler<dim>::active_cell_iterator c = father->child(0);
+
+ c->get_dof_indices(indices);
+
+ fine.reinit(c);
+
+ // Build a quadrature formula
+ // for the coarse cell
+ Quadrature<dim> q_coarse (fine.get_quadrature_points(),
+ fine.get_JxW_values());
+ FEValues<dim> coarse (mapping, fe_coarse, q_coarse,
+ update_values);
+ coarse.reinit(cell_coarse);
+ // Build RHS
+
+ for (unsigned int k=0;k<fine.n_quadrature_points;++k)
+ {
+ for (unsigned int i=0;i<fe_coarse.dofs_per_cell;++i)
+ {
+ double v = coarse.shape_value(i,k);
+ double f = fine.shape_value(n,k);
+ rhs(i) += fine.JxW(k) * v*f;
+ }
+ }
+ inverse.vmult(u,rhs);
+ for (unsigned int i=0;i<fe_coarse.dofs_per_cell;++i)
+ result[0][i][n] = u(i);
+ }
+
+ for (unsigned int cell=0;cell<1;++cell)
+ {
+ cout << "static double "
+ << name << "_"
+ << dim << "d[] =" << endl << '{' << endl;
+ for (unsigned int i=0;i<fe_coarse.dofs_per_cell;++i)
+ {
+ for (unsigned int j=0;j<fe_fine.dofs_per_cell;++j)
+ {
+ double a = result[cell][i][j];
+ if ((a<1.e-14) && (a>-1.e-14)) a = 0.;
+ cout << ' ' << setprecision(8) << a << ",";
+ }
+ cout << endl;
+ }
+ cout << "};" << endl << endl;
+ }
+}
+
+#define PUSH_Q(k) FE_Q<dim> q ## k(k);\
+ elements.push_back (ElPair(&q ## k, "dgq" # k))
+#define PUSH_DGQ(k) FE_DGQ<dim> dgq ## k(k);\
+ dgelements.push_back (ElPair(&dgq ## k, "dgq" # k))
+
+template <int dim>
+void loop ()
+{
+ Triangulation<dim> tr_coarse;
+ Triangulation<dim> tr_fine;
+ GridGenerator::hyper_cube (tr_coarse, 0, 1);
+ GridGenerator::hyper_cube (tr_fine, 0, 1);
+ tr_fine.refine_global(1);
+
+ typedef pair<const FiniteElement<dim>*, const char*> ElPair ;
+ vector <ElPair> elements(0);
+ vector <ElPair> dgelements(0);
+
+ PUSH_Q(1);
+ PUSH_Q(2);
+ PUSH_Q(3);
+ PUSH_Q(4);
+ PUSH_DGQ(0);
+ PUSH_DGQ(1);
+ PUSH_DGQ(2);
+ PUSH_DGQ(3);
+ PUSH_DGQ(4);
+ PUSH_DGQ(5);
+ PUSH_DGQ(6);
+ PUSH_DGQ(7);
+
+ char* name = new char[100];
+
+ unsigned int n = elements.size();
+ for (unsigned int i=0;i<n;++i)
+ for (unsigned int j=i;j<=i;++j)
+ {
+// ostrstream os (name, 99);
+// os << "interpolate " << elements[j].second << " refined onto "
+// << elements[i].second << ends;
+
+// compute_interpolation (tr_coarse, tr_fine,
+// *elements[i].first,
+// *elements[j].first,
+// name);
+ }
+
+ n = dgelements.size();
+ for (unsigned int i=0;i<n;++i)
+ for (unsigned int j=i;j<=i;++j)
+ {
+ ostrstream os (name, 99);
+ os << "project_" << dgelements[j].second << "_refined_onto_"
+ << dgelements[i].second << ends;
+
+ compute_projection (tr_coarse, tr_fine,
+ *dgelements[i].first,
+ *dgelements[j].first,
+ name);
+ }
+
+ delete [] name;
+}
+
+int main ()
+{
+ loop<1> ();
+ loop<2> ();
+ loop<3> ();
+}
projects / dealii-svn.git / commitdiff