]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
new utilities directory
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 5 Jun 2001 05:37:51 +0000 (05:37 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 5 Jun 2001 05:37:51 +0000 (05:37 +0000)
git-svn-id: https://svn.dealii.org/trunk@4741 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/contrib/configure
deal.II/contrib/configure.in
deal.II/contrib/utilities/.cvsignore [new file with mode: 0644]
deal.II/contrib/utilities/Makefile.in [new file with mode: 0644]
deal.II/contrib/utilities/configure [new file with mode: 0755]
deal.II/contrib/utilities/configure.in [new file with mode: 0644]
deal.II/contrib/utilities/embedding.cc [new file with mode: 0644]
deal.II/contrib/utilities/interpolation.cc [new file with mode: 0644]

index d0baa8ee378d9e20f9b50489c5b9c93be09003be..74ba8bdfb8591b10c637c14204e43273e1b1c500 100755 (executable)
@@ -527,6 +527,28 @@ cd $DEAL2_CONTRIB_DIR
 
 
 
+ac_aux_dir=
+for ac_dir in $srcdir $srcdir/.. $srcdir/../..; do
+  if test -f $ac_dir/install-sh; then
+    ac_aux_dir=$ac_dir
+    ac_install_sh="$ac_aux_dir/install-sh -c"
+    break
+  elif test -f $ac_dir/install.sh; then
+    ac_aux_dir=$ac_dir
+    ac_install_sh="$ac_aux_dir/install.sh -c"
+    break
+  fi
+done
+if test -z "$ac_aux_dir"; then
+  { echo "configure: error: can not find install-sh or install.sh in $srcdir $srcdir/.. $srcdir/../.." 1>&2; exit 1; }
+fi
+ac_config_guess=$ac_aux_dir/config.guess
+ac_config_sub=$ac_aux_dir/config.sub
+ac_configure=$ac_aux_dir/configure # This should be Cygnus configure.
+
+subdirs="utilities"
+
+
 
 
 trap '' 1 2 15
@@ -627,6 +649,7 @@ s%@includedir@%$includedir%g
 s%@oldincludedir@%$oldincludedir%g
 s%@infodir@%$infodir%g
 s%@mandir@%$mandir%g
+s%@subdirs@%$subdirs%g
 
 CEOF
 EOF
@@ -734,3 +757,96 @@ chmod +x $CONFIG_STATUS
 rm -fr confdefs* $ac_clean_files
 test "$no_create" = yes || ${CONFIG_SHELL-/bin/sh} $CONFIG_STATUS || exit 1
 
+if test "$no_recursion" != yes; then
+
+  # Remove --cache-file and --srcdir arguments so they do not pile up.
+  ac_sub_configure_args=
+  ac_prev=
+  for ac_arg in $ac_configure_args; do
+    if test -n "$ac_prev"; then
+      ac_prev=
+      continue
+    fi
+    case "$ac_arg" in
+    -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=*)
+      ;;
+    -srcdir | --srcdir | --srcdi | --srcd | --src | --sr)
+      ac_prev=srcdir ;;
+    -srcdir=* | --srcdir=* | --srcdi=* | --srcd=* | --src=* | --sr=*)
+      ;;
+    *) ac_sub_configure_args="$ac_sub_configure_args $ac_arg" ;;
+    esac
+  done
+
+  for ac_config_dir in utilities; do
+
+    # Do not complain, so a configure script can configure whichever
+    # parts of a large source tree are present.
+    if test ! -d $srcdir/$ac_config_dir; then
+      continue
+    fi
+
+    echo configuring in $ac_config_dir
+
+    case "$srcdir" in
+    .) ;;
+    *)
+      if test -d ./$ac_config_dir || mkdir ./$ac_config_dir; then :;
+      else
+        { echo "configure: error: can not create `pwd`/$ac_config_dir" 1>&2; exit 1; }
+      fi
+      ;;
+    esac
+
+    ac_popdir=`pwd`
+    cd $ac_config_dir
+
+      # A "../" for each directory in /$ac_config_dir.
+      ac_dots=`echo $ac_config_dir|sed -e 's%^\./%%' -e 's%[^/]$%&/%' -e 's%[^/]*/%../%g'`
+
+    case "$srcdir" in
+    .) # No --srcdir option.  We are building in place.
+      ac_sub_srcdir=$srcdir ;;
+    /*) # Absolute path.
+      ac_sub_srcdir=$srcdir/$ac_config_dir ;;
+    *) # Relative path.
+      ac_sub_srcdir=$ac_dots$srcdir/$ac_config_dir ;;
+    esac
+
+    # Check for guested configure; otherwise get Cygnus style configure.
+    if test -f $ac_sub_srcdir/configure; then
+      ac_sub_configure=$ac_sub_srcdir/configure
+    elif test -f $ac_sub_srcdir/configure.in; then
+      ac_sub_configure=$ac_configure
+    else
+      echo "configure: warning: no configuration information is in $ac_config_dir" 1>&2
+      ac_sub_configure=
+    fi
+
+    # The recursion is here.
+    if test -n "$ac_sub_configure"; then
+
+      # Make the cache file name correct relative to the subdirectory.
+      case "$cache_file" in
+      /*) ac_sub_cache_file=$cache_file ;;
+      *) # Relative path.
+        ac_sub_cache_file="$ac_dots$cache_file" ;;
+      esac
+
+      echo "running ${CONFIG_SHELL-/bin/sh} $ac_sub_configure $ac_sub_configure_args --cache-file=$ac_sub_cache_file --srcdir=$ac_sub_srcdir"
+      # The eval makes quoting arguments work.
+      if eval ${CONFIG_SHELL-/bin/sh} $ac_sub_configure $ac_sub_configure_args --cache-file=$ac_sub_cache_file --srcdir=$ac_sub_srcdir
+      then :
+      else
+        { echo "configure: error: $ac_sub_configure failed for $ac_config_dir" 1>&2; exit 1; }
+      fi
+    fi
+
+    cd $ac_popdir
+  done
+fi
+
index eedf13915a857b462337d49e6dbf52dbb0feaa99..b2ab0a4c7dfad0f10477b056d1dc305c67498a64 100644 (file)
@@ -27,6 +27,11 @@ cd kdoc/src
 cd $DEAL2_CONTRIB_DIR
 
 
+dnl -------------------------------------------------------------
+dnl Other subdirectories
+dnl -------------------------------------------------------------
+
+AC_CONFIG_SUBDIRS(utilities)
 
 
 dnl -------------------------------------------------------------
diff --git a/deal.II/contrib/utilities/.cvsignore b/deal.II/contrib/utilities/.cvsignore
new file mode 100644 (file)
index 0000000..b831c38
--- /dev/null
@@ -0,0 +1 @@
+config.status config.log Makefile
diff --git a/deal.II/contrib/utilities/Makefile.in b/deal.II/contrib/utilities/Makefile.in
new file mode 100644 (file)
index 0000000..72c1a97
--- /dev/null
@@ -0,0 +1,58 @@
+############################################################
+# $Id$
+# Copyright (C) 2000, 2001 by the deal.II authors
+############################################################
+
+############################################################
+# Include general settings for including DEAL libraries
+############################################################
+
+D = @DEAL2_DIR@
+include $D/common/Make.global_options
+debug-mode = on
+
+libraries = $(lib-deal2-1d.g) \
+            $(lib-deal2-2d.g) \
+            $(lib-deal2-3d.g) \
+            $(lib-lac.g)  \
+            $(lib-base.g)
+
+############################################################
+# First how to create executables, including all necessary
+# flags:
+############################################################
+
+flags     = $(CXXFLAGS.g)
+
+ifeq ($(findstring gcc,$(GXX_VERSION)),gcc)
+flags += -Wno-missing-noreturn
+endif
+
+%.o : %.cc Makefile
+       @echo =====optimized===== $<
+       @$(CXX) $(CXXFLAGS.o) -c $< -o $@
+%.g.o : %.cc Makefile
+       @echo =====debug========= $<
+       @$(CXX) $(CXXFLAGS.g) -c $< -o $@
+%.exe :
+       @echo =====linking======= $@
+       @$(CXX) $(LDFLAGS) -o $@ $^ $(LIBS)
+
+
+default: embedding.exe interpolation.exe
+
+############################################################
+
+embedding.exe    : embedding.o  $(libraries)
+interpolation.exe: interpolation.o $(libraries)
+
+############################################################
+# Automatic generation of dependencies
+############################################################
+
+Makefile.depend: $(shell echo *.cc)
+       @echo =====Dependencies== $@
+       @$(CXX) $(CXXFLAGS.g) $^ -M > $@
+       @perl -pi -e 's/(^[^.]+)\.o:/\1.o \1.g.o:/;' $@
+
+include Makefile.depend
diff --git a/deal.II/contrib/utilities/configure b/deal.II/contrib/utilities/configure
new file mode 100755 (executable)
index 0000000..359a936
--- /dev/null
@@ -0,0 +1,743 @@
+#! /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
+
diff --git a/deal.II/contrib/utilities/configure.in b/deal.II/contrib/utilities/configure.in
new file mode 100644 (file)
index 0000000..d35505a
--- /dev/null
@@ -0,0 +1,25 @@
+dnl
+dnl configure script to find the main deal library
+dnl in .. or ../deal.II or ../deal
+dnl
+
+define([AC_CACHE_LOAD], )dnl
+define([AC_CACHE_SAVE], )dnl
+AC_INIT(embedding.cc)
+
+dnl Supply deal-directory in command line or find it
+AC_ARG_WITH(deal,
+[  --with-deal=DIR where DIR is the root of the deal.II library],
+    DEAL2_DIR=$withval,
+    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)
+AC_SUBST(DEAL2_DIR)
+
+AC_OUTPUT(Makefile)
diff --git a/deal.II/contrib/utilities/embedding.cc b/deal.II/contrib/utilities/embedding.cc
new file mode 100644 (file)
index 0000000..f04b20b
--- /dev/null
@@ -0,0 +1,192 @@
+// $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> ();
+}
diff --git a/deal.II/contrib/utilities/interpolation.cc b/deal.II/contrib/utilities/interpolation.cc
new file mode 100644 (file)
index 0000000..9874f0d
--- /dev/null
@@ -0,0 +1,331 @@
+// $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> ();
+}

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.