1 /* Perform type resolution on the various stuctures.
2 Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008
3 Free Software Foundation, Inc.
4 Contributed by Andy Vaught
6 This file is part of GCC.
8 GCC is free software; you can redistribute it and/or modify it under
9 the terms of the GNU General Public License as published by the Free
10 Software Foundation; either version 3, or (at your option) any later
13 GCC is distributed in the hope that it will be useful, but WITHOUT ANY
14 WARRANTY; without even the implied warranty of MERCHANTABILITY or
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
18 You should have received a copy of the GNU General Public License
19 along with GCC; see the file COPYING3. If not see
20 <http://www.gnu.org/licenses/>. */
28 #include "arith.h" /* For gfc_compare_expr(). */
29 #include "dependency.h"
31 #include "target-memory.h" /* for gfc_simplify_transfer */
33 /* Types used in equivalence statements. */
37 SEQ_NONDEFAULT, SEQ_NUMERIC, SEQ_CHARACTER, SEQ_MIXED
41 /* Stack to keep track of the nesting of blocks as we move through the
42 code. See resolve_branch() and resolve_code(). */
44 typedef struct code_stack
46 struct gfc_code *head, *current, *tail;
47 struct code_stack *prev;
49 /* This bitmap keeps track of the targets valid for a branch from
51 bitmap reachable_labels;
55 static code_stack *cs_base = NULL;
58 /* Nonzero if we're inside a FORALL block. */
60 static int forall_flag;
62 /* Nonzero if we're inside a OpenMP WORKSHARE or PARALLEL WORKSHARE block. */
64 static int omp_workshare_flag;
66 /* Nonzero if we are processing a formal arglist. The corresponding function
67 resets the flag each time that it is read. */
68 static int formal_arg_flag = 0;
70 /* True if we are resolving a specification expression. */
71 static int specification_expr = 0;
73 /* The id of the last entry seen. */
74 static int current_entry_id;
76 /* We use bitmaps to determine if a branch target is valid. */
77 static bitmap_obstack labels_obstack;
80 gfc_is_formal_arg (void)
82 return formal_arg_flag;
85 /* Resolve types of formal argument lists. These have to be done early so that
86 the formal argument lists of module procedures can be copied to the
87 containing module before the individual procedures are resolved
88 individually. We also resolve argument lists of procedures in interface
89 blocks because they are self-contained scoping units.
91 Since a dummy argument cannot be a non-dummy procedure, the only
92 resort left for untyped names are the IMPLICIT types. */
95 resolve_formal_arglist (gfc_symbol *proc)
97 gfc_formal_arglist *f;
101 if (proc->result != NULL)
106 if (gfc_elemental (proc)
107 || sym->attr.pointer || sym->attr.allocatable
108 || (sym->as && sym->as->rank > 0))
110 proc->attr.always_explicit = 1;
111 sym->attr.always_explicit = 1;
116 for (f = proc->formal; f; f = f->next)
122 /* Alternate return placeholder. */
123 if (gfc_elemental (proc))
124 gfc_error ("Alternate return specifier in elemental subroutine "
125 "'%s' at %L is not allowed", proc->name,
127 if (proc->attr.function)
128 gfc_error ("Alternate return specifier in function "
129 "'%s' at %L is not allowed", proc->name,
134 if (sym->attr.if_source != IFSRC_UNKNOWN)
135 resolve_formal_arglist (sym);
137 if (sym->attr.subroutine || sym->attr.external || sym->attr.intrinsic)
139 if (gfc_pure (proc) && !gfc_pure (sym))
141 gfc_error ("Dummy procedure '%s' of PURE procedure at %L must "
142 "also be PURE", sym->name, &sym->declared_at);
146 if (gfc_elemental (proc))
148 gfc_error ("Dummy procedure at %L not allowed in ELEMENTAL "
149 "procedure", &sym->declared_at);
153 if (sym->attr.function
154 && sym->ts.type == BT_UNKNOWN
155 && sym->attr.intrinsic)
157 gfc_intrinsic_sym *isym;
158 isym = gfc_find_function (sym->name);
159 if (isym == NULL || !isym->specific)
161 gfc_error ("Unable to find a specific INTRINSIC procedure "
162 "for the reference '%s' at %L", sym->name,
171 if (sym->ts.type == BT_UNKNOWN)
173 if (!sym->attr.function || sym->result == sym)
174 gfc_set_default_type (sym, 1, sym->ns);
177 gfc_resolve_array_spec (sym->as, 0);
179 /* We can't tell if an array with dimension (:) is assumed or deferred
180 shape until we know if it has the pointer or allocatable attributes.
182 if (sym->as && sym->as->rank > 0 && sym->as->type == AS_DEFERRED
183 && !(sym->attr.pointer || sym->attr.allocatable))
185 sym->as->type = AS_ASSUMED_SHAPE;
186 for (i = 0; i < sym->as->rank; i++)
187 sym->as->lower[i] = gfc_int_expr (1);
190 if ((sym->as && sym->as->rank > 0 && sym->as->type == AS_ASSUMED_SHAPE)
191 || sym->attr.pointer || sym->attr.allocatable || sym->attr.target
192 || sym->attr.optional)
194 proc->attr.always_explicit = 1;
196 proc->result->attr.always_explicit = 1;
199 /* If the flavor is unknown at this point, it has to be a variable.
200 A procedure specification would have already set the type. */
202 if (sym->attr.flavor == FL_UNKNOWN)
203 gfc_add_flavor (&sym->attr, FL_VARIABLE, sym->name, &sym->declared_at);
205 if (gfc_pure (proc) && !sym->attr.pointer
206 && sym->attr.flavor != FL_PROCEDURE)
208 if (proc->attr.function && sym->attr.intent != INTENT_IN)
209 gfc_error ("Argument '%s' of pure function '%s' at %L must be "
210 "INTENT(IN)", sym->name, proc->name,
213 if (proc->attr.subroutine && sym->attr.intent == INTENT_UNKNOWN)
214 gfc_error ("Argument '%s' of pure subroutine '%s' at %L must "
215 "have its INTENT specified", sym->name, proc->name,
219 if (gfc_elemental (proc))
223 gfc_error ("Argument '%s' of elemental procedure at %L must "
224 "be scalar", sym->name, &sym->declared_at);
228 if (sym->attr.pointer)
230 gfc_error ("Argument '%s' of elemental procedure at %L cannot "
231 "have the POINTER attribute", sym->name,
236 if (sym->attr.flavor == FL_PROCEDURE)
238 gfc_error ("Dummy procedure '%s' not allowed in elemental "
239 "procedure '%s' at %L", sym->name, proc->name,
245 /* Each dummy shall be specified to be scalar. */
246 if (proc->attr.proc == PROC_ST_FUNCTION)
250 gfc_error ("Argument '%s' of statement function at %L must "
251 "be scalar", sym->name, &sym->declared_at);
255 if (sym->ts.type == BT_CHARACTER)
257 gfc_charlen *cl = sym->ts.cl;
258 if (!cl || !cl->length || cl->length->expr_type != EXPR_CONSTANT)
260 gfc_error ("Character-valued argument '%s' of statement "
261 "function at %L must have constant length",
262 sym->name, &sym->declared_at);
272 /* Work function called when searching for symbols that have argument lists
273 associated with them. */
276 find_arglists (gfc_symbol *sym)
278 if (sym->attr.if_source == IFSRC_UNKNOWN || sym->ns != gfc_current_ns)
281 resolve_formal_arglist (sym);
285 /* Given a namespace, resolve all formal argument lists within the namespace.
289 resolve_formal_arglists (gfc_namespace *ns)
294 gfc_traverse_ns (ns, find_arglists);
299 resolve_contained_fntype (gfc_symbol *sym, gfc_namespace *ns)
303 /* If this namespace is not a function or an entry master function,
305 if (! sym || !(sym->attr.function || sym->attr.flavor == FL_VARIABLE)
306 || sym->attr.entry_master)
309 /* Try to find out of what the return type is. */
310 if (sym->result->ts.type == BT_UNKNOWN)
312 t = gfc_set_default_type (sym->result, 0, ns);
314 if (t == FAILURE && !sym->result->attr.untyped)
316 if (sym->result == sym)
317 gfc_error ("Contained function '%s' at %L has no IMPLICIT type",
318 sym->name, &sym->declared_at);
320 gfc_error ("Result '%s' of contained function '%s' at %L has "
321 "no IMPLICIT type", sym->result->name, sym->name,
322 &sym->result->declared_at);
323 sym->result->attr.untyped = 1;
327 /* Fortran 95 Draft Standard, page 51, Section 5.1.1.5, on the Character
328 type, lists the only ways a character length value of * can be used:
329 dummy arguments of procedures, named constants, and function results
330 in external functions. Internal function results are not on that list;
331 ergo, not permitted. */
333 if (sym->result->ts.type == BT_CHARACTER)
335 gfc_charlen *cl = sym->result->ts.cl;
336 if (!cl || !cl->length)
337 gfc_error ("Character-valued internal function '%s' at %L must "
338 "not be assumed length", sym->name, &sym->declared_at);
343 /* Add NEW_ARGS to the formal argument list of PROC, taking care not to
344 introduce duplicates. */
347 merge_argument_lists (gfc_symbol *proc, gfc_formal_arglist *new_args)
349 gfc_formal_arglist *f, *new_arglist;
352 for (; new_args != NULL; new_args = new_args->next)
354 new_sym = new_args->sym;
355 /* See if this arg is already in the formal argument list. */
356 for (f = proc->formal; f; f = f->next)
358 if (new_sym == f->sym)
365 /* Add a new argument. Argument order is not important. */
366 new_arglist = gfc_get_formal_arglist ();
367 new_arglist->sym = new_sym;
368 new_arglist->next = proc->formal;
369 proc->formal = new_arglist;
374 /* Flag the arguments that are not present in all entries. */
377 check_argument_lists (gfc_symbol *proc, gfc_formal_arglist *new_args)
379 gfc_formal_arglist *f, *head;
382 for (f = proc->formal; f; f = f->next)
387 for (new_args = head; new_args; new_args = new_args->next)
389 if (new_args->sym == f->sym)
396 f->sym->attr.not_always_present = 1;
401 /* Resolve alternate entry points. If a symbol has multiple entry points we
402 create a new master symbol for the main routine, and turn the existing
403 symbol into an entry point. */
406 resolve_entries (gfc_namespace *ns)
408 gfc_namespace *old_ns;
412 char name[GFC_MAX_SYMBOL_LEN + 1];
413 static int master_count = 0;
415 if (ns->proc_name == NULL)
418 /* No need to do anything if this procedure doesn't have alternate entry
423 /* We may already have resolved alternate entry points. */
424 if (ns->proc_name->attr.entry_master)
427 /* If this isn't a procedure something has gone horribly wrong. */
428 gcc_assert (ns->proc_name->attr.flavor == FL_PROCEDURE);
430 /* Remember the current namespace. */
431 old_ns = gfc_current_ns;
435 /* Add the main entry point to the list of entry points. */
436 el = gfc_get_entry_list ();
437 el->sym = ns->proc_name;
439 el->next = ns->entries;
441 ns->proc_name->attr.entry = 1;
443 /* If it is a module function, it needs to be in the right namespace
444 so that gfc_get_fake_result_decl can gather up the results. The
445 need for this arose in get_proc_name, where these beasts were
446 left in their own namespace, to keep prior references linked to
447 the entry declaration.*/
448 if (ns->proc_name->attr.function
449 && ns->parent && ns->parent->proc_name->attr.flavor == FL_MODULE)
452 /* Do the same for entries where the master is not a module
453 procedure. These are retained in the module namespace because
454 of the module procedure declaration. */
455 for (el = el->next; el; el = el->next)
456 if (el->sym->ns->proc_name->attr.flavor == FL_MODULE
457 && el->sym->attr.mod_proc)
461 /* Add an entry statement for it. */
468 /* Create a new symbol for the master function. */
469 /* Give the internal function a unique name (within this file).
470 Also include the function name so the user has some hope of figuring
471 out what is going on. */
472 snprintf (name, GFC_MAX_SYMBOL_LEN, "master.%d.%s",
473 master_count++, ns->proc_name->name);
474 gfc_get_ha_symbol (name, &proc);
475 gcc_assert (proc != NULL);
477 gfc_add_procedure (&proc->attr, PROC_INTERNAL, proc->name, NULL);
478 if (ns->proc_name->attr.subroutine)
479 gfc_add_subroutine (&proc->attr, proc->name, NULL);
483 gfc_typespec *ts, *fts;
484 gfc_array_spec *as, *fas;
485 gfc_add_function (&proc->attr, proc->name, NULL);
487 fas = ns->entries->sym->as;
488 fas = fas ? fas : ns->entries->sym->result->as;
489 fts = &ns->entries->sym->result->ts;
490 if (fts->type == BT_UNKNOWN)
491 fts = gfc_get_default_type (ns->entries->sym->result, NULL);
492 for (el = ns->entries->next; el; el = el->next)
494 ts = &el->sym->result->ts;
496 as = as ? as : el->sym->result->as;
497 if (ts->type == BT_UNKNOWN)
498 ts = gfc_get_default_type (el->sym->result, NULL);
500 if (! gfc_compare_types (ts, fts)
501 || (el->sym->result->attr.dimension
502 != ns->entries->sym->result->attr.dimension)
503 || (el->sym->result->attr.pointer
504 != ns->entries->sym->result->attr.pointer))
506 else if (as && fas && ns->entries->sym->result != el->sym->result
507 && gfc_compare_array_spec (as, fas) == 0)
508 gfc_error ("Function %s at %L has entries with mismatched "
509 "array specifications", ns->entries->sym->name,
510 &ns->entries->sym->declared_at);
511 /* The characteristics need to match and thus both need to have
512 the same string length, i.e. both len=*, or both len=4.
513 Having both len=<variable> is also possible, but difficult to
514 check at compile time. */
515 else if (ts->type == BT_CHARACTER && ts->cl && fts->cl
516 && (((ts->cl->length && !fts->cl->length)
517 ||(!ts->cl->length && fts->cl->length))
519 && ts->cl->length->expr_type
520 != fts->cl->length->expr_type)
522 && ts->cl->length->expr_type == EXPR_CONSTANT
523 && mpz_cmp (ts->cl->length->value.integer,
524 fts->cl->length->value.integer) != 0)))
525 gfc_notify_std (GFC_STD_GNU, "Extension: Function %s at %L with "
526 "entries returning variables of different "
527 "string lengths", ns->entries->sym->name,
528 &ns->entries->sym->declared_at);
533 sym = ns->entries->sym->result;
534 /* All result types the same. */
536 if (sym->attr.dimension)
537 gfc_set_array_spec (proc, gfc_copy_array_spec (sym->as), NULL);
538 if (sym->attr.pointer)
539 gfc_add_pointer (&proc->attr, NULL);
543 /* Otherwise the result will be passed through a union by
545 proc->attr.mixed_entry_master = 1;
546 for (el = ns->entries; el; el = el->next)
548 sym = el->sym->result;
549 if (sym->attr.dimension)
551 if (el == ns->entries)
552 gfc_error ("FUNCTION result %s can't be an array in "
553 "FUNCTION %s at %L", sym->name,
554 ns->entries->sym->name, &sym->declared_at);
556 gfc_error ("ENTRY result %s can't be an array in "
557 "FUNCTION %s at %L", sym->name,
558 ns->entries->sym->name, &sym->declared_at);
560 else if (sym->attr.pointer)
562 if (el == ns->entries)
563 gfc_error ("FUNCTION result %s can't be a POINTER in "
564 "FUNCTION %s at %L", sym->name,
565 ns->entries->sym->name, &sym->declared_at);
567 gfc_error ("ENTRY result %s can't be a POINTER in "
568 "FUNCTION %s at %L", sym->name,
569 ns->entries->sym->name, &sym->declared_at);
574 if (ts->type == BT_UNKNOWN)
575 ts = gfc_get_default_type (sym, NULL);
579 if (ts->kind == gfc_default_integer_kind)
583 if (ts->kind == gfc_default_real_kind
584 || ts->kind == gfc_default_double_kind)
588 if (ts->kind == gfc_default_complex_kind)
592 if (ts->kind == gfc_default_logical_kind)
596 /* We will issue error elsewhere. */
604 if (el == ns->entries)
605 gfc_error ("FUNCTION result %s can't be of type %s "
606 "in FUNCTION %s at %L", sym->name,
607 gfc_typename (ts), ns->entries->sym->name,
610 gfc_error ("ENTRY result %s can't be of type %s "
611 "in FUNCTION %s at %L", sym->name,
612 gfc_typename (ts), ns->entries->sym->name,
619 proc->attr.access = ACCESS_PRIVATE;
620 proc->attr.entry_master = 1;
622 /* Merge all the entry point arguments. */
623 for (el = ns->entries; el; el = el->next)
624 merge_argument_lists (proc, el->sym->formal);
626 /* Check the master formal arguments for any that are not
627 present in all entry points. */
628 for (el = ns->entries; el; el = el->next)
629 check_argument_lists (proc, el->sym->formal);
631 /* Use the master function for the function body. */
632 ns->proc_name = proc;
634 /* Finalize the new symbols. */
635 gfc_commit_symbols ();
637 /* Restore the original namespace. */
638 gfc_current_ns = old_ns;
643 has_default_initializer (gfc_symbol *der)
647 gcc_assert (der->attr.flavor == FL_DERIVED);
648 for (c = der->components; c; c = c->next)
649 if ((c->ts.type != BT_DERIVED && c->initializer)
650 || (c->ts.type == BT_DERIVED
651 && (!c->pointer && has_default_initializer (c->ts.derived))))
657 /* Resolve common variables. */
659 resolve_common_vars (gfc_symbol *sym, bool named_common)
661 gfc_symbol *csym = sym;
663 for (; csym; csym = csym->common_next)
665 if (csym->value || csym->attr.data)
667 if (!csym->ns->is_block_data)
668 gfc_notify_std (GFC_STD_GNU, "Variable '%s' at %L is in COMMON "
669 "but only in BLOCK DATA initialization is "
670 "allowed", csym->name, &csym->declared_at);
671 else if (!named_common)
672 gfc_notify_std (GFC_STD_GNU, "Initialized variable '%s' at %L is "
673 "in a blank COMMON but initialization is only "
674 "allowed in named common blocks", csym->name,
678 if (csym->ts.type != BT_DERIVED)
681 if (!(csym->ts.derived->attr.sequence
682 || csym->ts.derived->attr.is_bind_c))
683 gfc_error_now ("Derived type variable '%s' in COMMON at %L "
684 "has neither the SEQUENCE nor the BIND(C) "
685 "attribute", csym->name, &csym->declared_at);
686 if (csym->ts.derived->attr.alloc_comp)
687 gfc_error_now ("Derived type variable '%s' in COMMON at %L "
688 "has an ultimate component that is "
689 "allocatable", csym->name, &csym->declared_at);
690 if (has_default_initializer (csym->ts.derived))
691 gfc_error_now ("Derived type variable '%s' in COMMON at %L "
692 "may not have default initializer", csym->name,
697 /* Resolve common blocks. */
699 resolve_common_blocks (gfc_symtree *common_root)
703 if (common_root == NULL)
706 if (common_root->left)
707 resolve_common_blocks (common_root->left);
708 if (common_root->right)
709 resolve_common_blocks (common_root->right);
711 resolve_common_vars (common_root->n.common->head, true);
713 gfc_find_symbol (common_root->name, gfc_current_ns, 0, &sym);
717 if (sym->attr.flavor == FL_PARAMETER)
718 gfc_error ("COMMON block '%s' at %L is used as PARAMETER at %L",
719 sym->name, &common_root->n.common->where, &sym->declared_at);
721 if (sym->attr.intrinsic)
722 gfc_error ("COMMON block '%s' at %L is also an intrinsic procedure",
723 sym->name, &common_root->n.common->where);
724 else if (sym->attr.result
725 ||(sym->attr.function && gfc_current_ns->proc_name == sym))
726 gfc_notify_std (GFC_STD_F2003, "Fortran 2003: COMMON block '%s' at %L "
727 "that is also a function result", sym->name,
728 &common_root->n.common->where);
729 else if (sym->attr.flavor == FL_PROCEDURE && sym->attr.proc != PROC_INTERNAL
730 && sym->attr.proc != PROC_ST_FUNCTION)
731 gfc_notify_std (GFC_STD_F2003, "Fortran 2003: COMMON block '%s' at %L "
732 "that is also a global procedure", sym->name,
733 &common_root->n.common->where);
737 /* Resolve contained function types. Because contained functions can call one
738 another, they have to be worked out before any of the contained procedures
741 The good news is that if a function doesn't already have a type, the only
742 way it can get one is through an IMPLICIT type or a RESULT variable, because
743 by definition contained functions are contained namespace they're contained
744 in, not in a sibling or parent namespace. */
747 resolve_contained_functions (gfc_namespace *ns)
749 gfc_namespace *child;
752 resolve_formal_arglists (ns);
754 for (child = ns->contained; child; child = child->sibling)
756 /* Resolve alternate entry points first. */
757 resolve_entries (child);
759 /* Then check function return types. */
760 resolve_contained_fntype (child->proc_name, child);
761 for (el = child->entries; el; el = el->next)
762 resolve_contained_fntype (el->sym, child);
767 /* Resolve all of the elements of a structure constructor and make sure that
768 the types are correct. */
771 resolve_structure_cons (gfc_expr *expr)
773 gfc_constructor *cons;
779 cons = expr->value.constructor;
780 /* A constructor may have references if it is the result of substituting a
781 parameter variable. In this case we just pull out the component we
784 comp = expr->ref->u.c.sym->components;
786 comp = expr->ts.derived->components;
788 /* See if the user is trying to invoke a structure constructor for one of
789 the iso_c_binding derived types. */
790 if (expr->ts.derived && expr->ts.derived->ts.is_iso_c && cons
791 && cons->expr != NULL)
793 gfc_error ("Components of structure constructor '%s' at %L are PRIVATE",
794 expr->ts.derived->name, &(expr->where));
798 for (; comp; comp = comp->next, cons = cons->next)
805 if (gfc_resolve_expr (cons->expr) == FAILURE)
811 rank = comp->as ? comp->as->rank : 0;
812 if (cons->expr->expr_type != EXPR_NULL && rank != cons->expr->rank
813 && (comp->allocatable || cons->expr->rank))
815 gfc_error ("The rank of the element in the derived type "
816 "constructor at %L does not match that of the "
817 "component (%d/%d)", &cons->expr->where,
818 cons->expr->rank, rank);
822 /* If we don't have the right type, try to convert it. */
824 if (!gfc_compare_types (&cons->expr->ts, &comp->ts))
827 if (comp->pointer && cons->expr->ts.type != BT_UNKNOWN)
828 gfc_error ("The element in the derived type constructor at %L, "
829 "for pointer component '%s', is %s but should be %s",
830 &cons->expr->where, comp->name,
831 gfc_basic_typename (cons->expr->ts.type),
832 gfc_basic_typename (comp->ts.type));
834 t = gfc_convert_type (cons->expr, &comp->ts, 1);
837 if (cons->expr->expr_type == EXPR_NULL
838 && !(comp->pointer || comp->allocatable))
841 gfc_error ("The NULL in the derived type constructor at %L is "
842 "being applied to component '%s', which is neither "
843 "a POINTER nor ALLOCATABLE", &cons->expr->where,
847 if (!comp->pointer || cons->expr->expr_type == EXPR_NULL)
850 a = gfc_expr_attr (cons->expr);
852 if (!a.pointer && !a.target)
855 gfc_error ("The element in the derived type constructor at %L, "
856 "for pointer component '%s' should be a POINTER or "
857 "a TARGET", &cons->expr->where, comp->name);
865 /****************** Expression name resolution ******************/
867 /* Returns 0 if a symbol was not declared with a type or
868 attribute declaration statement, nonzero otherwise. */
871 was_declared (gfc_symbol *sym)
877 if (!a.implicit_type && sym->ts.type != BT_UNKNOWN)
880 if (a.allocatable || a.dimension || a.dummy || a.external || a.intrinsic
881 || a.optional || a.pointer || a.save || a.target || a.volatile_
882 || a.value || a.access != ACCESS_UNKNOWN || a.intent != INTENT_UNKNOWN)
889 /* Determine if a symbol is generic or not. */
892 generic_sym (gfc_symbol *sym)
896 if (sym->attr.generic ||
897 (sym->attr.intrinsic && gfc_generic_intrinsic (sym->name)))
900 if (was_declared (sym) || sym->ns->parent == NULL)
903 gfc_find_symbol (sym->name, sym->ns->parent, 1, &s);
910 return generic_sym (s);
917 /* Determine if a symbol is specific or not. */
920 specific_sym (gfc_symbol *sym)
924 if (sym->attr.if_source == IFSRC_IFBODY
925 || sym->attr.proc == PROC_MODULE
926 || sym->attr.proc == PROC_INTERNAL
927 || sym->attr.proc == PROC_ST_FUNCTION
928 || (sym->attr.intrinsic && gfc_specific_intrinsic (sym->name))
929 || sym->attr.external)
932 if (was_declared (sym) || sym->ns->parent == NULL)
935 gfc_find_symbol (sym->name, sym->ns->parent, 1, &s);
937 return (s == NULL) ? 0 : specific_sym (s);
941 /* Figure out if the procedure is specific, generic or unknown. */
944 { PTYPE_GENERIC = 1, PTYPE_SPECIFIC, PTYPE_UNKNOWN }
948 procedure_kind (gfc_symbol *sym)
950 if (generic_sym (sym))
951 return PTYPE_GENERIC;
953 if (specific_sym (sym))
954 return PTYPE_SPECIFIC;
956 return PTYPE_UNKNOWN;
959 /* Check references to assumed size arrays. The flag need_full_assumed_size
960 is nonzero when matching actual arguments. */
962 static int need_full_assumed_size = 0;
965 check_assumed_size_reference (gfc_symbol *sym, gfc_expr *e)
967 if (need_full_assumed_size || !(sym->as && sym->as->type == AS_ASSUMED_SIZE))
970 if ((e->ref->u.ar.end[e->ref->u.ar.as->rank - 1] == NULL)
971 && (e->ref->u.ar.as->type == AS_ASSUMED_SIZE)
972 && (e->ref->u.ar.type == DIMEN_ELEMENT))
974 gfc_error ("The upper bound in the last dimension must "
975 "appear in the reference to the assumed size "
976 "array '%s' at %L", sym->name, &e->where);
983 /* Look for bad assumed size array references in argument expressions
984 of elemental and array valued intrinsic procedures. Since this is
985 called from procedure resolution functions, it only recurses at
989 resolve_assumed_size_actual (gfc_expr *e)
994 switch (e->expr_type)
997 if (e->symtree && check_assumed_size_reference (e->symtree->n.sym, e))
1002 if (resolve_assumed_size_actual (e->value.op.op1)
1003 || resolve_assumed_size_actual (e->value.op.op2))
1014 /* Resolve an actual argument list. Most of the time, this is just
1015 resolving the expressions in the list.
1016 The exception is that we sometimes have to decide whether arguments
1017 that look like procedure arguments are really simple variable
1021 resolve_actual_arglist (gfc_actual_arglist *arg, procedure_type ptype)
1024 gfc_symtree *parent_st;
1026 int save_need_full_assumed_size;
1028 for (; arg; arg = arg->next)
1033 /* Check the label is a valid branching target. */
1036 if (arg->label->defined == ST_LABEL_UNKNOWN)
1038 gfc_error ("Label %d referenced at %L is never defined",
1039 arg->label->value, &arg->label->where);
1046 if (e->expr_type == FL_VARIABLE && e->symtree->ambiguous)
1048 gfc_error ("'%s' at %L is ambiguous", e->symtree->n.sym->name,
1053 if (e->ts.type != BT_PROCEDURE)
1055 save_need_full_assumed_size = need_full_assumed_size;
1056 if (e->expr_type != FL_VARIABLE)
1057 need_full_assumed_size = 0;
1058 if (gfc_resolve_expr (e) != SUCCESS)
1060 need_full_assumed_size = save_need_full_assumed_size;
1064 /* See if the expression node should really be a variable reference. */
1066 sym = e->symtree->n.sym;
1068 if (sym->attr.flavor == FL_PROCEDURE
1069 || sym->attr.intrinsic
1070 || sym->attr.external)
1074 /* If a procedure is not already determined to be something else
1075 check if it is intrinsic. */
1076 if (!sym->attr.intrinsic
1077 && !(sym->attr.external || sym->attr.use_assoc
1078 || sym->attr.if_source == IFSRC_IFBODY)
1079 && gfc_intrinsic_name (sym->name, sym->attr.subroutine))
1080 sym->attr.intrinsic = 1;
1082 if (sym->attr.proc == PROC_ST_FUNCTION)
1084 gfc_error ("Statement function '%s' at %L is not allowed as an "
1085 "actual argument", sym->name, &e->where);
1088 actual_ok = gfc_intrinsic_actual_ok (sym->name,
1089 sym->attr.subroutine);
1090 if (sym->attr.intrinsic && actual_ok == 0)
1092 gfc_error ("Intrinsic '%s' at %L is not allowed as an "
1093 "actual argument", sym->name, &e->where);
1096 if (sym->attr.contained && !sym->attr.use_assoc
1097 && sym->ns->proc_name->attr.flavor != FL_MODULE)
1099 gfc_error ("Internal procedure '%s' is not allowed as an "
1100 "actual argument at %L", sym->name, &e->where);
1103 if (sym->attr.elemental && !sym->attr.intrinsic)
1105 gfc_error ("ELEMENTAL non-INTRINSIC procedure '%s' is not "
1106 "allowed as an actual argument at %L", sym->name,
1110 /* Check if a generic interface has a specific procedure
1111 with the same name before emitting an error. */
1112 if (sym->attr.generic)
1115 for (p = sym->generic; p; p = p->next)
1116 if (strcmp (sym->name, p->sym->name) == 0)
1118 e->symtree = gfc_find_symtree
1119 (p->sym->ns->sym_root, sym->name);
1124 if (p == NULL || e->symtree == NULL)
1125 gfc_error ("GENERIC procedure '%s' is not "
1126 "allowed as an actual argument at %L", sym->name,
1130 /* If the symbol is the function that names the current (or
1131 parent) scope, then we really have a variable reference. */
1133 if (sym->attr.function && sym->result == sym
1134 && (sym->ns->proc_name == sym
1135 || (sym->ns->parent != NULL
1136 && sym->ns->parent->proc_name == sym)))
1139 /* If all else fails, see if we have a specific intrinsic. */
1140 if (sym->ts.type == BT_UNKNOWN && sym->attr.intrinsic)
1142 gfc_intrinsic_sym *isym;
1144 isym = gfc_find_function (sym->name);
1145 if (isym == NULL || !isym->specific)
1147 gfc_error ("Unable to find a specific INTRINSIC procedure "
1148 "for the reference '%s' at %L", sym->name,
1153 sym->attr.intrinsic = 1;
1154 sym->attr.function = 1;
1159 /* See if the name is a module procedure in a parent unit. */
1161 if (was_declared (sym) || sym->ns->parent == NULL)
1164 if (gfc_find_sym_tree (sym->name, sym->ns->parent, 1, &parent_st))
1166 gfc_error ("Symbol '%s' at %L is ambiguous", sym->name, &e->where);
1170 if (parent_st == NULL)
1173 sym = parent_st->n.sym;
1174 e->symtree = parent_st; /* Point to the right thing. */
1176 if (sym->attr.flavor == FL_PROCEDURE
1177 || sym->attr.intrinsic
1178 || sym->attr.external)
1184 e->expr_type = EXPR_VARIABLE;
1186 if (sym->as != NULL)
1188 e->rank = sym->as->rank;
1189 e->ref = gfc_get_ref ();
1190 e->ref->type = REF_ARRAY;
1191 e->ref->u.ar.type = AR_FULL;
1192 e->ref->u.ar.as = sym->as;
1195 /* Expressions are assigned a default ts.type of BT_PROCEDURE in
1196 primary.c (match_actual_arg). If above code determines that it
1197 is a variable instead, it needs to be resolved as it was not
1198 done at the beginning of this function. */
1199 save_need_full_assumed_size = need_full_assumed_size;
1200 if (e->expr_type != FL_VARIABLE)
1201 need_full_assumed_size = 0;
1202 if (gfc_resolve_expr (e) != SUCCESS)
1204 need_full_assumed_size = save_need_full_assumed_size;
1207 /* Check argument list functions %VAL, %LOC and %REF. There is
1208 nothing to do for %REF. */
1209 if (arg->name && arg->name[0] == '%')
1211 if (strncmp ("%VAL", arg->name, 4) == 0)
1213 if (e->ts.type == BT_CHARACTER || e->ts.type == BT_DERIVED)
1215 gfc_error ("By-value argument at %L is not of numeric "
1222 gfc_error ("By-value argument at %L cannot be an array or "
1223 "an array section", &e->where);
1227 /* Intrinsics are still PROC_UNKNOWN here. However,
1228 since same file external procedures are not resolvable
1229 in gfortran, it is a good deal easier to leave them to
1231 if (ptype != PROC_UNKNOWN
1232 && ptype != PROC_DUMMY
1233 && ptype != PROC_EXTERNAL
1234 && ptype != PROC_MODULE)
1236 gfc_error ("By-value argument at %L is not allowed "
1237 "in this context", &e->where);
1242 /* Statement functions have already been excluded above. */
1243 else if (strncmp ("%LOC", arg->name, 4) == 0
1244 && e->ts.type == BT_PROCEDURE)
1246 if (e->symtree->n.sym->attr.proc == PROC_INTERNAL)
1248 gfc_error ("Passing internal procedure at %L by location "
1249 "not allowed", &e->where);
1260 /* Do the checks of the actual argument list that are specific to elemental
1261 procedures. If called with c == NULL, we have a function, otherwise if
1262 expr == NULL, we have a subroutine. */
1265 resolve_elemental_actual (gfc_expr *expr, gfc_code *c)
1267 gfc_actual_arglist *arg0;
1268 gfc_actual_arglist *arg;
1269 gfc_symbol *esym = NULL;
1270 gfc_intrinsic_sym *isym = NULL;
1272 gfc_intrinsic_arg *iformal = NULL;
1273 gfc_formal_arglist *eformal = NULL;
1274 bool formal_optional = false;
1275 bool set_by_optional = false;
1279 /* Is this an elemental procedure? */
1280 if (expr && expr->value.function.actual != NULL)
1282 if (expr->value.function.esym != NULL
1283 && expr->value.function.esym->attr.elemental)
1285 arg0 = expr->value.function.actual;
1286 esym = expr->value.function.esym;
1288 else if (expr->value.function.isym != NULL
1289 && expr->value.function.isym->elemental)
1291 arg0 = expr->value.function.actual;
1292 isym = expr->value.function.isym;
1297 else if (c && c->ext.actual != NULL && c->symtree->n.sym->attr.elemental)
1299 arg0 = c->ext.actual;
1300 esym = c->symtree->n.sym;
1305 /* The rank of an elemental is the rank of its array argument(s). */
1306 for (arg = arg0; arg; arg = arg->next)
1308 if (arg->expr != NULL && arg->expr->rank > 0)
1310 rank = arg->expr->rank;
1311 if (arg->expr->expr_type == EXPR_VARIABLE
1312 && arg->expr->symtree->n.sym->attr.optional)
1313 set_by_optional = true;
1315 /* Function specific; set the result rank and shape. */
1319 if (!expr->shape && arg->expr->shape)
1321 expr->shape = gfc_get_shape (rank);
1322 for (i = 0; i < rank; i++)
1323 mpz_init_set (expr->shape[i], arg->expr->shape[i]);
1330 /* If it is an array, it shall not be supplied as an actual argument
1331 to an elemental procedure unless an array of the same rank is supplied
1332 as an actual argument corresponding to a nonoptional dummy argument of
1333 that elemental procedure(12.4.1.5). */
1334 formal_optional = false;
1336 iformal = isym->formal;
1338 eformal = esym->formal;
1340 for (arg = arg0; arg; arg = arg->next)
1344 if (eformal->sym && eformal->sym->attr.optional)
1345 formal_optional = true;
1346 eformal = eformal->next;
1348 else if (isym && iformal)
1350 if (iformal->optional)
1351 formal_optional = true;
1352 iformal = iformal->next;
1355 formal_optional = true;
1357 if (pedantic && arg->expr != NULL
1358 && arg->expr->expr_type == EXPR_VARIABLE
1359 && arg->expr->symtree->n.sym->attr.optional
1362 && (set_by_optional || arg->expr->rank != rank)
1363 && !(isym && isym->id == GFC_ISYM_CONVERSION))
1365 gfc_warning ("'%s' at %L is an array and OPTIONAL; IF IT IS "
1366 "MISSING, it cannot be the actual argument of an "
1367 "ELEMENTAL procedure unless there is a non-optional "
1368 "argument with the same rank (12.4.1.5)",
1369 arg->expr->symtree->n.sym->name, &arg->expr->where);
1374 for (arg = arg0; arg; arg = arg->next)
1376 if (arg->expr == NULL || arg->expr->rank == 0)
1379 /* Being elemental, the last upper bound of an assumed size array
1380 argument must be present. */
1381 if (resolve_assumed_size_actual (arg->expr))
1384 /* Elemental procedure's array actual arguments must conform. */
1387 if (gfc_check_conformance ("elemental procedure", arg->expr, e)
1395 /* INTENT(OUT) is only allowed for subroutines; if any actual argument
1396 is an array, the intent inout/out variable needs to be also an array. */
1397 if (rank > 0 && esym && expr == NULL)
1398 for (eformal = esym->formal, arg = arg0; arg && eformal;
1399 arg = arg->next, eformal = eformal->next)
1400 if ((eformal->sym->attr.intent == INTENT_OUT
1401 || eformal->sym->attr.intent == INTENT_INOUT)
1402 && arg->expr && arg->expr->rank == 0)
1404 gfc_error ("Actual argument at %L for INTENT(%s) dummy '%s' of "
1405 "ELEMENTAL subroutine '%s' is a scalar, but another "
1406 "actual argument is an array", &arg->expr->where,
1407 (eformal->sym->attr.intent == INTENT_OUT) ? "OUT"
1408 : "INOUT", eformal->sym->name, esym->name);
1415 /* Go through each actual argument in ACTUAL and see if it can be
1416 implemented as an inlined, non-copying intrinsic. FNSYM is the
1417 function being called, or NULL if not known. */
1420 find_noncopying_intrinsics (gfc_symbol *fnsym, gfc_actual_arglist *actual)
1422 gfc_actual_arglist *ap;
1425 for (ap = actual; ap; ap = ap->next)
1427 && (expr = gfc_get_noncopying_intrinsic_argument (ap->expr))
1428 && !gfc_check_fncall_dependency (expr, INTENT_IN, fnsym, actual))
1429 ap->expr->inline_noncopying_intrinsic = 1;
1433 /* This function does the checking of references to global procedures
1434 as defined in sections 18.1 and 14.1, respectively, of the Fortran
1435 77 and 95 standards. It checks for a gsymbol for the name, making
1436 one if it does not already exist. If it already exists, then the
1437 reference being resolved must correspond to the type of gsymbol.
1438 Otherwise, the new symbol is equipped with the attributes of the
1439 reference. The corresponding code that is called in creating
1440 global entities is parse.c. */
1443 resolve_global_procedure (gfc_symbol *sym, locus *where, int sub)
1448 type = sub ? GSYM_SUBROUTINE : GSYM_FUNCTION;
1450 gsym = gfc_get_gsymbol (sym->name);
1452 if ((gsym->type != GSYM_UNKNOWN && gsym->type != type))
1453 gfc_global_used (gsym, where);
1455 if (gsym->type == GSYM_UNKNOWN)
1458 gsym->where = *where;
1465 /************* Function resolution *************/
1467 /* Resolve a function call known to be generic.
1468 Section 14.1.2.4.1. */
1471 resolve_generic_f0 (gfc_expr *expr, gfc_symbol *sym)
1475 if (sym->attr.generic)
1477 s = gfc_search_interface (sym->generic, 0, &expr->value.function.actual);
1480 expr->value.function.name = s->name;
1481 expr->value.function.esym = s;
1483 if (s->ts.type != BT_UNKNOWN)
1485 else if (s->result != NULL && s->result->ts.type != BT_UNKNOWN)
1486 expr->ts = s->result->ts;
1489 expr->rank = s->as->rank;
1490 else if (s->result != NULL && s->result->as != NULL)
1491 expr->rank = s->result->as->rank;
1493 gfc_set_sym_referenced (expr->value.function.esym);
1498 /* TODO: Need to search for elemental references in generic
1502 if (sym->attr.intrinsic)
1503 return gfc_intrinsic_func_interface (expr, 0);
1510 resolve_generic_f (gfc_expr *expr)
1515 sym = expr->symtree->n.sym;
1519 m = resolve_generic_f0 (expr, sym);
1522 else if (m == MATCH_ERROR)
1526 if (sym->ns->parent == NULL)
1528 gfc_find_symbol (sym->name, sym->ns->parent, 1, &sym);
1532 if (!generic_sym (sym))
1536 /* Last ditch attempt. See if the reference is to an intrinsic
1537 that possesses a matching interface. 14.1.2.4 */
1538 if (sym && !gfc_intrinsic_name (sym->name, 0))
1540 gfc_error ("There is no specific function for the generic '%s' at %L",
1541 expr->symtree->n.sym->name, &expr->where);
1545 m = gfc_intrinsic_func_interface (expr, 0);
1549 gfc_error ("Generic function '%s' at %L is not consistent with a "
1550 "specific intrinsic interface", expr->symtree->n.sym->name,
1557 /* Resolve a function call known to be specific. */
1560 resolve_specific_f0 (gfc_symbol *sym, gfc_expr *expr)
1564 /* See if we have an intrinsic interface. */
1566 if (sym->ts.interface != NULL && sym->ts.interface->attr.intrinsic)
1568 gfc_intrinsic_sym *isym;
1569 isym = gfc_find_function (sym->ts.interface->name);
1571 /* Existance of isym should be checked already. */
1575 sym->attr.function = 1;
1576 sym->attr.proc = PROC_EXTERNAL;
1580 if (sym->attr.external || sym->attr.if_source == IFSRC_IFBODY)
1582 if (sym->attr.dummy)
1584 sym->attr.proc = PROC_DUMMY;
1588 sym->attr.proc = PROC_EXTERNAL;
1592 if (sym->attr.proc == PROC_MODULE
1593 || sym->attr.proc == PROC_ST_FUNCTION
1594 || sym->attr.proc == PROC_INTERNAL)
1597 if (sym->attr.intrinsic)
1599 m = gfc_intrinsic_func_interface (expr, 1);
1603 gfc_error ("Function '%s' at %L is INTRINSIC but is not compatible "
1604 "with an intrinsic", sym->name, &expr->where);
1612 gfc_procedure_use (sym, &expr->value.function.actual, &expr->where);
1615 expr->value.function.name = sym->name;
1616 expr->value.function.esym = sym;
1617 if (sym->as != NULL)
1618 expr->rank = sym->as->rank;
1625 resolve_specific_f (gfc_expr *expr)
1630 sym = expr->symtree->n.sym;
1634 m = resolve_specific_f0 (sym, expr);
1637 if (m == MATCH_ERROR)
1640 if (sym->ns->parent == NULL)
1643 gfc_find_symbol (sym->name, sym->ns->parent, 1, &sym);
1649 gfc_error ("Unable to resolve the specific function '%s' at %L",
1650 expr->symtree->n.sym->name, &expr->where);
1656 /* Resolve a procedure call not known to be generic nor specific. */
1659 resolve_unknown_f (gfc_expr *expr)
1664 sym = expr->symtree->n.sym;
1666 if (sym->attr.dummy)
1668 sym->attr.proc = PROC_DUMMY;
1669 expr->value.function.name = sym->name;
1673 /* See if we have an intrinsic function reference. */
1675 if (gfc_intrinsic_name (sym->name, 0))
1677 if (gfc_intrinsic_func_interface (expr, 1) == MATCH_YES)
1682 /* The reference is to an external name. */
1684 sym->attr.proc = PROC_EXTERNAL;
1685 expr->value.function.name = sym->name;
1686 expr->value.function.esym = expr->symtree->n.sym;
1688 if (sym->as != NULL)
1689 expr->rank = sym->as->rank;
1691 /* Type of the expression is either the type of the symbol or the
1692 default type of the symbol. */
1695 gfc_procedure_use (sym, &expr->value.function.actual, &expr->where);
1697 if (sym->ts.type != BT_UNKNOWN)
1701 ts = gfc_get_default_type (sym, sym->ns);
1703 if (ts->type == BT_UNKNOWN)
1705 gfc_error ("Function '%s' at %L has no IMPLICIT type",
1706 sym->name, &expr->where);
1717 /* Return true, if the symbol is an external procedure. */
1719 is_external_proc (gfc_symbol *sym)
1721 if (!sym->attr.dummy && !sym->attr.contained
1722 && !(sym->attr.intrinsic
1723 || gfc_intrinsic_name (sym->name, sym->attr.subroutine))
1724 && sym->attr.proc != PROC_ST_FUNCTION
1725 && !sym->attr.use_assoc
1733 /* Figure out if a function reference is pure or not. Also set the name
1734 of the function for a potential error message. Return nonzero if the
1735 function is PURE, zero if not. */
1737 pure_stmt_function (gfc_expr *, gfc_symbol *);
1740 pure_function (gfc_expr *e, const char **name)
1746 if (e->symtree != NULL
1747 && e->symtree->n.sym != NULL
1748 && e->symtree->n.sym->attr.proc == PROC_ST_FUNCTION)
1749 return pure_stmt_function (e, e->symtree->n.sym);
1751 if (e->value.function.esym)
1753 pure = gfc_pure (e->value.function.esym);
1754 *name = e->value.function.esym->name;
1756 else if (e->value.function.isym)
1758 pure = e->value.function.isym->pure
1759 || e->value.function.isym->elemental;
1760 *name = e->value.function.isym->name;
1764 /* Implicit functions are not pure. */
1766 *name = e->value.function.name;
1774 impure_stmt_fcn (gfc_expr *e, gfc_symbol *sym,
1775 int *f ATTRIBUTE_UNUSED)
1779 /* Don't bother recursing into other statement functions
1780 since they will be checked individually for purity. */
1781 if (e->expr_type != EXPR_FUNCTION
1783 || e->symtree->n.sym == sym
1784 || e->symtree->n.sym->attr.proc == PROC_ST_FUNCTION)
1787 return pure_function (e, &name) ? false : true;
1792 pure_stmt_function (gfc_expr *e, gfc_symbol *sym)
1794 return gfc_traverse_expr (e, sym, impure_stmt_fcn, 0) ? 0 : 1;
1799 is_scalar_expr_ptr (gfc_expr *expr)
1801 try retval = SUCCESS;
1806 /* See if we have a gfc_ref, which means we have a substring, array
1807 reference, or a component. */
1808 if (expr->ref != NULL)
1811 while (ref->next != NULL)
1817 if (ref->u.ss.length != NULL
1818 && ref->u.ss.length->length != NULL
1820 && ref->u.ss.start->expr_type == EXPR_CONSTANT
1822 && ref->u.ss.end->expr_type == EXPR_CONSTANT)
1824 start = (int) mpz_get_si (ref->u.ss.start->value.integer);
1825 end = (int) mpz_get_si (ref->u.ss.end->value.integer);
1826 if (end - start + 1 != 1)
1833 if (ref->u.ar.type == AR_ELEMENT)
1835 else if (ref->u.ar.type == AR_FULL)
1837 /* The user can give a full array if the array is of size 1. */
1838 if (ref->u.ar.as != NULL
1839 && ref->u.ar.as->rank == 1
1840 && ref->u.ar.as->type == AS_EXPLICIT
1841 && ref->u.ar.as->lower[0] != NULL
1842 && ref->u.ar.as->lower[0]->expr_type == EXPR_CONSTANT
1843 && ref->u.ar.as->upper[0] != NULL
1844 && ref->u.ar.as->upper[0]->expr_type == EXPR_CONSTANT)
1846 /* If we have a character string, we need to check if
1847 its length is one. */
1848 if (expr->ts.type == BT_CHARACTER)
1850 if (expr->ts.cl == NULL
1851 || expr->ts.cl->length == NULL
1852 || mpz_cmp_si (expr->ts.cl->length->value.integer, 1)
1858 /* We have constant lower and upper bounds. If the
1859 difference between is 1, it can be considered a
1861 start = (int) mpz_get_si
1862 (ref->u.ar.as->lower[0]->value.integer);
1863 end = (int) mpz_get_si
1864 (ref->u.ar.as->upper[0]->value.integer);
1865 if (end - start + 1 != 1)
1880 else if (expr->ts.type == BT_CHARACTER && expr->rank == 0)
1882 /* Character string. Make sure it's of length 1. */
1883 if (expr->ts.cl == NULL
1884 || expr->ts.cl->length == NULL
1885 || mpz_cmp_si (expr->ts.cl->length->value.integer, 1) != 0)
1888 else if (expr->rank != 0)
1895 /* Match one of the iso_c_binding functions (c_associated or c_loc)
1896 and, in the case of c_associated, set the binding label based on
1900 gfc_iso_c_func_interface (gfc_symbol *sym, gfc_actual_arglist *args,
1901 gfc_symbol **new_sym)
1903 char name[GFC_MAX_SYMBOL_LEN + 1];
1904 char binding_label[GFC_MAX_BINDING_LABEL_LEN + 1];
1905 int optional_arg = 0;
1906 try retval = SUCCESS;
1907 gfc_symbol *args_sym;
1908 gfc_typespec *arg_ts;
1909 gfc_ref *parent_ref;
1912 if (args->expr->expr_type == EXPR_CONSTANT
1913 || args->expr->expr_type == EXPR_OP
1914 || args->expr->expr_type == EXPR_NULL)
1916 gfc_error ("Argument to '%s' at %L is not a variable",
1917 sym->name, &(args->expr->where));
1921 args_sym = args->expr->symtree->n.sym;
1923 /* The typespec for the actual arg should be that stored in the expr
1924 and not necessarily that of the expr symbol (args_sym), because
1925 the actual expression could be a part-ref of the expr symbol. */
1926 arg_ts = &(args->expr->ts);
1928 /* Get the parent reference (if any) for the expression. This happens for
1929 cases such as a%b%c. */
1930 parent_ref = args->expr->ref;
1932 if (parent_ref != NULL)
1934 curr_ref = parent_ref->next;
1935 while (curr_ref != NULL && curr_ref->next != NULL)
1937 parent_ref = curr_ref;
1938 curr_ref = curr_ref->next;
1942 /* If curr_ref is non-NULL, we had a part-ref expression. If the curr_ref
1943 is for a REF_COMPONENT, then we need to use it as the parent_ref for
1944 the name, etc. Otherwise, the current parent_ref should be correct. */
1945 if (curr_ref != NULL && curr_ref->type == REF_COMPONENT)
1946 parent_ref = curr_ref;
1948 if (parent_ref == args->expr->ref)
1950 else if (parent_ref != NULL && parent_ref->type != REF_COMPONENT)
1951 gfc_internal_error ("Unexpected expression reference type in "
1952 "gfc_iso_c_func_interface");
1954 if (sym->intmod_sym_id == ISOCBINDING_ASSOCIATED)
1956 /* If the user gave two args then they are providing something for
1957 the optional arg (the second cptr). Therefore, set the name and
1958 binding label to the c_associated for two cptrs. Otherwise,
1959 set c_associated to expect one cptr. */
1963 sprintf (name, "%s_2", sym->name);
1964 sprintf (binding_label, "%s_2", sym->binding_label);
1970 sprintf (name, "%s_1", sym->name);
1971 sprintf (binding_label, "%s_1", sym->binding_label);
1975 /* Get a new symbol for the version of c_associated that
1977 *new_sym = get_iso_c_sym (sym, name, binding_label, optional_arg);
1979 else if (sym->intmod_sym_id == ISOCBINDING_LOC
1980 || sym->intmod_sym_id == ISOCBINDING_FUNLOC)
1982 sprintf (name, "%s", sym->name);
1983 sprintf (binding_label, "%s", sym->binding_label);
1985 /* Error check the call. */
1986 if (args->next != NULL)
1988 gfc_error_now ("More actual than formal arguments in '%s' "
1989 "call at %L", name, &(args->expr->where));
1992 else if (sym->intmod_sym_id == ISOCBINDING_LOC)
1994 /* Make sure we have either the target or pointer attribute. */
1995 if (!(args_sym->attr.target)
1996 && !(args_sym->attr.pointer)
1997 && (parent_ref == NULL ||
1998 !parent_ref->u.c.component->pointer))
2000 gfc_error_now ("Parameter '%s' to '%s' at %L must be either "
2001 "a TARGET or an associated pointer",
2003 sym->name, &(args->expr->where));
2007 /* See if we have interoperable type and type param. */
2008 if (verify_c_interop (arg_ts,
2009 (parent_ref ? parent_ref->u.c.component->name
2011 &(args->expr->where)) == SUCCESS
2012 || gfc_check_any_c_kind (arg_ts) == SUCCESS)
2014 if (args_sym->attr.target == 1)
2016 /* Case 1a, section 15.1.2.5, J3/04-007: variable that
2017 has the target attribute and is interoperable. */
2018 /* Case 1b, section 15.1.2.5, J3/04-007: allocated
2019 allocatable variable that has the TARGET attribute and
2020 is not an array of zero size. */
2021 if (args_sym->attr.allocatable == 1)
2023 if (args_sym->attr.dimension != 0
2024 && (args_sym->as && args_sym->as->rank == 0))
2026 gfc_error_now ("Allocatable variable '%s' used as a "
2027 "parameter to '%s' at %L must not be "
2028 "an array of zero size",
2029 args_sym->name, sym->name,
2030 &(args->expr->where));
2036 /* A non-allocatable target variable with C
2037 interoperable type and type parameters must be
2039 if (args_sym && args_sym->attr.dimension)
2041 if (args_sym->as->type == AS_ASSUMED_SHAPE)
2043 gfc_error ("Assumed-shape array '%s' at %L "
2044 "cannot be an argument to the "
2045 "procedure '%s' because "
2046 "it is not C interoperable",
2048 &(args->expr->where), sym->name);
2051 else if (args_sym->as->type == AS_DEFERRED)
2053 gfc_error ("Deferred-shape array '%s' at %L "
2054 "cannot be an argument to the "
2055 "procedure '%s' because "
2056 "it is not C interoperable",
2058 &(args->expr->where), sym->name);
2063 /* Make sure it's not a character string. Arrays of
2064 any type should be ok if the variable is of a C
2065 interoperable type. */
2066 if (arg_ts->type == BT_CHARACTER)
2067 if (arg_ts->cl != NULL
2068 && (arg_ts->cl->length == NULL
2069 || arg_ts->cl->length->expr_type
2072 (arg_ts->cl->length->value.integer, 1)
2074 && is_scalar_expr_ptr (args->expr) != SUCCESS)
2076 gfc_error_now ("CHARACTER argument '%s' to '%s' "
2077 "at %L must have a length of 1",
2078 args_sym->name, sym->name,
2079 &(args->expr->where));
2084 else if ((args_sym->attr.pointer == 1 ||
2086 && parent_ref->u.c.component->pointer))
2087 && is_scalar_expr_ptr (args->expr) != SUCCESS)
2089 /* Case 1c, section 15.1.2.5, J3/04-007: an associated
2091 gfc_error_now ("Argument '%s' to '%s' at %L must be an "
2092 "associated scalar POINTER", args_sym->name,
2093 sym->name, &(args->expr->where));
2099 /* The parameter is not required to be C interoperable. If it
2100 is not C interoperable, it must be a nonpolymorphic scalar
2101 with no length type parameters. It still must have either
2102 the pointer or target attribute, and it can be
2103 allocatable (but must be allocated when c_loc is called). */
2104 if (args->expr->rank != 0
2105 && is_scalar_expr_ptr (args->expr) != SUCCESS)
2107 gfc_error_now ("Parameter '%s' to '%s' at %L must be a "
2108 "scalar", args_sym->name, sym->name,
2109 &(args->expr->where));
2112 else if (arg_ts->type == BT_CHARACTER
2113 && is_scalar_expr_ptr (args->expr) != SUCCESS)
2115 gfc_error_now ("CHARACTER argument '%s' to '%s' at "
2116 "%L must have a length of 1",
2117 args_sym->name, sym->name,
2118 &(args->expr->where));
2123 else if (sym->intmod_sym_id == ISOCBINDING_FUNLOC)
2125 if (args_sym->attr.flavor != FL_PROCEDURE)
2127 /* TODO: Update this error message to allow for procedure
2128 pointers once they are implemented. */
2129 gfc_error_now ("Parameter '%s' to '%s' at %L must be a "
2131 args_sym->name, sym->name,
2132 &(args->expr->where));
2135 else if (args_sym->attr.is_bind_c != 1)
2137 gfc_error_now ("Parameter '%s' to '%s' at %L must be "
2139 args_sym->name, sym->name,
2140 &(args->expr->where));
2145 /* for c_loc/c_funloc, the new symbol is the same as the old one */
2150 gfc_internal_error ("gfc_iso_c_func_interface(): Unhandled "
2151 "iso_c_binding function: '%s'!\n", sym->name);
2158 /* Resolve a function call, which means resolving the arguments, then figuring
2159 out which entity the name refers to. */
2160 /* TODO: Check procedure arguments so that an INTENT(IN) isn't passed
2161 to INTENT(OUT) or INTENT(INOUT). */
2164 resolve_function (gfc_expr *expr)
2166 gfc_actual_arglist *arg;
2171 procedure_type p = PROC_INTRINSIC;
2175 sym = expr->symtree->n.sym;
2177 if (sym && sym->attr.flavor == FL_VARIABLE)
2179 gfc_error ("'%s' at %L is not a function", sym->name, &expr->where);
2183 if (sym && sym->attr.abstract)
2185 gfc_error ("ABSTRACT INTERFACE '%s' must not be referenced at %L",
2186 sym->name, &expr->where);
2190 /* If the procedure is external, check for usage. */
2191 if (sym && is_external_proc (sym))
2192 resolve_global_procedure (sym, &expr->where, 0);
2194 /* Switch off assumed size checking and do this again for certain kinds
2195 of procedure, once the procedure itself is resolved. */
2196 need_full_assumed_size++;
2198 if (expr->symtree && expr->symtree->n.sym)
2199 p = expr->symtree->n.sym->attr.proc;
2201 if (resolve_actual_arglist (expr->value.function.actual, p) == FAILURE)
2204 /* Need to setup the call to the correct c_associated, depending on
2205 the number of cptrs to user gives to compare. */
2206 if (sym && sym->attr.is_iso_c == 1)
2208 if (gfc_iso_c_func_interface (sym, expr->value.function.actual, &sym)
2212 /* Get the symtree for the new symbol (resolved func).
2213 the old one will be freed later, when it's no longer used. */
2214 gfc_find_sym_tree (sym->name, sym->ns, 1, &(expr->symtree));
2217 /* Resume assumed_size checking. */
2218 need_full_assumed_size--;
2220 if (sym && sym->ts.type == BT_CHARACTER
2222 && sym->ts.cl->length == NULL
2224 && expr->value.function.esym == NULL
2225 && !sym->attr.contained)
2227 /* Internal procedures are taken care of in resolve_contained_fntype. */
2228 gfc_error ("Function '%s' is declared CHARACTER(*) and cannot "
2229 "be used at %L since it is not a dummy argument",
2230 sym->name, &expr->where);
2234 /* See if function is already resolved. */
2236 if (expr->value.function.name != NULL)
2238 if (expr->ts.type == BT_UNKNOWN)
2244 /* Apply the rules of section 14.1.2. */
2246 switch (procedure_kind (sym))
2249 t = resolve_generic_f (expr);
2252 case PTYPE_SPECIFIC:
2253 t = resolve_specific_f (expr);
2257 t = resolve_unknown_f (expr);
2261 gfc_internal_error ("resolve_function(): bad function type");
2265 /* If the expression is still a function (it might have simplified),
2266 then we check to see if we are calling an elemental function. */
2268 if (expr->expr_type != EXPR_FUNCTION)
2271 temp = need_full_assumed_size;
2272 need_full_assumed_size = 0;
2274 if (resolve_elemental_actual (expr, NULL) == FAILURE)
2277 if (omp_workshare_flag
2278 && expr->value.function.esym
2279 && ! gfc_elemental (expr->value.function.esym))
2281 gfc_error ("User defined non-ELEMENTAL function '%s' at %L not allowed "
2282 "in WORKSHARE construct", expr->value.function.esym->name,
2287 #define GENERIC_ID expr->value.function.isym->id
2288 else if (expr->value.function.actual != NULL
2289 && expr->value.function.isym != NULL
2290 && GENERIC_ID != GFC_ISYM_LBOUND
2291 && GENERIC_ID != GFC_ISYM_LEN
2292 && GENERIC_ID != GFC_ISYM_LOC
2293 && GENERIC_ID != GFC_ISYM_PRESENT)
2295 /* Array intrinsics must also have the last upper bound of an
2296 assumed size array argument. UBOUND and SIZE have to be
2297 excluded from the check if the second argument is anything
2300 inquiry = GENERIC_ID == GFC_ISYM_UBOUND
2301 || GENERIC_ID == GFC_ISYM_SIZE;
2303 for (arg = expr->value.function.actual; arg; arg = arg->next)
2305 if (inquiry && arg->next != NULL && arg->next->expr)
2307 if (arg->next->expr->expr_type != EXPR_CONSTANT)
2310 if ((int)mpz_get_si (arg->next->expr->value.integer)
2315 if (arg->expr != NULL
2316 && arg->expr->rank > 0
2317 && resolve_assumed_size_actual (arg->expr))
2323 need_full_assumed_size = temp;
2326 if (!pure_function (expr, &name) && name)
2330 gfc_error ("reference to non-PURE function '%s' at %L inside a "
2331 "FORALL %s", name, &expr->where,
2332 forall_flag == 2 ? "mask" : "block");
2335 else if (gfc_pure (NULL))
2337 gfc_error ("Function reference to '%s' at %L is to a non-PURE "
2338 "procedure within a PURE procedure", name, &expr->where);
2343 /* Functions without the RECURSIVE attribution are not allowed to
2344 * call themselves. */
2345 if (expr->value.function.esym && !expr->value.function.esym->attr.recursive)
2347 gfc_symbol *esym, *proc;
2348 esym = expr->value.function.esym;
2349 proc = gfc_current_ns->proc_name;
2352 gfc_error ("Function '%s' at %L cannot call itself, as it is not "
2353 "RECURSIVE", name, &expr->where);
2357 if (esym->attr.entry && esym->ns->entries && proc->ns->entries
2358 && esym->ns->entries->sym == proc->ns->entries->sym)
2360 gfc_error ("Call to ENTRY '%s' at %L is recursive, but function "
2361 "'%s' is not declared as RECURSIVE",
2362 esym->name, &expr->where, esym->ns->entries->sym->name);
2367 /* Character lengths of use associated functions may contains references to
2368 symbols not referenced from the current program unit otherwise. Make sure
2369 those symbols are marked as referenced. */
2371 if (expr->ts.type == BT_CHARACTER && expr->value.function.esym
2372 && expr->value.function.esym->attr.use_assoc)
2374 gfc_expr_set_symbols_referenced (expr->ts.cl->length);
2378 && !((expr->value.function.esym
2379 && expr->value.function.esym->attr.elemental)
2381 (expr->value.function.isym
2382 && expr->value.function.isym->elemental)))
2383 find_noncopying_intrinsics (expr->value.function.esym,
2384 expr->value.function.actual);
2386 /* Make sure that the expression has a typespec that works. */
2387 if (expr->ts.type == BT_UNKNOWN)
2389 if (expr->symtree->n.sym->result
2390 && expr->symtree->n.sym->result->ts.type != BT_UNKNOWN)
2391 expr->ts = expr->symtree->n.sym->result->ts;
2398 /************* Subroutine resolution *************/
2401 pure_subroutine (gfc_code *c, gfc_symbol *sym)
2407 gfc_error ("Subroutine call to '%s' in FORALL block at %L is not PURE",
2408 sym->name, &c->loc);
2409 else if (gfc_pure (NULL))
2410 gfc_error ("Subroutine call to '%s' at %L is not PURE", sym->name,
2416 resolve_generic_s0 (gfc_code *c, gfc_symbol *sym)
2420 if (sym->attr.generic)
2422 s = gfc_search_interface (sym->generic, 1, &c->ext.actual);
2425 c->resolved_sym = s;
2426 pure_subroutine (c, s);
2430 /* TODO: Need to search for elemental references in generic interface. */
2433 if (sym->attr.intrinsic)
2434 return gfc_intrinsic_sub_interface (c, 0);
2441 resolve_generic_s (gfc_code *c)
2446 sym = c->symtree->n.sym;
2450 m = resolve_generic_s0 (c, sym);
2453 else if (m == MATCH_ERROR)
2457 if (sym->ns->parent == NULL)
2459 gfc_find_symbol (sym->name, sym->ns->parent, 1, &sym);
2463 if (!generic_sym (sym))
2467 /* Last ditch attempt. See if the reference is to an intrinsic
2468 that possesses a matching interface. 14.1.2.4 */
2469 sym = c->symtree->n.sym;
2471 if (!gfc_intrinsic_name (sym->name, 1))
2473 gfc_error ("There is no specific subroutine for the generic '%s' at %L",
2474 sym->name, &c->loc);
2478 m = gfc_intrinsic_sub_interface (c, 0);
2482 gfc_error ("Generic subroutine '%s' at %L is not consistent with an "
2483 "intrinsic subroutine interface", sym->name, &c->loc);
2489 /* Set the name and binding label of the subroutine symbol in the call
2490 expression represented by 'c' to include the type and kind of the
2491 second parameter. This function is for resolving the appropriate
2492 version of c_f_pointer() and c_f_procpointer(). For example, a
2493 call to c_f_pointer() for a default integer pointer could have a
2494 name of c_f_pointer_i4. If no second arg exists, which is an error
2495 for these two functions, it defaults to the generic symbol's name
2496 and binding label. */
2499 set_name_and_label (gfc_code *c, gfc_symbol *sym,
2500 char *name, char *binding_label)
2502 gfc_expr *arg = NULL;
2506 /* The second arg of c_f_pointer and c_f_procpointer determines
2507 the type and kind for the procedure name. */
2508 arg = c->ext.actual->next->expr;
2512 /* Set up the name to have the given symbol's name,
2513 plus the type and kind. */
2514 /* a derived type is marked with the type letter 'u' */
2515 if (arg->ts.type == BT_DERIVED)
2518 kind = 0; /* set the kind as 0 for now */
2522 type = gfc_type_letter (arg->ts.type);
2523 kind = arg->ts.kind;
2526 if (arg->ts.type == BT_CHARACTER)
2527 /* Kind info for character strings not needed. */
2530 sprintf (name, "%s_%c%d", sym->name, type, kind);
2531 /* Set up the binding label as the given symbol's label plus
2532 the type and kind. */
2533 sprintf (binding_label, "%s_%c%d", sym->binding_label, type, kind);
2537 /* If the second arg is missing, set the name and label as
2538 was, cause it should at least be found, and the missing
2539 arg error will be caught by compare_parameters(). */
2540 sprintf (name, "%s", sym->name);
2541 sprintf (binding_label, "%s", sym->binding_label);
2548 /* Resolve a generic version of the iso_c_binding procedure given
2549 (sym) to the specific one based on the type and kind of the
2550 argument(s). Currently, this function resolves c_f_pointer() and
2551 c_f_procpointer based on the type and kind of the second argument
2552 (FPTR). Other iso_c_binding procedures aren't specially handled.
2553 Upon successfully exiting, c->resolved_sym will hold the resolved
2554 symbol. Returns MATCH_ERROR if an error occurred; MATCH_YES
2558 gfc_iso_c_sub_interface (gfc_code *c, gfc_symbol *sym)
2560 gfc_symbol *new_sym;
2561 /* this is fine, since we know the names won't use the max */
2562 char name[GFC_MAX_SYMBOL_LEN + 1];
2563 char binding_label[GFC_MAX_BINDING_LABEL_LEN + 1];
2564 /* default to success; will override if find error */
2565 match m = MATCH_YES;
2567 /* Make sure the actual arguments are in the necessary order (based on the
2568 formal args) before resolving. */
2569 gfc_procedure_use (sym, &c->ext.actual, &(c->loc));
2571 if ((sym->intmod_sym_id == ISOCBINDING_F_POINTER) ||
2572 (sym->intmod_sym_id == ISOCBINDING_F_PROCPOINTER))
2574 set_name_and_label (c, sym, name, binding_label);
2576 if (sym->intmod_sym_id == ISOCBINDING_F_POINTER)
2578 if (c->ext.actual != NULL && c->ext.actual->next != NULL)
2580 /* Make sure we got a third arg if the second arg has non-zero
2581 rank. We must also check that the type and rank are
2582 correct since we short-circuit this check in
2583 gfc_procedure_use() (called above to sort actual args). */
2584 if (c->ext.actual->next->expr->rank != 0)
2586 if(c->ext.actual->next->next == NULL
2587 || c->ext.actual->next->next->expr == NULL)
2590 gfc_error ("Missing SHAPE parameter for call to %s "
2591 "at %L", sym->name, &(c->loc));
2593 else if (c->ext.actual->next->next->expr->ts.type
2595 || c->ext.actual->next->next->expr->rank != 1)
2598 gfc_error ("SHAPE parameter for call to %s at %L must "
2599 "be a rank 1 INTEGER array", sym->name,
2606 if (m != MATCH_ERROR)
2608 /* the 1 means to add the optional arg to formal list */
2609 new_sym = get_iso_c_sym (sym, name, binding_label, 1);
2611 /* for error reporting, say it's declared where the original was */
2612 new_sym->declared_at = sym->declared_at;
2617 /* no differences for c_loc or c_funloc */
2621 /* set the resolved symbol */
2622 if (m != MATCH_ERROR)
2623 c->resolved_sym = new_sym;
2625 c->resolved_sym = sym;
2631 /* Resolve a subroutine call known to be specific. */
2634 resolve_specific_s0 (gfc_code *c, gfc_symbol *sym)
2638 /* See if we have an intrinsic interface. */
2639 if (sym->ts.interface != NULL && !sym->ts.interface->attr.abstract
2640 && !sym->ts.interface->attr.subroutine)
2642 gfc_intrinsic_sym *isym;
2644 isym = gfc_find_function (sym->ts.interface->name);
2646 /* Existance of isym should be checked already. */
2650 sym->attr.function = 1;
2654 if(sym->attr.is_iso_c)
2656 m = gfc_iso_c_sub_interface (c,sym);
2660 if (sym->attr.external || sym->attr.if_source == IFSRC_IFBODY)
2662 if (sym->attr.dummy)
2664 sym->attr.proc = PROC_DUMMY;
2668 sym->attr.proc = PROC_EXTERNAL;
2672 if (sym->attr.proc == PROC_MODULE || sym->attr.proc == PROC_INTERNAL)
2675 if (sym->attr.intrinsic)
2677 m = gfc_intrinsic_sub_interface (c, 1);
2681 gfc_error ("Subroutine '%s' at %L is INTRINSIC but is not compatible "
2682 "with an intrinsic", sym->name, &c->loc);
2690 gfc_procedure_use (sym, &c->ext.actual, &c->loc);
2692 c->resolved_sym = sym;
2693 pure_subroutine (c, sym);
2700 resolve_specific_s (gfc_code *c)
2705 sym = c->symtree->n.sym;
2709 m = resolve_specific_s0 (c, sym);
2712 if (m == MATCH_ERROR)
2715 if (sym->ns->parent == NULL)
2718 gfc_find_symbol (sym->name, sym->ns->parent, 1, &sym);
2724 sym = c->symtree->n.sym;
2725 gfc_error ("Unable to resolve the specific subroutine '%s' at %L",
2726 sym->name, &c->loc);
2732 /* Resolve a subroutine call not known to be generic nor specific. */
2735 resolve_unknown_s (gfc_code *c)
2739 sym = c->symtree->n.sym;
2741 if (sym->attr.dummy)
2743 sym->attr.proc = PROC_DUMMY;
2747 /* See if we have an intrinsic function reference. */
2749 if (gfc_intrinsic_name (sym->name, 1))
2751 if (gfc_intrinsic_sub_interface (c, 1) == MATCH_YES)
2756 /* The reference is to an external name. */
2759 gfc_procedure_use (sym, &c->ext.actual, &c->loc);
2761 c->resolved_sym = sym;
2763 pure_subroutine (c, sym);
2769 /* Resolve a subroutine call. Although it was tempting to use the same code
2770 for functions, subroutines and functions are stored differently and this
2771 makes things awkward. */
2774 resolve_call (gfc_code *c)
2777 procedure_type ptype = PROC_INTRINSIC;
2779 if (c->symtree && c->symtree->n.sym
2780 && c->symtree->n.sym->ts.type != BT_UNKNOWN)
2782 gfc_error ("'%s' at %L has a type, which is not consistent with "
2783 "the CALL at %L", c->symtree->n.sym->name,
2784 &c->symtree->n.sym->declared_at, &c->loc);
2788 /* If external, check for usage. */
2789 if (c->symtree && is_external_proc (c->symtree->n.sym))
2790 resolve_global_procedure (c->symtree->n.sym, &c->loc, 1);
2792 /* Subroutines without the RECURSIVE attribution are not allowed to
2793 * call themselves. */
2794 if (c->symtree && c->symtree->n.sym && !c->symtree->n.sym->attr.recursive)
2796 gfc_symbol *csym, *proc;
2797 csym = c->symtree->n.sym;
2798 proc = gfc_current_ns->proc_name;
2801 gfc_error ("SUBROUTINE '%s' at %L cannot call itself, as it is not "
2802 "RECURSIVE", csym->name, &c->loc);
2806 if (csym->attr.entry && csym->ns->entries && proc->ns->entries
2807 && csym->ns->entries->sym == proc->ns->entries->sym)
2809 gfc_error ("Call to ENTRY '%s' at %L is recursive, but subroutine "
2810 "'%s' is not declared as RECURSIVE",
2811 csym->name, &c->loc, csym->ns->entries->sym->name);
2816 /* Switch off assumed size checking and do this again for certain kinds
2817 of procedure, once the procedure itself is resolved. */
2818 need_full_assumed_size++;
2820 if (c->symtree && c->symtree->n.sym)
2821 ptype = c->symtree->n.sym->attr.proc;
2823 if (resolve_actual_arglist (c->ext.actual, ptype) == FAILURE)
2826 /* Resume assumed_size checking. */
2827 need_full_assumed_size--;
2830 if (c->resolved_sym == NULL)
2831 switch (procedure_kind (c->symtree->n.sym))
2834 t = resolve_generic_s (c);
2837 case PTYPE_SPECIFIC:
2838 t = resolve_specific_s (c);
2842 t = resolve_unknown_s (c);
2846 gfc_internal_error ("resolve_subroutine(): bad function type");
2849 /* Some checks of elemental subroutine actual arguments. */
2850 if (resolve_elemental_actual (NULL, c) == FAILURE)
2853 if (t == SUCCESS && !(c->resolved_sym && c->resolved_sym->attr.elemental))
2854 find_noncopying_intrinsics (c->resolved_sym, c->ext.actual);
2859 /* Compare the shapes of two arrays that have non-NULL shapes. If both
2860 op1->shape and op2->shape are non-NULL return SUCCESS if their shapes
2861 match. If both op1->shape and op2->shape are non-NULL return FAILURE
2862 if their shapes do not match. If either op1->shape or op2->shape is
2863 NULL, return SUCCESS. */
2866 compare_shapes (gfc_expr *op1, gfc_expr *op2)
2873 if (op1->shape != NULL && op2->shape != NULL)
2875 for (i = 0; i < op1->rank; i++)
2877 if (mpz_cmp (op1->shape[i], op2->shape[i]) != 0)
2879 gfc_error ("Shapes for operands at %L and %L are not conformable",
2880 &op1->where, &op2->where);
2891 /* Resolve an operator expression node. This can involve replacing the
2892 operation with a user defined function call. */
2895 resolve_operator (gfc_expr *e)
2897 gfc_expr *op1, *op2;
2899 bool dual_locus_error;
2902 /* Resolve all subnodes-- give them types. */
2904 switch (e->value.op.operator)
2907 if (gfc_resolve_expr (e->value.op.op2) == FAILURE)
2910 /* Fall through... */
2913 case INTRINSIC_UPLUS:
2914 case INTRINSIC_UMINUS:
2915 case INTRINSIC_PARENTHESES:
2916 if (gfc_resolve_expr (e->value.op.op1) == FAILURE)
2921 /* Typecheck the new node. */
2923 op1 = e->value.op.op1;
2924 op2 = e->value.op.op2;
2925 dual_locus_error = false;
2927 if ((op1 && op1->expr_type == EXPR_NULL)
2928 || (op2 && op2->expr_type == EXPR_NULL))
2930 sprintf (msg, _("Invalid context for NULL() pointer at %%L"));
2934 switch (e->value.op.operator)
2936 case INTRINSIC_UPLUS:
2937 case INTRINSIC_UMINUS:
2938 if (op1->ts.type == BT_INTEGER
2939 || op1->ts.type == BT_REAL
2940 || op1->ts.type == BT_COMPLEX)
2946 sprintf (msg, _("Operand of unary numeric operator '%s' at %%L is %s"),
2947 gfc_op2string (e->value.op.operator), gfc_typename (&e->ts));
2950 case INTRINSIC_PLUS:
2951 case INTRINSIC_MINUS:
2952 case INTRINSIC_TIMES:
2953 case INTRINSIC_DIVIDE:
2954 case INTRINSIC_POWER:
2955 if (gfc_numeric_ts (&op1->ts) && gfc_numeric_ts (&op2->ts))
2957 gfc_type_convert_binary (e);
2962 _("Operands of binary numeric operator '%s' at %%L are %s/%s"),
2963 gfc_op2string (e->value.op.operator), gfc_typename (&op1->ts),
2964 gfc_typename (&op2->ts));
2967 case INTRINSIC_CONCAT:
2968 if (op1->ts.type == BT_CHARACTER && op2->ts.type == BT_CHARACTER)
2970 e->ts.type = BT_CHARACTER;
2971 e->ts.kind = op1->ts.kind;
2976 _("Operands of string concatenation operator at %%L are %s/%s"),
2977 gfc_typename (&op1->ts), gfc_typename (&op2->ts));
2983 case INTRINSIC_NEQV:
2984 if (op1->ts.type == BT_LOGICAL && op2->ts.type == BT_LOGICAL)
2986 e->ts.type = BT_LOGICAL;
2987 e->ts.kind = gfc_kind_max (op1, op2);
2988 if (op1->ts.kind < e->ts.kind)
2989 gfc_convert_type (op1, &e->ts, 2);
2990 else if (op2->ts.kind < e->ts.kind)
2991 gfc_convert_type (op2, &e->ts, 2);
2995 sprintf (msg, _("Operands of logical operator '%s' at %%L are %s/%s"),
2996 gfc_op2string (e->value.op.operator), gfc_typename (&op1->ts),
2997 gfc_typename (&op2->ts));
3002 if (op1->ts.type == BT_LOGICAL)
3004 e->ts.type = BT_LOGICAL;
3005 e->ts.kind = op1->ts.kind;
3009 sprintf (msg, _("Operand of .not. operator at %%L is %s"),
3010 gfc_typename (&op1->ts));
3014 case INTRINSIC_GT_OS:
3016 case INTRINSIC_GE_OS:
3018 case INTRINSIC_LT_OS:
3020 case INTRINSIC_LE_OS:
3021 if (op1->ts.type == BT_COMPLEX || op2->ts.type == BT_COMPLEX)
3023 strcpy (msg, _("COMPLEX quantities cannot be compared at %L"));
3027 /* Fall through... */
3030 case INTRINSIC_EQ_OS:
3032 case INTRINSIC_NE_OS:
3033 if (op1->ts.type == BT_CHARACTER && op2->ts.type == BT_CHARACTER)
3035 e->ts.type = BT_LOGICAL;
3036 e->ts.kind = gfc_default_logical_kind;
3040 if (gfc_numeric_ts (&op1->ts) && gfc_numeric_ts (&op2->ts))
3042 gfc_type_convert_binary (e);
3044 e->ts.type = BT_LOGICAL;
3045 e->ts.kind = gfc_default_logical_kind;
3049 if (op1->ts.type == BT_LOGICAL && op2->ts.type == BT_LOGICAL)
3051 _("Logicals at %%L must be compared with %s instead of %s"),
3052 (e->value.op.operator == INTRINSIC_EQ
3053 || e->value.op.operator == INTRINSIC_EQ_OS)
3054 ? ".eqv." : ".neqv.", gfc_op2string (e->value.op.operator));
3057 _("Operands of comparison operator '%s' at %%L are %s/%s"),
3058 gfc_op2string (e->value.op.operator), gfc_typename (&op1->ts),
3059 gfc_typename (&op2->ts));
3063 case INTRINSIC_USER:
3064 if (e->value.op.uop->operator == NULL)
3065 sprintf (msg, _("Unknown operator '%s' at %%L"), e->value.op.uop->name);
3066 else if (op2 == NULL)
3067 sprintf (msg, _("Operand of user operator '%s' at %%L is %s"),
3068 e->value.op.uop->name, gfc_typename (&op1->ts));
3070 sprintf (msg, _("Operands of user operator '%s' at %%L are %s/%s"),
3071 e->value.op.uop->name, gfc_typename (&op1->ts),
3072 gfc_typename (&op2->ts));
3076 case INTRINSIC_PARENTHESES:
3078 if (e->ts.type == BT_CHARACTER)
3079 e->ts.cl = op1->ts.cl;
3083 gfc_internal_error ("resolve_operator(): Bad intrinsic");
3086 /* Deal with arrayness of an operand through an operator. */
3090 switch (e->value.op.operator)
3092 case INTRINSIC_PLUS:
3093 case INTRINSIC_MINUS:
3094 case INTRINSIC_TIMES:
3095 case INTRINSIC_DIVIDE:
3096 case INTRINSIC_POWER:
3097 case INTRINSIC_CONCAT:
3101 case INTRINSIC_NEQV:
3103 case INTRINSIC_EQ_OS:
3105 case INTRINSIC_NE_OS:
3107 case INTRINSIC_GT_OS:
3109 case INTRINSIC_GE_OS:
3111 case INTRINSIC_LT_OS:
3113 case INTRINSIC_LE_OS:
3115 if (op1->rank == 0 && op2->rank == 0)
3118 if (op1->rank == 0 && op2->rank != 0)
3120 e->rank = op2->rank;
3122 if (e->shape == NULL)
3123 e->shape = gfc_copy_shape (op2->shape, op2->rank);
3126 if (op1->rank != 0 && op2->rank == 0)
3128 e->rank = op1->rank;
3130 if (e->shape == NULL)
3131 e->shape = gfc_copy_shape (op1->shape, op1->rank);
3134 if (op1->rank != 0 && op2->rank != 0)
3136 if (op1->rank == op2->rank)
3138 e->rank = op1->rank;
3139 if (e->shape == NULL)
3141 t = compare_shapes(op1, op2);
3145 e->shape = gfc_copy_shape (op1->shape, op1->rank);
3150 /* Allow higher level expressions to work. */
3153 /* Try user-defined operators, and otherwise throw an error. */
3154 dual_locus_error = true;
3156 _("Inconsistent ranks for operator at %%L and %%L"));
3163 case INTRINSIC_PARENTHESES:
3165 case INTRINSIC_UPLUS:
3166 case INTRINSIC_UMINUS:
3167 /* Simply copy arrayness attribute */
3168 e->rank = op1->rank;
3170 if (e->shape == NULL)
3171 e->shape = gfc_copy_shape (op1->shape, op1->rank);
3179 /* Attempt to simplify the expression. */
3182 t = gfc_simplify_expr (e, 0);
3183 /* Some calls do not succeed in simplification and return FAILURE
3184 even though there is no error; eg. variable references to
3185 PARAMETER arrays. */
3186 if (!gfc_is_constant_expr (e))
3193 if (gfc_extend_expr (e) == SUCCESS)
3196 if (dual_locus_error)
3197 gfc_error (msg, &op1->where, &op2->where);
3199 gfc_error (msg, &e->where);
3205 /************** Array resolution subroutines **************/
3208 { CMP_LT, CMP_EQ, CMP_GT, CMP_UNKNOWN }
3211 /* Compare two integer expressions. */
3214 compare_bound (gfc_expr *a, gfc_expr *b)
3218 if (a == NULL || a->expr_type != EXPR_CONSTANT
3219 || b == NULL || b->expr_type != EXPR_CONSTANT)
3222 /* If either of the types isn't INTEGER, we must have
3223 raised an error earlier. */
3225 if (a->ts.type != BT_INTEGER || b->ts.type != BT_INTEGER)
3228 i = mpz_cmp (a->value.integer, b->value.integer);
3238 /* Compare an integer expression with an integer. */
3241 compare_bound_int (gfc_expr *a, int b)
3245 if (a == NULL || a->expr_type != EXPR_CONSTANT)
3248 if (a->ts.type != BT_INTEGER)
3249 gfc_internal_error ("compare_bound_int(): Bad expression");
3251 i = mpz_cmp_si (a->value.integer, b);
3261 /* Compare an integer expression with a mpz_t. */
3264 compare_bound_mpz_t (gfc_expr *a, mpz_t b)
3268 if (a == NULL || a->expr_type != EXPR_CONSTANT)
3271 if (a->ts.type != BT_INTEGER)
3272 gfc_internal_error ("compare_bound_int(): Bad expression");
3274 i = mpz_cmp (a->value.integer, b);
3284 /* Compute the last value of a sequence given by a triplet.
3285 Return 0 if it wasn't able to compute the last value, or if the
3286 sequence if empty, and 1 otherwise. */
3289 compute_last_value_for_triplet (gfc_expr *start, gfc_expr *end,
3290 gfc_expr *stride, mpz_t last)
3294 if (start == NULL || start->expr_type != EXPR_CONSTANT
3295 || end == NULL || end->expr_type != EXPR_CONSTANT
3296 || (stride != NULL && stride->expr_type != EXPR_CONSTANT))
3299 if (start->ts.type != BT_INTEGER || end->ts.type != BT_INTEGER
3300 || (stride != NULL && stride->ts.type != BT_INTEGER))
3303 if (stride == NULL || compare_bound_int(stride, 1) == CMP_EQ)
3305 if (compare_bound (start, end) == CMP_GT)
3307 mpz_set (last, end->value.integer);
3311 if (compare_bound_int (stride, 0) == CMP_GT)
3313 /* Stride is positive */
3314 if (mpz_cmp (start->value.integer, end->value.integer) > 0)
3319 /* Stride is negative */
3320 if (mpz_cmp (start->value.integer, end->value.integer) < 0)
3325 mpz_sub (rem, end->value.integer, start->value.integer);
3326 mpz_tdiv_r (rem, rem, stride->value.integer);
3327 mpz_sub (last, end->value.integer, rem);
3334 /* Compare a single dimension of an array reference to the array
3338 check_dimension (int i, gfc_array_ref *ar, gfc_array_spec *as)
3342 /* Given start, end and stride values, calculate the minimum and
3343 maximum referenced indexes. */
3345 switch (ar->dimen_type[i])
3351 if (compare_bound (ar->start[i], as->lower[i]) == CMP_LT)
3353 gfc_warning ("Array reference at %L is out of bounds "
3354 "(%ld < %ld) in dimension %d", &ar->c_where[i],
3355 mpz_get_si (ar->start[i]->value.integer),
3356 mpz_get_si (as->lower[i]->value.integer), i+1);
3359 if (compare_bound (ar->start[i], as->upper[i]) == CMP_GT)
3361 gfc_warning ("Array reference at %L is out of bounds "
3362 "(%ld > %ld) in dimension %d", &ar->c_where[i],
3363 mpz_get_si (ar->start[i]->value.integer),
3364 mpz_get_si (as->upper[i]->value.integer), i+1);
3372 #define AR_START (ar->start[i] ? ar->start[i] : as->lower[i])
3373 #define AR_END (ar->end[i] ? ar->end[i] : as->upper[i])
3375 comparison comp_start_end = compare_bound (AR_START, AR_END);
3377 /* Check for zero stride, which is not allowed. */
3378 if (compare_bound_int (ar->stride[i], 0) == CMP_EQ)
3380 gfc_error ("Illegal stride of zero at %L", &ar->c_where[i]);
3384 /* if start == len || (stride > 0 && start < len)
3385 || (stride < 0 && start > len),
3386 then the array section contains at least one element. In this
3387 case, there is an out-of-bounds access if
3388 (start < lower || start > upper). */
3389 if (compare_bound (AR_START, AR_END) == CMP_EQ
3390 || ((compare_bound_int (ar->stride[i], 0) == CMP_GT
3391 || ar->stride[i] == NULL) && comp_start_end == CMP_LT)
3392 || (compare_bound_int (ar->stride[i], 0) == CMP_LT
3393 && comp_start_end == CMP_GT))
3395 if (compare_bound (AR_START, as->lower[i]) == CMP_LT)
3397 gfc_warning ("Lower array reference at %L is out of bounds "
3398 "(%ld < %ld) in dimension %d", &ar->c_where[i],
3399 mpz_get_si (AR_START->value.integer),
3400 mpz_get_si (as->lower[i]->value.integer), i+1);
3403 if (compare_bound (AR_START, as->upper[i]) == CMP_GT)
3405 gfc_warning ("Lower array reference at %L is out of bounds "
3406 "(%ld > %ld) in dimension %d", &ar->c_where[i],
3407 mpz_get_si (AR_START->value.integer),
3408 mpz_get_si (as->upper[i]->value.integer), i+1);
3413 /* If we can compute the highest index of the array section,
3414 then it also has to be between lower and upper. */
3415 mpz_init (last_value);
3416 if (compute_last_value_for_triplet (AR_START, AR_END, ar->stride[i],
3419 if (compare_bound_mpz_t (as->lower[i], last_value) == CMP_GT)
3421 gfc_warning ("Upper array reference at %L is out of bounds "
3422 "(%ld < %ld) in dimension %d", &ar->c_where[i],
3423 mpz_get_si (last_value),
3424 mpz_get_si (as->lower[i]->value.integer), i+1);
3425 mpz_clear (last_value);
3428 if (compare_bound_mpz_t (as->upper[i], last_value) == CMP_LT)
3430 gfc_warning ("Upper array reference at %L is out of bounds "
3431 "(%ld > %ld) in dimension %d", &ar->c_where[i],
3432 mpz_get_si (last_value),
3433 mpz_get_si (as->upper[i]->value.integer), i+1);
3434 mpz_clear (last_value);
3438 mpz_clear (last_value);
3446 gfc_internal_error ("check_dimension(): Bad array reference");
3453 /* Compare an array reference with an array specification. */
3456 compare_spec_to_ref (gfc_array_ref *ar)
3463 /* TODO: Full array sections are only allowed as actual parameters. */
3464 if (as->type == AS_ASSUMED_SIZE
3465 && (/*ar->type == AR_FULL
3466 ||*/ (ar->type == AR_SECTION
3467 && ar->dimen_type[i] == DIMEN_RANGE && ar->end[i] == NULL)))
3469 gfc_error ("Rightmost upper bound of assumed size array section "
3470 "not specified at %L", &ar->where);
3474 if (ar->type == AR_FULL)
3477 if (as->rank != ar->dimen)
3479 gfc_error ("Rank mismatch in array reference at %L (%d/%d)",
3480 &ar->where, ar->dimen, as->rank);
3484 for (i = 0; i < as->rank; i++)
3485 if (check_dimension (i, ar, as) == FAILURE)
3492 /* Resolve one part of an array index. */
3495 gfc_resolve_index (gfc_expr *index, int check_scalar)
3502 if (gfc_resolve_expr (index) == FAILURE)
3505 if (check_scalar && index->rank != 0)
3507 gfc_error ("Array index at %L must be scalar", &index->where);
3511 if (index->ts.type != BT_INTEGER && index->ts.type != BT_REAL)
3513 gfc_error ("Array index at %L must be of INTEGER type",
3518 if (index->ts.type == BT_REAL)
3519 if (gfc_notify_std (GFC_STD_LEGACY, "Extension: REAL array index at %L",
3520 &index->where) == FAILURE)
3523 if (index->ts.kind != gfc_index_integer_kind
3524 || index->ts.type != BT_INTEGER)
3527 ts.type = BT_INTEGER;
3528 ts.kind = gfc_index_integer_kind;
3530 gfc_convert_type_warn (index, &ts, 2, 0);
3536 /* Resolve a dim argument to an intrinsic function. */
3539 gfc_resolve_dim_arg (gfc_expr *dim)
3544 if (gfc_resolve_expr (dim) == FAILURE)
3549 gfc_error ("Argument dim at %L must be scalar", &dim->where);
3554 if (dim->ts.type != BT_INTEGER)
3556 gfc_error ("Argument dim at %L must be of INTEGER type", &dim->where);
3560 if (dim->ts.kind != gfc_index_integer_kind)
3564 ts.type = BT_INTEGER;
3565 ts.kind = gfc_index_integer_kind;
3567 gfc_convert_type_warn (dim, &ts, 2, 0);
3573 /* Given an expression that contains array references, update those array
3574 references to point to the right array specifications. While this is
3575 filled in during matching, this information is difficult to save and load
3576 in a module, so we take care of it here.
3578 The idea here is that the original array reference comes from the
3579 base symbol. We traverse the list of reference structures, setting
3580 the stored reference to references. Component references can
3581 provide an additional array specification. */
3584 find_array_spec (gfc_expr *e)
3588 gfc_symbol *derived;
3591 as = e->symtree->n.sym->as;
3594 for (ref = e->ref; ref; ref = ref->next)
3599 gfc_internal_error ("find_array_spec(): Missing spec");
3606 if (derived == NULL)
3607 derived = e->symtree->n.sym->ts.derived;
3609 c = derived->components;
3611 for (; c; c = c->next)
3612 if (c == ref->u.c.component)
3614 /* Track the sequence of component references. */
3615 if (c->ts.type == BT_DERIVED)
3616 derived = c->ts.derived;
3621 gfc_internal_error ("find_array_spec(): Component not found");
3626 gfc_internal_error ("find_array_spec(): unused as(1)");
3637 gfc_internal_error ("find_array_spec(): unused as(2)");
3641 /* Resolve an array reference. */
3644 resolve_array_ref (gfc_array_ref *ar)
3646 int i, check_scalar;
3649 for (i = 0; i < ar->dimen; i++)
3651 check_scalar = ar->dimen_type[i] == DIMEN_RANGE;
3653 if (gfc_resolve_index (ar->start[i], check_scalar) == FAILURE)
3655 if (gfc_resolve_index (ar->end[i], check_scalar) == FAILURE)
3657 if (gfc_resolve_index (ar->stride[i], check_scalar) == FAILURE)
3662 if (ar->dimen_type[i] == DIMEN_UNKNOWN)
3666 ar->dimen_type[i] = DIMEN_ELEMENT;
3670 ar->dimen_type[i] = DIMEN_VECTOR;
3671 if (e->expr_type == EXPR_VARIABLE
3672 && e->symtree->n.sym->ts.type == BT_DERIVED)
3673 ar->start[i] = gfc_get_parentheses (e);
3677 gfc_error ("Array index at %L is an array of rank %d",
3678 &ar->c_where[i], e->rank);
3683 /* If the reference type is unknown, figure out what kind it is. */
3685 if (ar->type == AR_UNKNOWN)
3687 ar->type = AR_ELEMENT;
3688 for (i = 0; i < ar->dimen; i++)
3689 if (ar->dimen_type[i] == DIMEN_RANGE
3690 || ar->dimen_type[i] == DIMEN_VECTOR)
3692 ar->type = AR_SECTION;
3697 if (!ar->as->cray_pointee && compare_spec_to_ref (ar) == FAILURE)
3705 resolve_substring (gfc_ref *ref)
3707 if (ref->u.ss.start != NULL)
3709 if (gfc_resolve_expr (ref->u.ss.start) == FAILURE)
3712 if (ref->u.ss.start->ts.type != BT_INTEGER)
3714 gfc_error ("Substring start index at %L must be of type INTEGER",
3715 &ref->u.ss.start->where);
3719 if (ref->u.ss.start->rank != 0)
3721 gfc_error ("Substring start index at %L must be scalar",
3722 &ref->u.ss.start->where);
3726 if (compare_bound_int (ref->u.ss.start, 1) == CMP_LT
3727 && (compare_bound (ref->u.ss.end, ref->u.ss.start) == CMP_EQ
3728 || compare_bound (ref->u.ss.end, ref->u.ss.start) == CMP_GT))
3730 gfc_error ("Substring start index at %L is less than one",
3731 &ref->u.ss.start->where);
3736 if (ref->u.ss.end != NULL)
3738 if (gfc_resolve_expr (ref->u.ss.end) == FAILURE)
3741 if (ref->u.ss.end->ts.type != BT_INTEGER)
3743 gfc_error ("Substring end index at %L must be of type INTEGER",
3744 &ref->u.ss.end->where);
3748 if (ref->u.ss.end->rank != 0)
3750 gfc_error ("Substring end index at %L must be scalar",
3751 &ref->u.ss.end->where);
3755 if (ref->u.ss.length != NULL
3756 && compare_bound (ref->u.ss.end, ref->u.ss.length->length) == CMP_GT
3757 && (compare_bound (ref->u.ss.end, ref->u.ss.start) == CMP_EQ
3758 || compare_bound (ref->u.ss.end, ref->u.ss.start) == CMP_GT))
3760 gfc_error ("Substring end index at %L exceeds the string length",
3761 &ref->u.ss.start->where);
3770 /* This function supplies missing substring charlens. */
3773 gfc_resolve_substring_charlen (gfc_expr *e)
3776 gfc_expr *start, *end;
3778 for (char_ref = e->ref; char_ref; char_ref = char_ref->next)
3779 if (char_ref->type == REF_SUBSTRING)
3785 gcc_assert (char_ref->next == NULL);
3789 if (e->ts.cl->length)
3790 gfc_free_expr (e->ts.cl->length);
3791 else if (e->expr_type == EXPR_VARIABLE
3792 && e->symtree->n.sym->attr.dummy)
3796 e->ts.type = BT_CHARACTER;
3797 e->ts.kind = gfc_default_character_kind;
3801 e->ts.cl = gfc_get_charlen ();
3802 e->ts.cl->next = gfc_current_ns->cl_list;
3803 gfc_current_ns->cl_list = e->ts.cl;
3806 if (char_ref->u.ss.start)
3807 start = gfc_copy_expr (char_ref->u.ss.start);
3809 start = gfc_int_expr (1);
3811 if (char_ref->u.ss.end)
3812 end = gfc_copy_expr (char_ref->u.ss.end);
3813 else if (e->expr_type == EXPR_VARIABLE)
3814 end = gfc_copy_expr (e->symtree->n.sym->ts.cl->length);
3821 /* Length = (end - start +1). */
3822 e->ts.cl->length = gfc_subtract (end, start);
3823 e->ts.cl->length = gfc_add (e->ts.cl->length, gfc_int_expr (1));
3825 e->ts.cl->length->ts.type = BT_INTEGER;
3826 e->ts.cl->length->ts.kind = gfc_charlen_int_kind;;
3828 /* Make sure that the length is simplified. */
3829 gfc_simplify_expr (e->ts.cl->length, 1);
3830 gfc_resolve_expr (e->ts.cl->length);
3834 /* Resolve subtype references. */
3837 resolve_ref (gfc_expr *expr)
3839 int current_part_dimension, n_components, seen_part_dimension;
3842 for (ref = expr->ref; ref; ref = ref->next)
3843 if (ref->type == REF_ARRAY && ref->u.ar.as == NULL)
3845 find_array_spec (expr);
3849 for (ref = expr->ref; ref; ref = ref->next)
3853 if (resolve_array_ref (&ref->u.ar) == FAILURE)
3861 resolve_substring (ref);
3865 /* Check constraints on part references. */
3867 current_part_dimension = 0;
3868 seen_part_dimension = 0;
3871 for (ref = expr->ref; ref; ref = ref->next)
3876 switch (ref->u.ar.type)
3880 current_part_dimension = 1;
3884 current_part_dimension = 0;
3888 gfc_internal_error ("resolve_ref(): Bad array reference");
3894 if (current_part_dimension || seen_part_dimension)
3896 if (ref->u.c.component->pointer)
3898 gfc_error ("Component to the right of a part reference "
3899 "with nonzero rank must not have the POINTER "
3900 "attribute at %L", &expr->where);
3903 else if (ref->u.c.component->allocatable)
3905 gfc_error ("Component to the right of a part reference "
3906 "with nonzero rank must not have the ALLOCATABLE "
3907 "attribute at %L", &expr->where);
3919 if (((ref->type == REF_COMPONENT && n_components > 1)
3920 || ref->next == NULL)
3921 && current_part_dimension
3922 && seen_part_dimension)
3924 gfc_error ("Two or more part references with nonzero rank must "
3925 "not be specified at %L", &expr->where);
3929 if (ref->type == REF_COMPONENT)
3931 if (current_part_dimension)
3932 seen_part_dimension = 1;
3934 /* reset to make sure */
3935 current_part_dimension = 0;
3943 /* Given an expression, determine its shape. This is easier than it sounds.
3944 Leaves the shape array NULL if it is not possible to determine the shape. */
3947 expression_shape (gfc_expr *e)
3949 mpz_t array[GFC_MAX_DIMENSIONS];
3952 if (e->rank == 0 || e->shape != NULL)
3955 for (i = 0; i < e->rank; i++)
3956 if (gfc_array_dimen_size (e, i, &array[i]) == FAILURE)
3959 e->shape = gfc_get_shape (e->rank);
3961 memcpy (e->shape, array, e->rank * sizeof (mpz_t));
3966 for (i--; i >= 0; i--)
3967 mpz_clear (array[i]);
3971 /* Given a variable expression node, compute the rank of the expression by
3972 examining the base symbol and any reference structures it may have. */
3975 expression_rank (gfc_expr *e)
3982 if (e->expr_type == EXPR_ARRAY)
3984 /* Constructors can have a rank different from one via RESHAPE(). */
3986 if (e->symtree == NULL)
3992 e->rank = (e->symtree->n.sym->as == NULL)
3993 ? 0 : e->symtree->n.sym->as->rank;
3999 for (ref = e->ref; ref; ref = ref->next)
4001 if (ref->type != REF_ARRAY)
4004 if (ref->u.ar.type == AR_FULL)
4006 rank = ref->u.ar.as->rank;
4010 if (ref->u.ar.type == AR_SECTION)
4012 /* Figure out the rank of the section. */
4014 gfc_internal_error ("expression_rank(): Two array specs");
4016 for (i = 0; i < ref->u.ar.dimen; i++)
4017 if (ref->u.ar.dimen_type[i] == DIMEN_RANGE
4018 || ref->u.ar.dimen_type[i] == DIMEN_VECTOR)
4028 expression_shape (e);
4032 /* Resolve a variable expression. */
4035 resolve_variable (gfc_expr *e)
4042 if (e->symtree == NULL)
4045 if (e->ref && resolve_ref (e) == FAILURE)
4048 sym = e->symtree->n.sym;
4049 if (sym->attr.flavor == FL_PROCEDURE && !sym->attr.function)
4051 e->ts.type = BT_PROCEDURE;
4055 if (sym->ts.type != BT_UNKNOWN)
4056 gfc_variable_attr (e, &e->ts);
4059 /* Must be a simple variable reference. */
4060 if (gfc_set_default_type (sym, 1, sym->ns) == FAILURE)
4065 if (check_assumed_size_reference (sym, e))
4068 /* Deal with forward references to entries during resolve_code, to
4069 satisfy, at least partially, 12.5.2.5. */
4070 if (gfc_current_ns->entries
4071 && current_entry_id == sym->entry_id
4074 && cs_base->current->op != EXEC_ENTRY)
4076 gfc_entry_list *entry;
4077 gfc_formal_arglist *formal;
4081 /* If the symbol is a dummy... */
4082 if (sym->attr.dummy && sym->ns == gfc_current_ns)
4084 entry = gfc_current_ns->entries;
4087 /* ...test if the symbol is a parameter of previous entries. */
4088 for (; entry && entry->id <= current_entry_id; entry = entry->next)
4089 for (formal = entry->sym->formal; formal; formal = formal->next)
4091 if (formal->sym && sym->name == formal->sym->name)
4095 /* If it has not been seen as a dummy, this is an error. */
4098 if (specification_expr)
4099 gfc_error ("Variable '%s', used in a specification expression"
4100 ", is referenced at %L before the ENTRY statement "
4101 "in which it is a parameter",
4102 sym->name, &cs_base->current->loc);
4104 gfc_error ("Variable '%s' is used at %L before the ENTRY "
4105 "statement in which it is a parameter",
4106 sym->name, &cs_base->current->loc);
4111 /* Now do the same check on the specification expressions. */
4112 specification_expr = 1;
4113 if (sym->ts.type == BT_CHARACTER
4114 && gfc_resolve_expr (sym->ts.cl->length) == FAILURE)
4118 for (n = 0; n < sym->as->rank; n++)
4120 specification_expr = 1;
4121 if (gfc_resolve_expr (sym->as->lower[n]) == FAILURE)
4123 specification_expr = 1;
4124 if (gfc_resolve_expr (sym->as->upper[n]) == FAILURE)
4127 specification_expr = 0;
4130 /* Update the symbol's entry level. */
4131 sym->entry_id = current_entry_id + 1;
4138 /* Checks to see that the correct symbol has been host associated.
4139 The only situation where this arises is that in which a twice
4140 contained function is parsed after the host association is made.
4141 Therefore, on detecting this, the line is rematched, having got
4142 rid of the existing references and actual_arg_list. */
4144 check_host_association (gfc_expr *e)
4146 gfc_symbol *sym, *old_sym;
4150 bool retval = e->expr_type == EXPR_FUNCTION;
4152 if (e->symtree == NULL || e->symtree->n.sym == NULL)
4155 old_sym = e->symtree->n.sym;
4157 if (old_sym->attr.use_assoc)
4160 if (gfc_current_ns->parent
4161 && old_sym->ns != gfc_current_ns)
4163 gfc_find_symbol (old_sym->name, gfc_current_ns, 1, &sym);
4164 if (sym && old_sym != sym
4165 && sym->attr.flavor == FL_PROCEDURE
4166 && sym->attr.contained)
4168 temp_locus = gfc_current_locus;
4169 gfc_current_locus = e->where;
4171 gfc_buffer_error (1);
4173 gfc_free_ref_list (e->ref);
4178 gfc_free_actual_arglist (e->value.function.actual);
4179 e->value.function.actual = NULL;
4182 if (e->shape != NULL)
4184 for (n = 0; n < e->rank; n++)
4185 mpz_clear (e->shape[n]);
4187 gfc_free (e->shape);
4190 gfc_match_rvalue (&expr);
4192 gfc_buffer_error (0);
4194 gcc_assert (expr && sym == expr->symtree->n.sym);
4200 gfc_current_locus = temp_locus;
4203 /* This might have changed! */
4204 return e->expr_type == EXPR_FUNCTION;
4209 gfc_resolve_character_operator (gfc_expr *e)
4211 gfc_expr *op1 = e->value.op.op1;
4212 gfc_expr *op2 = e->value.op.op2;
4213 gfc_expr *e1 = NULL;
4214 gfc_expr *e2 = NULL;
4216 gcc_assert (e->value.op.operator == INTRINSIC_CONCAT);
4218 if (op1->ts.cl && op1->ts.cl->length)
4219 e1 = gfc_copy_expr (op1->ts.cl->length);
4220 else if (op1->expr_type == EXPR_CONSTANT)
4221 e1 = gfc_int_expr (op1->value.character.length);
4223 if (op2->ts.cl && op2->ts.cl->length)
4224 e2 = gfc_copy_expr (op2->ts.cl->length);
4225 else if (op2->expr_type == EXPR_CONSTANT)
4226 e2 = gfc_int_expr (op2->value.character.length);
4228 e->ts.cl = gfc_get_charlen ();
4229 e->ts.cl->next = gfc_current_ns->cl_list;
4230 gfc_current_ns->cl_list = e->ts.cl;
4235 e->ts.cl->length = gfc_add (e1, e2);
4236 e->ts.cl->length->ts.type = BT_INTEGER;
4237 e->ts.cl->length->ts.kind = gfc_charlen_int_kind;;
4238 gfc_simplify_expr (e->ts.cl->length, 0);
4239 gfc_resolve_expr (e->ts.cl->length);
4245 /* Ensure that an character expression has a charlen and, if possible, a
4246 length expression. */
4249 fixup_charlen (gfc_expr *e)
4251 /* The cases fall through so that changes in expression type and the need
4252 for multiple fixes are picked up. In all circumstances, a charlen should
4253 be available for the middle end to hang a backend_decl on. */
4254 switch (e->expr_type)
4257 gfc_resolve_character_operator (e);
4260 if (e->expr_type == EXPR_ARRAY)
4261 gfc_resolve_character_array_constructor (e);
4263 case EXPR_SUBSTRING:
4264 if (!e->ts.cl && e->ref)
4265 gfc_resolve_substring_charlen (e);
4270 e->ts.cl = gfc_get_charlen ();
4271 e->ts.cl->next = gfc_current_ns->cl_list;
4272 gfc_current_ns->cl_list = e->ts.cl;
4280 /* Resolve an expression. That is, make sure that types of operands agree
4281 with their operators, intrinsic operators are converted to function calls
4282 for overloaded types and unresolved function references are resolved. */
4285 gfc_resolve_expr (gfc_expr *e)
4292 switch (e->expr_type)
4295 t = resolve_operator (e);
4301 if (check_host_association (e))
4302 t = resolve_function (e);
4305 t = resolve_variable (e);
4307 expression_rank (e);
4310 if (e->ts.type == BT_CHARACTER && e->ts.cl == NULL && e->ref
4311 && e->ref->type != REF_SUBSTRING)
4312 gfc_resolve_substring_charlen (e);
4316 case EXPR_SUBSTRING:
4317 t = resolve_ref (e);
4327 if (resolve_ref (e) == FAILURE)
4330 t = gfc_resolve_array_constructor (e);
4331 /* Also try to expand a constructor. */
4334 expression_rank (e);
4335 gfc_expand_constructor (e);
4338 /* This provides the opportunity for the length of constructors with
4339 character valued function elements to propagate the string length
4340 to the expression. */
4341 if (e->ts.type == BT_CHARACTER)
4342 gfc_resolve_character_array_constructor (e);
4346 case EXPR_STRUCTURE:
4347 t = resolve_ref (e);
4351 t = resolve_structure_cons (e);
4355 t = gfc_simplify_expr (e, 0);
4359 gfc_internal_error ("gfc_resolve_expr(): Bad expression type");
4362 if (e->ts.type == BT_CHARACTER && t == SUCCESS && !e->ts.cl)
4369 /* Resolve an expression from an iterator. They must be scalar and have
4370 INTEGER or (optionally) REAL type. */
4373 gfc_resolve_iterator_expr (gfc_expr *expr, bool real_ok,
4374 const char *name_msgid)
4376 if (gfc_resolve_expr (expr) == FAILURE)
4379 if (expr->rank != 0)
4381 gfc_error ("%s at %L must be a scalar", _(name_msgid), &expr->where);
4385 if (expr->ts.type != BT_INTEGER)
4387 if (expr->ts.type == BT_REAL)
4390 return gfc_notify_std (GFC_STD_F95_DEL,
4391 "Deleted feature: %s at %L must be integer",
4392 _(name_msgid), &expr->where);
4395 gfc_error ("%s at %L must be INTEGER", _(name_msgid),
4402 gfc_error ("%s at %L must be INTEGER", _(name_msgid), &expr->where);
4410 /* Resolve the expressions in an iterator structure. If REAL_OK is
4411 false allow only INTEGER type iterators, otherwise allow REAL types. */
4414 gfc_resolve_iterator (gfc_iterator *iter, bool real_ok)
4416 if (gfc_resolve_iterator_expr (iter->var, real_ok, "Loop variable")
4420 if (gfc_pure (NULL) && gfc_impure_variable (iter->var->symtree->n.sym))
4422 gfc_error ("Cannot assign to loop variable in PURE procedure at %L",
4427 if (gfc_resolve_iterator_expr (iter->start, real_ok,
4428 "Start expression in DO loop") == FAILURE)
4431 if (gfc_resolve_iterator_expr (iter->end, real_ok,
4432 "End expression in DO loop") == FAILURE)
4435 if (gfc_resolve_iterator_expr (iter->step, real_ok,
4436 "Step expression in DO loop") == FAILURE)
4439 if (iter->step->expr_type == EXPR_CONSTANT)
4441 if ((iter->step->ts.type == BT_INTEGER
4442 && mpz_cmp_ui (iter->step->value.integer, 0) == 0)
4443 || (iter->step->ts.type == BT_REAL
4444 && mpfr_sgn (iter->step->value.real) == 0))
4446 gfc_error ("Step expression in DO loop at %L cannot be zero",
4447 &iter->step->where);
4452 /* Convert start, end, and step to the same type as var. */
4453 if (iter->start->ts.kind != iter->var->ts.kind
4454 || iter->start->ts.type != iter->var->ts.type)
4455 gfc_convert_type (iter->start, &iter->var->ts, 2);
4457 if (iter->end->ts.kind != iter->var->ts.kind
4458 || iter->end->ts.type != iter->var->ts.type)
4459 gfc_convert_type (iter->end, &iter->var->ts, 2);
4461 if (iter->step->ts.kind != iter->var->ts.kind
4462 || iter->step->ts.type != iter->var->ts.type)
4463 gfc_convert_type (iter->step, &iter->var->ts, 2);
4469 /* Traversal function for find_forall_index. f == 2 signals that
4470 that variable itself is not to be checked - only the references. */
4473 forall_index (gfc_expr *expr, gfc_symbol *sym, int *f)
4475 if (expr->expr_type != EXPR_VARIABLE)
4478 /* A scalar assignment */
4479 if (!expr->ref || *f == 1)
4481 if (expr->symtree->n.sym == sym)
4493 /* Check whether the FORALL index appears in the expression or not.
4494 Returns SUCCESS if SYM is found in EXPR. */
4497 find_forall_index (gfc_expr *expr, gfc_symbol *sym, int f)
4499 if (gfc_traverse_expr (expr, sym, forall_index, f))
4506 /* Resolve a list of FORALL iterators. The FORALL index-name is constrained
4507 to be a scalar INTEGER variable. The subscripts and stride are scalar
4508 INTEGERs, and if stride is a constant it must be nonzero.
4509 Furthermore "A subscript or stride in a forall-triplet-spec shall
4510 not contain a reference to any index-name in the
4511 forall-triplet-spec-list in which it appears." (7.5.4.1) */
4514 resolve_forall_iterators (gfc_forall_iterator *it)
4516 gfc_forall_iterator *iter, *iter2;
4518 for (iter = it; iter; iter = iter->next)
4520 if (gfc_resolve_expr (iter->var) == SUCCESS
4521 && (iter->var->ts.type != BT_INTEGER || iter->var->rank != 0))
4522 gfc_error ("FORALL index-name at %L must be a scalar INTEGER",
4525 if (gfc_resolve_expr (iter->start) == SUCCESS
4526 && (iter->start->ts.type != BT_INTEGER || iter->start->rank != 0))
4527 gfc_error ("FORALL start expression at %L must be a scalar INTEGER",
4528 &iter->start->where);
4529 if (iter->var->ts.kind != iter->start->ts.kind)
4530 gfc_convert_type (iter->start, &iter->var->ts, 2);
4532 if (gfc_resolve_expr (iter->end) == SUCCESS
4533 && (iter->end->ts.type != BT_INTEGER || iter->end->rank != 0))
4534 gfc_error ("FORALL end expression at %L must be a scalar INTEGER",
4536 if (iter->var->ts.kind != iter->end->ts.kind)
4537 gfc_convert_type (iter->end, &iter->var->ts, 2);
4539 if (gfc_resolve_expr (iter->stride) == SUCCESS)
4541 if (iter->stride->ts.type != BT_INTEGER || iter->stride->rank != 0)
4542 gfc_error ("FORALL stride expression at %L must be a scalar %s",
4543 &iter->stride->where, "INTEGER");
4545 if (iter->stride->expr_type == EXPR_CONSTANT
4546 && mpz_cmp_ui(iter->stride->value.integer, 0) == 0)
4547 gfc_error ("FORALL stride expression at %L cannot be zero",
4548 &iter->stride->where);
4550 if (iter->var->ts.kind != iter->stride->ts.kind)
4551 gfc_convert_type (iter->stride, &iter->var->ts, 2);
4554 for (iter = it; iter; iter = iter->next)
4555 for (iter2 = iter; iter2; iter2 = iter2->next)
4557 if (find_forall_index (iter2->start,
4558 iter->var->symtree->n.sym, 0) == SUCCESS
4559 || find_forall_index (iter2->end,
4560 iter->var->symtree->n.sym, 0) == SUCCESS
4561 || find_forall_index (iter2->stride,
4562 iter->var->symtree->n.sym, 0) == SUCCESS)
4563 gfc_error ("FORALL index '%s' may not appear in triplet "
4564 "specification at %L", iter->var->symtree->name,
4565 &iter2->start->where);
4570 /* Given a pointer to a symbol that is a derived type, see if it's
4571 inaccessible, i.e. if it's defined in another module and the components are
4572 PRIVATE. The search is recursive if necessary. Returns zero if no
4573 inaccessible components are found, nonzero otherwise. */
4576 derived_inaccessible (gfc_symbol *sym)
4580 if (sym->attr.use_assoc && sym->attr.private_comp)
4583 for (c = sym->components; c; c = c->next)
4585 if (c->ts.type == BT_DERIVED && derived_inaccessible (c->ts.derived))
4593 /* Resolve the argument of a deallocate expression. The expression must be
4594 a pointer or a full array. */
4597 resolve_deallocate_expr (gfc_expr *e)
4599 symbol_attribute attr;
4600 int allocatable, pointer, check_intent_in;
4603 /* Check INTENT(IN), unless the object is a sub-component of a pointer. */
4604 check_intent_in = 1;
4606 if (gfc_resolve_expr (e) == FAILURE)
4609 if (e->expr_type != EXPR_VARIABLE)
4612 allocatable = e->symtree->n.sym->attr.allocatable;
4613 pointer = e->symtree->n.sym->attr.pointer;
4614 for (ref = e->ref; ref; ref = ref->next)
4617 check_intent_in = 0;
4622 if (ref->u.ar.type != AR_FULL)
4627 allocatable = (ref->u.c.component->as != NULL
4628 && ref->u.c.component->as->type == AS_DEFERRED);
4629 pointer = ref->u.c.component->pointer;
4638 attr = gfc_expr_attr (e);
4640 if (allocatable == 0 && attr.pointer == 0)
4643 gfc_error ("Expression in DEALLOCATE statement at %L must be "
4644 "ALLOCATABLE or a POINTER", &e->where);
4648 && e->symtree->n.sym->attr.intent == INTENT_IN)
4650 gfc_error ("Cannot deallocate INTENT(IN) variable '%s' at %L",
4651 e->symtree->n.sym->name, &e->where);
4659 /* Returns true if the expression e contains a reference to the symbol sym. */
4661 sym_in_expr (gfc_expr *e, gfc_symbol *sym, int *f ATTRIBUTE_UNUSED)
4663 if (e->expr_type == EXPR_VARIABLE && e->symtree->n.sym == sym)
4670 find_sym_in_expr (gfc_symbol *sym, gfc_expr *e)
4672 return gfc_traverse_expr (e, sym, sym_in_expr, 0);
4676 /* Given the expression node e for an allocatable/pointer of derived type to be
4677 allocated, get the expression node to be initialized afterwards (needed for
4678 derived types with default initializers, and derived types with allocatable
4679 components that need nullification.) */
4682 expr_to_initialize (gfc_expr *e)
4688 result = gfc_copy_expr (e);
4690 /* Change the last array reference from AR_ELEMENT to AR_FULL. */
4691 for (ref = result->ref; ref; ref = ref->next)
4692 if (ref->type == REF_ARRAY && ref->next == NULL)
4694 ref->u.ar.type = AR_FULL;
4696 for (i = 0; i < ref->u.ar.dimen; i++)
4697 ref->u.ar.start[i] = ref->u.ar.end[i] = ref->u.ar.stride[i] = NULL;
4699 result->rank = ref->u.ar.dimen;
4707 /* Resolve the expression in an ALLOCATE statement, doing the additional
4708 checks to see whether the expression is OK or not. The expression must
4709 have a trailing array reference that gives the size of the array. */
4712 resolve_allocate_expr (gfc_expr *e, gfc_code *code)
4714 int i, pointer, allocatable, dimension, check_intent_in;
4715 symbol_attribute attr;
4716 gfc_ref *ref, *ref2;
4723 /* Check INTENT(IN), unless the object is a sub-component of a pointer. */
4724 check_intent_in = 1;
4726 if (gfc_resolve_expr (e) == FAILURE)
4729 if (code->expr && code->expr->expr_type == EXPR_VARIABLE)
4730 sym = code->expr->symtree->n.sym;
4734 /* Make sure the expression is allocatable or a pointer. If it is
4735 pointer, the next-to-last reference must be a pointer. */
4739 if (e->expr_type != EXPR_VARIABLE)
4742 attr = gfc_expr_attr (e);
4743 pointer = attr.pointer;
4744 dimension = attr.dimension;
4748 allocatable = e->symtree->n.sym->attr.allocatable;
4749 pointer = e->symtree->n.sym->attr.pointer;
4750 dimension = e->symtree->n.sym->attr.dimension;
4752 if (sym == e->symtree->n.sym && sym->ts.type != BT_DERIVED)
4754 gfc_error ("The STAT variable '%s' in an ALLOCATE statement must "
4755 "not be allocated in the same statement at %L",
4756 sym->name, &e->where);
4760 for (ref = e->ref; ref; ref2 = ref, ref = ref->next)
4763 check_intent_in = 0;
4768 if (ref->next != NULL)
4773 allocatable = (ref->u.c.component->as != NULL
4774 && ref->u.c.component->as->type == AS_DEFERRED);
4776 pointer = ref->u.c.component->pointer;
4777 dimension = ref->u.c.component->dimension;
4788 if (allocatable == 0 && pointer == 0)
4790 gfc_error ("Expression in ALLOCATE statement at %L must be "
4791 "ALLOCATABLE or a POINTER", &e->where);
4796 && e->symtree->n.sym->attr.intent == INTENT_IN)
4798 gfc_error ("Cannot allocate INTENT(IN) variable '%s' at %L",
4799 e->symtree->n.sym->name, &e->where);
4803 /* Add default initializer for those derived types that need them. */
4804 if (e->ts.type == BT_DERIVED && (init_e = gfc_default_initializer (&e->ts)))
4806 init_st = gfc_get_code ();
4807 init_st->loc = code->loc;
4808 init_st->op = EXEC_INIT_ASSIGN;
4809 init_st->expr = expr_to_initialize (e);
4810 init_st->expr2 = init_e;
4811 init_st->next = code->next;
4812 code->next = init_st;
4815 if (pointer && dimension == 0)
4818 /* Make sure the next-to-last reference node is an array specification. */
4820 if (ref2 == NULL || ref2->type != REF_ARRAY || ref2->u.ar.type == AR_FULL)
4822 gfc_error ("Array specification required in ALLOCATE statement "
4823 "at %L", &e->where);
4827 /* Make sure that the array section reference makes sense in the
4828 context of an ALLOCATE specification. */
4832 for (i = 0; i < ar->dimen; i++)
4834 if (ref2->u.ar.type == AR_ELEMENT)
4837 switch (ar->dimen_type[i])
4843 if (ar->start[i] != NULL
4844 && ar->end[i] != NULL
4845 && ar->stride[i] == NULL)
4848 /* Fall Through... */
4852 gfc_error ("Bad array specification in ALLOCATE statement at %L",
4859 for (a = code->ext.alloc_list; a; a = a->next)
4861 sym = a->expr->symtree->n.sym;
4863 /* TODO - check derived type components. */
4864 if (sym->ts.type == BT_DERIVED)
4867 if ((ar->start[i] != NULL && find_sym_in_expr (sym, ar->start[i]))
4868 || (ar->end[i] != NULL && find_sym_in_expr (sym, ar->end[i])))
4870 gfc_error ("'%s' must not appear an the array specification at "
4871 "%L in the same ALLOCATE statement where it is "
4872 "itself allocated", sym->name, &ar->where);
4882 resolve_allocate_deallocate (gfc_code *code, const char *fcn)
4884 gfc_symbol *s = NULL;
4888 s = code->expr->symtree->n.sym;
4892 if (s->attr.intent == INTENT_IN)
4893 gfc_error ("STAT variable '%s' of %s statement at %C cannot "
4894 "be INTENT(IN)", s->name, fcn);
4896 if (gfc_pure (NULL) && gfc_impure_variable (s))
4897 gfc_error ("Illegal STAT variable in %s statement at %C "
4898 "for a PURE procedure", fcn);
4901 if (s && code->expr->ts.type != BT_INTEGER)
4902 gfc_error ("STAT tag in %s statement at %L must be "
4903 "of type INTEGER", fcn, &code->expr->where);
4905 if (strcmp (fcn, "ALLOCATE") == 0)
4907 for (a = code->ext.alloc_list; a; a = a->next)
4908 resolve_allocate_expr (a->expr, code);
4912 for (a = code->ext.alloc_list; a; a = a->next)
4913 resolve_deallocate_expr (a->expr);
4917 /************ SELECT CASE resolution subroutines ************/
4919 /* Callback function for our mergesort variant. Determines interval
4920 overlaps for CASEs. Return <0 if op1 < op2, 0 for overlap, >0 for
4921 op1 > op2. Assumes we're not dealing with the default case.
4922 We have op1 = (:L), (K:L) or (K:) and op2 = (:N), (M:N) or (M:).
4923 There are nine situations to check. */
4926 compare_cases (const gfc_case *op1, const gfc_case *op2)
4930 if (op1->low == NULL) /* op1 = (:L) */
4932 /* op2 = (:N), so overlap. */
4934 /* op2 = (M:) or (M:N), L < M */
4935 if (op2->low != NULL
4936 && gfc_compare_expr (op1->high, op2->low, INTRINSIC_LT) < 0)
4939 else if (op1->high == NULL) /* op1 = (K:) */
4941 /* op2 = (M:), so overlap. */
4943 /* op2 = (:N) or (M:N), K > N */
4944 if (op2->high != NULL
4945 && gfc_compare_expr (op1->low, op2->high, INTRINSIC_GT) > 0)
4948 else /* op1 = (K:L) */
4950 if (op2->low == NULL) /* op2 = (:N), K > N */
4951 retval = (gfc_compare_expr (op1->low, op2->high, INTRINSIC_GT) > 0)
4953 else if (op2->high == NULL) /* op2 = (M:), L < M */
4954 retval = (gfc_compare_expr (op1->high, op2->low, INTRINSIC_LT) < 0)
4956 else /* op2 = (M:N) */
4960 if (gfc_compare_expr (op1->high, op2->low, INTRINSIC_LT) < 0)
4963 else if (gfc_compare_expr (op1->low, op2->high, INTRINSIC_GT) > 0)
4972 /* Merge-sort a double linked case list, detecting overlap in the
4973 process. LIST is the head of the double linked case list before it
4974 is sorted. Returns the head of the sorted list if we don't see any
4975 overlap, or NULL otherwise. */
4978 check_case_overlap (gfc_case *list)
4980 gfc_case *p, *q, *e, *tail;
4981 int insize, nmerges, psize, qsize, cmp, overlap_seen;
4983 /* If the passed list was empty, return immediately. */
4990 /* Loop unconditionally. The only exit from this loop is a return
4991 statement, when we've finished sorting the case list. */
4998 /* Count the number of merges we do in this pass. */
5001 /* Loop while there exists a merge to be done. */
5006 /* Count this merge. */
5009 /* Cut the list in two pieces by stepping INSIZE places
5010 forward in the list, starting from P. */
5013 for (i = 0; i < insize; i++)
5022 /* Now we have two lists. Merge them! */
5023 while (psize > 0 || (qsize > 0 && q != NULL))
5025 /* See from which the next case to merge comes from. */
5028 /* P is empty so the next case must come from Q. */
5033 else if (qsize == 0 || q == NULL)
5042 cmp = compare_cases (p, q);
5045 /* The whole case range for P is less than the
5053 /* The whole case range for Q is greater than
5054 the case range for P. */
5061 /* The cases overlap, or they are the same
5062 element in the list. Either way, we must
5063 issue an error and get the next case from P. */
5064 /* FIXME: Sort P and Q by line number. */
5065 gfc_error ("CASE label at %L overlaps with CASE "
5066 "label at %L", &p->where, &q->where);
5074 /* Add the next element to the merged list. */
5083 /* P has now stepped INSIZE places along, and so has Q. So
5084 they're the same. */
5089 /* If we have done only one merge or none at all, we've
5090 finished sorting the cases. */
5099 /* Otherwise repeat, merging lists twice the size. */
5105 /* Check to see if an expression is suitable for use in a CASE statement.
5106 Makes sure that all case expressions are scalar constants of the same
5107 type. Return FAILURE if anything is wrong. */
5110 validate_case_label_expr (gfc_expr *e, gfc_expr *case_expr)
5112 if (e == NULL) return SUCCESS;
5114 if (e->ts.type != case_expr->ts.type)
5116 gfc_error ("Expression in CASE statement at %L must be of type %s",
5117 &e->where, gfc_basic_typename (case_expr->ts.type));
5121 /* C805 (R808) For a given case-construct, each case-value shall be of
5122 the same type as case-expr. For character type, length differences
5123 are allowed, but the kind type parameters shall be the same. */
5125 if (case_expr->ts.type == BT_CHARACTER && e->ts.kind != case_expr->ts.kind)
5127 gfc_error("Expression in CASE statement at %L must be kind %d",
5128 &e->where, case_expr->ts.kind);
5132 /* Convert the case value kind to that of case expression kind, if needed.
5133 FIXME: Should a warning be issued? */
5134 if (e->ts.kind != case_expr->ts.kind)
5135 gfc_convert_type_warn (e, &case_expr->ts, 2, 0);
5139 gfc_error ("Expression in CASE statement at %L must be scalar",
5148 /* Given a completely parsed select statement, we:
5150 - Validate all expressions and code within the SELECT.
5151 - Make sure that the selection expression is not of the wrong type.
5152 - Make sure that no case ranges overlap.
5153 - Eliminate unreachable cases and unreachable code resulting from
5154 removing case labels.
5156 The standard does allow unreachable cases, e.g. CASE (5:3). But
5157 they are a hassle for code generation, and to prevent that, we just
5158 cut them out here. This is not necessary for overlapping cases
5159 because they are illegal and we never even try to generate code.
5161 We have the additional caveat that a SELECT construct could have
5162 been a computed GOTO in the source code. Fortunately we can fairly
5163 easily work around that here: The case_expr for a "real" SELECT CASE
5164 is in code->expr1, but for a computed GOTO it is in code->expr2. All
5165 we have to do is make sure that the case_expr is a scalar integer
5169 resolve_select (gfc_code *code)
5172 gfc_expr *case_expr;
5173 gfc_case *cp, *default_case, *tail, *head;
5174 int seen_unreachable;
5180 if (code->expr == NULL)
5182 /* This was actually a computed GOTO statement. */
5183 case_expr = code->expr2;
5184 if (case_expr->ts.type != BT_INTEGER|| case_expr->rank != 0)
5185 gfc_error ("Selection expression in computed GOTO statement "
5186 "at %L must be a scalar integer expression",
5189 /* Further checking is not necessary because this SELECT was built
5190 by the compiler, so it should always be OK. Just move the
5191 case_expr from expr2 to expr so that we can handle computed
5192 GOTOs as normal SELECTs from here on. */
5193 code->expr = code->expr2;
5198 case_expr = code->expr;
5200 type = case_expr->ts.type;
5201 if (type != BT_LOGICAL && type != BT_INTEGER && type != BT_CHARACTER)
5203 gfc_error ("Argument of SELECT statement at %L cannot be %s",
5204 &case_expr->where, gfc_typename (&case_expr->ts));
5206 /* Punt. Going on here just produce more garbage error messages. */
5210 if (case_expr->rank != 0)
5212 gfc_error ("Argument of SELECT statement at %L must be a scalar "
5213 "expression", &case_expr->where);
5219 /* PR 19168 has a long discussion concerning a mismatch of the kinds
5220 of the SELECT CASE expression and its CASE values. Walk the lists
5221 of case values, and if we find a mismatch, promote case_expr to
5222 the appropriate kind. */
5224 if (type == BT_LOGICAL || type == BT_INTEGER)
5226 for (body = code->block; body; body = body->block)
5228 /* Walk the case label list. */
5229 for (cp = body->ext.case_list; cp; cp = cp->next)
5231 /* Intercept the DEFAULT case. It does not have a kind. */
5232 if (cp->low == NULL && cp->high == NULL)
5235 /* Unreachable case ranges are discarded, so ignore. */
5236 if (cp->low != NULL && cp->high != NULL
5237 && cp->low != cp->high
5238 && gfc_compare_expr (cp->low, cp->high, INTRINSIC_GT) > 0)
5241 /* FIXME: Should a warning be issued? */
5243 && case_expr->ts.kind != gfc_kind_max(case_expr, cp->low))
5244 gfc_convert_type_warn (case_expr, &cp->low->ts, 2, 0);
5246 if (cp->high != NULL
5247 && case_expr->ts.kind != gfc_kind_max(case_expr, cp->high))
5248 gfc_convert_type_warn (case_expr, &cp->high->ts, 2, 0);
5253 /* Assume there is no DEFAULT case. */
5254 default_case = NULL;
5259 for (body = code->block; body; body = body->block)
5261 /* Assume the CASE list is OK, and all CASE labels can be matched. */
5263 seen_unreachable = 0;
5265 /* Walk the case label list, making sure that all case labels
5267 for (cp = body->ext.case_list; cp; cp = cp->next)
5269 /* Count the number of cases in the whole construct. */
5272 /* Intercept the DEFAULT case. */
5273 if (cp->low == NULL && cp->high == NULL)
5275 if (default_case != NULL)
5277 gfc_error ("The DEFAULT CASE at %L cannot be followed "
5278 "by a second DEFAULT CASE at %L",
5279 &default_case->where, &cp->where);
5290 /* Deal with single value cases and case ranges. Errors are
5291 issued from the validation function. */
5292 if(validate_case_label_expr (cp->low, case_expr) != SUCCESS
5293 || validate_case_label_expr (cp->high, case_expr) != SUCCESS)
5299 if (type == BT_LOGICAL
5300 && ((cp->low == NULL || cp->high == NULL)
5301 || cp->low != cp->high))
5303 gfc_error ("Logical range in CASE statement at %L is not "
5304 "allowed", &cp->low->where);
5309 if (type == BT_LOGICAL && cp->low->expr_type == EXPR_CONSTANT)
5312 value = cp->low->value.logical == 0 ? 2 : 1;
5313 if (value & seen_logical)
5315 gfc_error ("constant logical value in CASE statement "
5316 "is repeated at %L",
5321 seen_logical |= value;
5324 if (cp->low != NULL && cp->high != NULL
5325 && cp->low != cp->high
5326 && gfc_compare_expr (cp->low, cp->high, INTRINSIC_GT) > 0)
5328 if (gfc_option.warn_surprising)
5329 gfc_warning ("Range specification at %L can never "
5330 "be matched", &cp->where);
5332 cp->unreachable = 1;
5333 seen_unreachable = 1;
5337 /* If the case range can be matched, it can also overlap with
5338 other cases. To make sure it does not, we put it in a
5339 double linked list here. We sort that with a merge sort
5340 later on to detect any overlapping cases. */
5344 head->right = head->left = NULL;
5349 tail->right->left = tail;
5356 /* It there was a failure in the previous case label, give up
5357 for this case label list. Continue with the next block. */
5361 /* See if any case labels that are unreachable have been seen.
5362 If so, we eliminate them. This is a bit of a kludge because
5363 the case lists for a single case statement (label) is a
5364 single forward linked lists. */
5365 if (seen_unreachable)
5367 /* Advance until the first case in the list is reachable. */
5368 while (body->ext.case_list != NULL
5369 && body->ext.case_list->unreachable)
5371 gfc_case *n = body->ext.case_list;
5372 body->ext.case_list = body->ext.case_list->next;
5374 gfc_free_case_list (n);
5377 /* Strip all other unreachable cases. */
5378 if (body->ext.case_list)
5380 for (cp = body->ext.case_list; cp->next; cp = cp->next)
5382 if (cp->next->unreachable)
5384 gfc_case *n = cp->next;
5385 cp->next = cp->next->next;
5387 gfc_free_case_list (n);
5394 /* See if there were overlapping cases. If the check returns NULL,
5395 there was overlap. In that case we don't do anything. If head
5396 is non-NULL, we prepend the DEFAULT case. The sorted list can
5397 then used during code generation for SELECT CASE constructs with
5398 a case expression of a CHARACTER type. */
5401 head = check_case_overlap (head);
5403 /* Prepend the default_case if it is there. */
5404 if (head != NULL && default_case)
5406 default_case->left = NULL;
5407 default_case->right = head;
5408 head->left = default_case;
5412 /* Eliminate dead blocks that may be the result if we've seen
5413 unreachable case labels for a block. */
5414 for (body = code; body && body->block; body = body->block)
5416 if (body->block->ext.case_list == NULL)
5418 /* Cut the unreachable block from the code chain. */
5419 gfc_code *c = body->block;
5420 body->block = c->block;
5422 /* Kill the dead block, but not the blocks below it. */
5424 gfc_free_statements (c);
5428 /* More than two cases is legal but insane for logical selects.
5429 Issue a warning for it. */
5430 if (gfc_option.warn_surprising && type == BT_LOGICAL
5432 gfc_warning ("Logical SELECT CASE block at %L has more that two cases",
5437 /* Resolve a transfer statement. This is making sure that:
5438 -- a derived type being transferred has only non-pointer components
5439 -- a derived type being transferred doesn't have private components, unless
5440 it's being transferred from the module where the type was defined
5441 -- we're not trying to transfer a whole assumed size array. */
5444 resolve_transfer (gfc_code *code)
5453 if (exp->expr_type != EXPR_VARIABLE && exp->expr_type != EXPR_FUNCTION)
5456 sym = exp->symtree->n.sym;
5459 /* Go to actual component transferred. */
5460 for (ref = code->expr->ref; ref; ref = ref->next)
5461 if (ref->type == REF_COMPONENT)
5462 ts = &ref->u.c.component->ts;
5464 if (ts->type == BT_DERIVED)
5466 /* Check that transferred derived type doesn't contain POINTER
5468 if (ts->derived->attr.pointer_comp)
5470 gfc_error ("Data transfer element at %L cannot have "
5471 "POINTER components", &code->loc);
5475 if (ts->derived->attr.alloc_comp)
5477 gfc_error ("Data transfer element at %L cannot have "
5478 "ALLOCATABLE components", &code->loc);
5482 if (derived_inaccessible (ts->derived))
5484 gfc_error ("Data transfer element at %L cannot have "
5485 "PRIVATE components",&code->loc);
5490 if (sym->as != NULL && sym->as->type == AS_ASSUMED_SIZE
5491 && exp->ref->type == REF_ARRAY && exp->ref->u.ar.type == AR_FULL)
5493 gfc_error ("Data transfer element at %L cannot be a full reference to "
5494 "an assumed-size array", &code->loc);
5500 /*********** Toplevel code resolution subroutines ***********/
5502 /* Find the set of labels that are reachable from this block. We also
5503 record the last statement in each block so that we don't have to do
5504 a linear search to find the END DO statements of the blocks. */
5507 reachable_labels (gfc_code *block)
5514 cs_base->reachable_labels = bitmap_obstack_alloc (&labels_obstack);
5516 /* Collect labels in this block. */
5517 for (c = block; c; c = c->next)
5520 bitmap_set_bit (cs_base->reachable_labels, c->here->value);
5522 if (!c->next && cs_base->prev)
5523 cs_base->prev->tail = c;
5526 /* Merge with labels from parent block. */
5529 gcc_assert (cs_base->prev->reachable_labels);
5530 bitmap_ior_into (cs_base->reachable_labels,
5531 cs_base->prev->reachable_labels);
5535 /* Given a branch to a label and a namespace, if the branch is conforming.
5536 The code node describes where the branch is located. */
5539 resolve_branch (gfc_st_label *label, gfc_code *code)
5546 /* Step one: is this a valid branching target? */
5548 if (label->defined == ST_LABEL_UNKNOWN)
5550 gfc_error ("Label %d referenced at %L is never defined", label->value,
5555 if (label->defined != ST_LABEL_TARGET)
5557 gfc_error ("Statement at %L is not a valid branch target statement "
5558 "for the branch statement at %L", &label->where, &code->loc);
5562 /* Step two: make sure this branch is not a branch to itself ;-) */
5564 if (code->here == label)
5566 gfc_warning ("Branch at %L may result in an infinite loop", &code->loc);
5570 /* Step three: See if the label is in the same block as the
5571 branching statement. The hard work has been done by setting up
5572 the bitmap reachable_labels. */
5574 if (!bitmap_bit_p (cs_base->reachable_labels, label->value))
5576 /* The label is not in an enclosing block, so illegal. This was
5577 allowed in Fortran 66, so we allow it as extension. No
5578 further checks are necessary in this case. */
5579 gfc_notify_std (GFC_STD_LEGACY, "Label at %L is not in the same block "
5580 "as the GOTO statement at %L", &label->where,
5585 /* Step four: Make sure that the branching target is legal if
5586 the statement is an END {SELECT,IF}. */
5588 for (stack = cs_base; stack; stack = stack->prev)
5589 if (stack->current->next && stack->current->next->here == label)
5592 if (stack && stack->current->next->op == EXEC_NOP)
5594 gfc_notify_std (GFC_STD_F95_DEL, "Deleted feature: GOTO at %L jumps to "
5595 "END of construct at %L", &code->loc,
5596 &stack->current->next->loc);
5597 return; /* We know this is not an END DO. */
5600 /* Step five: Make sure that we're not jumping to the end of a DO
5601 loop from within the loop. */
5603 for (stack = cs_base; stack; stack = stack->prev)
5604 if ((stack->current->op == EXEC_DO
5605 || stack->current->op == EXEC_DO_WHILE)
5606 && stack->tail->here == label && stack->tail->op == EXEC_NOP)
5608 gfc_notify_std (GFC_STD_F95_DEL, "Deleted feature: GOTO at %L jumps "
5609 "to END of construct at %L", &code->loc,
5617 /* Check whether EXPR1 has the same shape as EXPR2. */
5620 resolve_where_shape (gfc_expr *expr1, gfc_expr *expr2)
5622 mpz_t shape[GFC_MAX_DIMENSIONS];
5623 mpz_t shape2[GFC_MAX_DIMENSIONS];
5624 try result = FAILURE;
5627 /* Compare the rank. */
5628 if (expr1->rank != expr2->rank)
5631 /* Compare the size of each dimension. */
5632 for (i=0; i<expr1->rank; i++)
5634 if (gfc_array_dimen_size (expr1, i, &shape[i]) == FAILURE)
5637 if (gfc_array_dimen_size (expr2, i, &shape2[i]) == FAILURE)
5640 if (mpz_cmp (shape[i], shape2[i]))
5644 /* When either of the two expression is an assumed size array, we
5645 ignore the comparison of dimension sizes. */
5650 for (i--; i >= 0; i--)
5652 mpz_clear (shape[i]);
5653 mpz_clear (shape2[i]);
5659 /* Check whether a WHERE assignment target or a WHERE mask expression
5660 has the same shape as the outmost WHERE mask expression. */
5663 resolve_where (gfc_code *code, gfc_expr *mask)
5669 cblock = code->block;
5671 /* Store the first WHERE mask-expr of the WHERE statement or construct.
5672 In case of nested WHERE, only the outmost one is stored. */
5673 if (mask == NULL) /* outmost WHERE */
5675 else /* inner WHERE */
5682 /* Check if the mask-expr has a consistent shape with the
5683 outmost WHERE mask-expr. */
5684 if (resolve_where_shape (cblock->expr, e) == FAILURE)
5685 gfc_error ("WHERE mask at %L has inconsistent shape",
5686 &cblock->expr->where);
5689 /* the assignment statement of a WHERE statement, or the first
5690 statement in where-body-construct of a WHERE construct */
5691 cnext = cblock->next;
5696 /* WHERE assignment statement */
5699 /* Check shape consistent for WHERE assignment target. */
5700 if (e && resolve_where_shape (cnext->expr, e) == FAILURE)
5701 gfc_error ("WHERE assignment target at %L has "
5702 "inconsistent shape", &cnext->expr->where);
5706 case EXEC_ASSIGN_CALL:
5707 resolve_call (cnext);
5708 if (!cnext->resolved_sym->attr.elemental)
5709 gfc_error("Non-ELEMENTAL user-defined assignment in WHERE at %L",
5710 &cnext->ext.actual->expr->where);
5713 /* WHERE or WHERE construct is part of a where-body-construct */
5715 resolve_where (cnext, e);
5719 gfc_error ("Unsupported statement inside WHERE at %L",
5722 /* the next statement within the same where-body-construct */
5723 cnext = cnext->next;
5725 /* the next masked-elsewhere-stmt, elsewhere-stmt, or end-where-stmt */
5726 cblock = cblock->block;
5731 /* Resolve assignment in FORALL construct.
5732 NVAR is the number of FORALL index variables, and VAR_EXPR records the
5733 FORALL index variables. */
5736 gfc_resolve_assign_in_forall (gfc_code *code, int nvar, gfc_expr **var_expr)
5740 for (n = 0; n < nvar; n++)
5742 gfc_symbol *forall_index;
5744 forall_index = var_expr[n]->symtree->n.sym;
5746 /* Check whether the assignment target is one of the FORALL index
5748 if ((code->expr->expr_type == EXPR_VARIABLE)
5749 && (code->expr->symtree->n.sym == forall_index))
5750 gfc_error ("Assignment to a FORALL index variable at %L",
5751 &code->expr->where);
5754 /* If one of the FORALL index variables doesn't appear in the
5755 assignment target, then there will be a many-to-one
5757 if (find_forall_index (code->expr, forall_index, 0) == FAILURE)
5758 gfc_error ("The FORALL with index '%s' cause more than one "
5759 "assignment to this object at %L",
5760 var_expr[n]->symtree->name, &code->expr->where);
5766 /* Resolve WHERE statement in FORALL construct. */
5769 gfc_resolve_where_code_in_forall (gfc_code *code, int nvar,
5770 gfc_expr **var_expr)
5775 cblock = code->block;
5778 /* the assignment statement of a WHERE statement, or the first
5779 statement in where-body-construct of a WHERE construct */
5780 cnext = cblock->next;
5785 /* WHERE assignment statement */
5787 gfc_resolve_assign_in_forall (cnext, nvar, var_expr);
5790 /* WHERE operator assignment statement */
5791 case EXEC_ASSIGN_CALL:
5792 resolve_call (cnext);
5793 if (!cnext->resolved_sym->attr.elemental)
5794 gfc_error("Non-ELEMENTAL user-defined assignment in WHERE at %L",
5795 &cnext->ext.actual->expr->where);
5798 /* WHERE or WHERE construct is part of a where-body-construct */
5800 gfc_resolve_where_code_in_forall (cnext, nvar, var_expr);
5804 gfc_error ("Unsupported statement inside WHERE at %L",
5807 /* the next statement within the same where-body-construct */
5808 cnext = cnext->next;
5810 /* the next masked-elsewhere-stmt, elsewhere-stmt, or end-where-stmt */
5811 cblock = cblock->block;
5816 /* Traverse the FORALL body to check whether the following errors exist:
5817 1. For assignment, check if a many-to-one assignment happens.
5818 2. For WHERE statement, check the WHERE body to see if there is any
5819 many-to-one assignment. */
5822 gfc_resolve_forall_body (gfc_code *code, int nvar, gfc_expr **var_expr)
5826 c = code->block->next;
5832 case EXEC_POINTER_ASSIGN:
5833 gfc_resolve_assign_in_forall (c, nvar, var_expr);
5836 case EXEC_ASSIGN_CALL:
5840 /* Because the gfc_resolve_blocks() will handle the nested FORALL,
5841 there is no need to handle it here. */
5845 gfc_resolve_where_code_in_forall(c, nvar, var_expr);
5850 /* The next statement in the FORALL body. */
5856 /* Given a FORALL construct, first resolve the FORALL iterator, then call
5857 gfc_resolve_forall_body to resolve the FORALL body. */
5860 gfc_resolve_forall (gfc_code *code, gfc_namespace *ns, int forall_save)
5862 static gfc_expr **var_expr;
5863 static int total_var = 0;
5864 static int nvar = 0;
5865 gfc_forall_iterator *fa;
5869 /* Start to resolve a FORALL construct */
5870 if (forall_save == 0)
5872 /* Count the total number of FORALL index in the nested FORALL
5873 construct in order to allocate the VAR_EXPR with proper size. */
5875 while ((next != NULL) && (next->op == EXEC_FORALL))
5877 for (fa = next->ext.forall_iterator; fa; fa = fa->next)
5879 next = next->block->next;
5882 /* Allocate VAR_EXPR with NUMBER_OF_FORALL_INDEX elements. */
5883 var_expr = (gfc_expr **) gfc_getmem (total_var * sizeof (gfc_expr *));
5886 /* The information about FORALL iterator, including FORALL index start, end
5887 and stride. The FORALL index can not appear in start, end or stride. */
5888 for (fa = code->ext.forall_iterator; fa; fa = fa->next)
5890 /* Check if any outer FORALL index name is the same as the current
5892 for (i = 0; i < nvar; i++)
5894 if (fa->var->symtree->n.sym == var_expr[i]->symtree->n.sym)
5896 gfc_error ("An outer FORALL construct already has an index "
5897 "with this name %L", &fa->var->where);
5901 /* Record the current FORALL index. */
5902 var_expr[nvar] = gfc_copy_expr (fa->var);
5907 /* Resolve the FORALL body. */
5908 gfc_resolve_forall_body (code, nvar, var_expr);
5910 /* May call gfc_resolve_forall to resolve the inner FORALL loop. */
5911 gfc_resolve_blocks (code->block, ns);
5913 /* Free VAR_EXPR after the whole FORALL construct resolved. */
5914 for (i = 0; i < total_var; i++)
5915 gfc_free_expr (var_expr[i]);
5917 /* Reset the counters. */
5923 /* Resolve lists of blocks found in IF, SELECT CASE, WHERE, FORALL ,GOTO and
5926 static void resolve_code (gfc_code *, gfc_namespace *);
5929 gfc_resolve_blocks (gfc_code *b, gfc_namespace *ns)
5933 for (; b; b = b->block)
5935 t = gfc_resolve_expr (b->expr);
5936 if (gfc_resolve_expr (b->expr2) == FAILURE)
5942 if (t == SUCCESS && b->expr != NULL
5943 && (b->expr->ts.type != BT_LOGICAL || b->expr->rank != 0))
5944 gfc_error ("IF clause at %L requires a scalar LOGICAL expression",
5951 && (b->expr->ts.type != BT_LOGICAL || b->expr->rank == 0))
5952 gfc_error ("WHERE/ELSEWHERE clause at %L requires a LOGICAL array",
5957 resolve_branch (b->label, b);
5970 case EXEC_OMP_ATOMIC:
5971 case EXEC_OMP_CRITICAL:
5973 case EXEC_OMP_MASTER:
5974 case EXEC_OMP_ORDERED:
5975 case EXEC_OMP_PARALLEL:
5976 case EXEC_OMP_PARALLEL_DO:
5977 case EXEC_OMP_PARALLEL_SECTIONS:
5978 case EXEC_OMP_PARALLEL_WORKSHARE:
5979 case EXEC_OMP_SECTIONS:
5980 case EXEC_OMP_SINGLE:
5981 case EXEC_OMP_WORKSHARE:
5985 gfc_internal_error ("resolve_block(): Bad block type");
5988 resolve_code (b->next, ns);
5993 /* Does everything to resolve an ordinary assignment. Returns true
5994 if this is an interface asignment. */
5996 resolve_ordinary_assign (gfc_code *code, gfc_namespace *ns)
6006 if (gfc_extend_assign (code, ns) == SUCCESS)
6008 lhs = code->ext.actual->expr;
6009 rhs = code->ext.actual->next->expr;
6010 if (gfc_pure (NULL) && !gfc_pure (code->symtree->n.sym))
6012 gfc_error ("Subroutine '%s' called instead of assignment at "
6013 "%L must be PURE", code->symtree->n.sym->name,
6018 /* Make a temporary rhs when there is a default initializer
6019 and rhs is the same symbol as the lhs. */
6020 if (rhs->expr_type == EXPR_VARIABLE
6021 && rhs->symtree->n.sym->ts.type == BT_DERIVED
6022 && has_default_initializer (rhs->symtree->n.sym->ts.derived)
6023 && (lhs->symtree->n.sym == rhs->symtree->n.sym))
6024 code->ext.actual->next->expr = gfc_get_parentheses (rhs);
6033 && gfc_notify_std (GFC_STD_GNU, "Extension: BOZ literal at %L outside "
6034 "a DATA statement and outside INT/REAL/DBLE/CMPLX",
6035 &code->loc) == FAILURE)
6038 /* Handle the case of a BOZ literal on the RHS. */
6039 if (rhs->is_boz && lhs->ts.type != BT_INTEGER)
6042 if (gfc_option.warn_surprising)
6043 gfc_warning ("BOZ literal at %L is bitwise transferred "
6044 "non-integer symbol '%s'", &code->loc,
6045 lhs->symtree->n.sym->name);
6047 if (!gfc_convert_boz (rhs, &lhs->ts))
6049 if ((rc = gfc_range_check (rhs)) != ARITH_OK)
6051 if (rc == ARITH_UNDERFLOW)
6052 gfc_error ("Arithmetic underflow of bit-wise transferred BOZ at %L"
6053 ". This check can be disabled with the option "
6054 "-fno-range-check", &rhs->where);
6055 else if (rc == ARITH_OVERFLOW)
6056 gfc_error ("Arithmetic overflow of bit-wise transferred BOZ at %L"
6057 ". This check can be disabled with the option "
6058 "-fno-range-check", &rhs->where);
6059 else if (rc == ARITH_NAN)
6060 gfc_error ("Arithmetic NaN of bit-wise transferred BOZ at %L"
6061 ". This check can be disabled with the option "
6062 "-fno-range-check", &rhs->where);
6068 if (lhs->ts.type == BT_CHARACTER
6069 && gfc_option.warn_character_truncation)
6071 if (lhs->ts.cl != NULL
6072 && lhs->ts.cl->length != NULL
6073 && lhs->ts.cl->length->expr_type == EXPR_CONSTANT)
6074 llen = mpz_get_si (lhs->ts.cl->length->value.integer);
6076 if (rhs->expr_type == EXPR_CONSTANT)
6077 rlen = rhs->value.character.length;
6079 else if (rhs->ts.cl != NULL
6080 && rhs->ts.cl->length != NULL
6081 && rhs->ts.cl->length->expr_type == EXPR_CONSTANT)
6082 rlen = mpz_get_si (rhs->ts.cl->length->value.integer);
6084 if (rlen && llen && rlen > llen)
6085 gfc_warning_now ("CHARACTER expression will be truncated "
6086 "in assignment (%d/%d) at %L",
6087 llen, rlen, &code->loc);
6090 /* Ensure that a vector index expression for the lvalue is evaluated
6091 to a temporary if the lvalue symbol is referenced in it. */
6094 for (ref = lhs->ref; ref; ref= ref->next)
6095 if (ref->type == REF_ARRAY)
6097 for (n = 0; n < ref->u.ar.dimen; n++)
6098 if (ref->u.ar.dimen_type[n] == DIMEN_VECTOR
6099 && find_sym_in_expr (lhs->symtree->n.sym,
6100 ref->u.ar.start[n]))
6102 = gfc_get_parentheses (ref->u.ar.start[n]);
6106 if (gfc_pure (NULL))
6108 if (gfc_impure_variable (lhs->symtree->n.sym))
6110 gfc_error ("Cannot assign to variable '%s' in PURE "
6112 lhs->symtree->n.sym->name,
6117 if (lhs->ts.type == BT_DERIVED
6118 && lhs->expr_type == EXPR_VARIABLE
6119 && lhs->ts.derived->attr.pointer_comp
6120 && gfc_impure_variable (rhs->symtree->n.sym))
6122 gfc_error ("The impure variable at %L is assigned to "
6123 "a derived type variable with a POINTER "
6124 "component in a PURE procedure (12.6)",
6130 gfc_check_assign (lhs, rhs, 1);
6134 /* Given a block of code, recursively resolve everything pointed to by this
6138 resolve_code (gfc_code *code, gfc_namespace *ns)
6140 int omp_workshare_save;
6145 frame.prev = cs_base;
6149 reachable_labels (code);
6151 for (; code; code = code->next)
6153 frame.current = code;
6154 forall_save = forall_flag;
6156 if (code->op == EXEC_FORALL)
6159 gfc_resolve_forall (code, ns, forall_save);
6162 else if (code->block)
6164 omp_workshare_save = -1;
6167 case EXEC_OMP_PARALLEL_WORKSHARE:
6168 omp_workshare_save = omp_workshare_flag;
6169 omp_workshare_flag = 1;
6170 gfc_resolve_omp_parallel_blocks (code, ns);
6172 case EXEC_OMP_PARALLEL:
6173 case EXEC_OMP_PARALLEL_DO:
6174 case EXEC_OMP_PARALLEL_SECTIONS:
6175 omp_workshare_save = omp_workshare_flag;
6176 omp_workshare_flag = 0;
6177 gfc_resolve_omp_parallel_blocks (code, ns);
6180 gfc_resolve_omp_do_blocks (code, ns);
6182 case EXEC_OMP_WORKSHARE:
6183 omp_workshare_save = omp_workshare_flag;
6184 omp_workshare_flag = 1;
6187 gfc_resolve_blocks (code->block, ns);
6191 if (omp_workshare_save != -1)
6192 omp_workshare_flag = omp_workshare_save;
6195 t = gfc_resolve_expr (code->expr);
6196 forall_flag = forall_save;
6198 if (gfc_resolve_expr (code->expr2) == FAILURE)
6213 /* Keep track of which entry we are up to. */
6214 current_entry_id = code->ext.entry->id;
6218 resolve_where (code, NULL);
6222 if (code->expr != NULL)
6224 if (code->expr->ts.type != BT_INTEGER)
6225 gfc_error ("ASSIGNED GOTO statement at %L requires an "
6226 "INTEGER variable", &code->expr->where);
6227 else if (code->expr->symtree->n.sym->attr.assign != 1)
6228 gfc_error ("Variable '%s' has not been assigned a target "
6229 "label at %L", code->expr->symtree->n.sym->name,
6230 &code->expr->where);
6233 resolve_branch (code->label, code);
6237 if (code->expr != NULL
6238 && (code->expr->ts.type != BT_INTEGER || code->expr->rank))
6239 gfc_error ("Alternate RETURN statement at %L requires a SCALAR-"
6240 "INTEGER return specifier", &code->expr->where);
6243 case EXEC_INIT_ASSIGN:
6250 if (resolve_ordinary_assign (code, ns))
6255 case EXEC_LABEL_ASSIGN:
6256 if (code->label->defined == ST_LABEL_UNKNOWN)
6257 gfc_error ("Label %d referenced at %L is never defined",
6258 code->label->value, &code->label->where);
6260 && (code->expr->expr_type != EXPR_VARIABLE
6261 || code->expr->symtree->n.sym->ts.type != BT_INTEGER
6262 || code->expr->symtree->n.sym->ts.kind
6263 != gfc_default_integer_kind
6264 || code->expr->symtree->n.sym->as != NULL))
6265 gfc_error ("ASSIGN statement at %L requires a scalar "
6266 "default INTEGER variable", &code->expr->where);
6269 case EXEC_POINTER_ASSIGN:
6273 gfc_check_pointer_assign (code->expr, code->expr2);
6276 case EXEC_ARITHMETIC_IF:
6278 && code->expr->ts.type != BT_INTEGER
6279 && code->expr->ts.type != BT_REAL)
6280 gfc_error ("Arithmetic IF statement at %L requires a numeric "
6281 "expression", &code->expr->where);
6283 resolve_branch (code->label, code);
6284 resolve_branch (code->label2, code);
6285 resolve_branch (code->label3, code);
6289 if (t == SUCCESS && code->expr != NULL
6290 && (code->expr->ts.type != BT_LOGICAL
6291 || code->expr->rank != 0))
6292 gfc_error ("IF clause at %L requires a scalar LOGICAL expression",
6293 &code->expr->where);
6298 resolve_call (code);
6302 /* Select is complicated. Also, a SELECT construct could be
6303 a transformed computed GOTO. */
6304 resolve_select (code);
6308 if (code->ext.iterator != NULL)
6310 gfc_iterator *iter = code->ext.iterator;
6311 if (gfc_resolve_iterator (iter, true) != FAILURE)
6312 gfc_resolve_do_iterator (code, iter->var->symtree->n.sym);
6317 if (code->expr == NULL)
6318 gfc_internal_error ("resolve_code(): No expression on DO WHILE");
6320 && (code->expr->rank != 0
6321 || code->expr->ts.type != BT_LOGICAL))
6322 gfc_error ("Exit condition of DO WHILE loop at %L must be "
6323 "a scalar LOGICAL expression", &code->expr->where);
6328 resolve_allocate_deallocate (code, "ALLOCATE");
6332 case EXEC_DEALLOCATE:
6334 resolve_allocate_deallocate (code, "DEALLOCATE");
6339 if (gfc_resolve_open (code->ext.open) == FAILURE)
6342 resolve_branch (code->ext.open->err, code);
6346 if (gfc_resolve_close (code->ext.close) == FAILURE)
6349 resolve_branch (code->ext.close->err, code);
6352 case EXEC_BACKSPACE:
6356 if (gfc_resolve_filepos (code->ext.filepos) == FAILURE)
6359 resolve_branch (code->ext.filepos->err, code);
6363 if (gfc_resolve_inquire (code->ext.inquire) == FAILURE)
6366 resolve_branch (code->ext.inquire->err, code);
6370 gcc_assert (code->ext.inquire != NULL);
6371 if (gfc_resolve_inquire (code->ext.inquire) == FAILURE)
6374 resolve_branch (code->ext.inquire->err, code);
6378 if (gfc_resolve_wait (code->ext.wait) == FAILURE)
6381 resolve_branch (code->ext.wait->err, code);
6382 resolve_branch (code->ext.wait->end, code);
6383 resolve_branch (code->ext.wait->eor, code);
6388 if (gfc_resolve_dt (code->ext.dt) == FAILURE)
6391 resolve_branch (code->ext.dt->err, code);
6392 resolve_branch (code->ext.dt->end, code);
6393 resolve_branch (code->ext.dt->eor, code);
6397 resolve_transfer (code);
6401 resolve_forall_iterators (code->ext.forall_iterator);
6403 if (code->expr != NULL && code->expr->ts.type != BT_LOGICAL)
6404 gfc_error ("FORALL mask clause at %L requires a LOGICAL "
6405 "expression", &code->expr->where);
6408 case EXEC_OMP_ATOMIC:
6409 case EXEC_OMP_BARRIER:
6410 case EXEC_OMP_CRITICAL:
6411 case EXEC_OMP_FLUSH:
6413 case EXEC_OMP_MASTER:
6414 case EXEC_OMP_ORDERED:
6415 case EXEC_OMP_SECTIONS:
6416 case EXEC_OMP_SINGLE:
6417 case EXEC_OMP_WORKSHARE:
6418 gfc_resolve_omp_directive (code, ns);
6421 case EXEC_OMP_PARALLEL:
6422 case EXEC_OMP_PARALLEL_DO:
6423 case EXEC_OMP_PARALLEL_SECTIONS:
6424 case EXEC_OMP_PARALLEL_WORKSHARE:
6425 omp_workshare_save = omp_workshare_flag;
6426 omp_workshare_flag = 0;
6427 gfc_resolve_omp_directive (code, ns);
6428 omp_workshare_flag = omp_workshare_save;
6432 gfc_internal_error ("resolve_code(): Bad statement code");
6436 cs_base = frame.prev;
6440 /* Resolve initial values and make sure they are compatible with
6444 resolve_values (gfc_symbol *sym)
6446 if (sym->value == NULL)
6449 if (gfc_resolve_expr (sym->value) == FAILURE)
6452 gfc_check_assign_symbol (sym, sym->value);
6456 /* Verify the binding labels for common blocks that are BIND(C). The label
6457 for a BIND(C) common block must be identical in all scoping units in which
6458 the common block is declared. Further, the binding label can not collide
6459 with any other global entity in the program. */
6462 resolve_bind_c_comms (gfc_symtree *comm_block_tree)
6464 if (comm_block_tree->n.common->is_bind_c == 1)
6466 gfc_gsymbol *binding_label_gsym;
6467 gfc_gsymbol *comm_name_gsym;
6469 /* See if a global symbol exists by the common block's name. It may
6470 be NULL if the common block is use-associated. */
6471 comm_name_gsym = gfc_find_gsymbol (gfc_gsym_root,
6472 comm_block_tree->n.common->name);
6473 if (comm_name_gsym != NULL && comm_name_gsym->type != GSYM_COMMON)
6474 gfc_error ("Binding label '%s' for common block '%s' at %L collides "
6475 "with the global entity '%s' at %L",
6476 comm_block_tree->n.common->binding_label,
6477 comm_block_tree->n.common->name,
6478 &(comm_block_tree->n.common->where),
6479 comm_name_gsym->name, &(comm_name_gsym->where));
6480 else if (comm_name_gsym != NULL
6481 && strcmp (comm_name_gsym->name,
6482 comm_block_tree->n.common->name) == 0)
6484 /* TODO: Need to make sure the fields of gfc_gsymbol are initialized
6486 if (comm_name_gsym->binding_label == NULL)
6487 /* No binding label for common block stored yet; save this one. */
6488 comm_name_gsym->binding_label =
6489 comm_block_tree->n.common->binding_label;
6491 if (strcmp (comm_name_gsym->binding_label,
6492 comm_block_tree->n.common->binding_label) != 0)
6494 /* Common block names match but binding labels do not. */
6495 gfc_error ("Binding label '%s' for common block '%s' at %L "
6496 "does not match the binding label '%s' for common "
6498 comm_block_tree->n.common->binding_label,
6499 comm_block_tree->n.common->name,
6500 &(comm_block_tree->n.common->where),
6501 comm_name_gsym->binding_label,
6502 comm_name_gsym->name,
6503 &(comm_name_gsym->where));
6508 /* There is no binding label (NAME="") so we have nothing further to
6509 check and nothing to add as a global symbol for the label. */
6510 if (comm_block_tree->n.common->binding_label[0] == '\0' )
6513 binding_label_gsym =
6514 gfc_find_gsymbol (gfc_gsym_root,
6515 comm_block_tree->n.common->binding_label);
6516 if (binding_label_gsym == NULL)
6518 /* Need to make a global symbol for the binding label to prevent
6519 it from colliding with another. */
6520 binding_label_gsym =
6521 gfc_get_gsymbol (comm_block_tree->n.common->binding_label);
6522 binding_label_gsym->sym_name = comm_block_tree->n.common->name;
6523 binding_label_gsym->type = GSYM_COMMON;
6527 /* If comm_name_gsym is NULL, the name common block is use
6528 associated and the name could be colliding. */
6529 if (binding_label_gsym->type != GSYM_COMMON)
6530 gfc_error ("Binding label '%s' for common block '%s' at %L "
6531 "collides with the global entity '%s' at %L",
6532 comm_block_tree->n.common->binding_label,
6533 comm_block_tree->n.common->name,
6534 &(comm_block_tree->n.common->where),
6535 binding_label_gsym->name,
6536 &(binding_label_gsym->where));
6537 else if (comm_name_gsym != NULL
6538 && (strcmp (binding_label_gsym->name,
6539 comm_name_gsym->binding_label) != 0)
6540 && (strcmp (binding_label_gsym->sym_name,
6541 comm_name_gsym->name) != 0))
6542 gfc_error ("Binding label '%s' for common block '%s' at %L "
6543 "collides with global entity '%s' at %L",
6544 binding_label_gsym->name, binding_label_gsym->sym_name,
6545 &(comm_block_tree->n.common->where),
6546 comm_name_gsym->name, &(comm_name_gsym->where));
6554 /* Verify any BIND(C) derived types in the namespace so we can report errors
6555 for them once, rather than for each variable declared of that type. */
6558 resolve_bind_c_derived_types (gfc_symbol *derived_sym)
6560 if (derived_sym != NULL && derived_sym->attr.flavor == FL_DERIVED
6561 && derived_sym->attr.is_bind_c == 1)
6562 verify_bind_c_derived_type (derived_sym);
6568 /* Verify that any binding labels used in a given namespace do not collide
6569 with the names or binding labels of any global symbols. */
6572 gfc_verify_binding_labels (gfc_symbol *sym)
6576 if (sym != NULL && sym->attr.is_bind_c && sym->attr.is_iso_c == 0
6577 && sym->attr.flavor != FL_DERIVED && sym->binding_label[0] != '\0')
6579 gfc_gsymbol *bind_c_sym;
6581 bind_c_sym = gfc_find_gsymbol (gfc_gsym_root, sym->binding_label);
6582 if (bind_c_sym != NULL
6583 && strcmp (bind_c_sym->name, sym->binding_label) == 0)
6585 if (sym->attr.if_source == IFSRC_DECL
6586 && (bind_c_sym->type != GSYM_SUBROUTINE
6587 && bind_c_sym->type != GSYM_FUNCTION)
6588 && ((sym->attr.contained == 1
6589 && strcmp (bind_c_sym->sym_name, sym->name) != 0)
6590 || (sym->attr.use_assoc == 1
6591 && (strcmp (bind_c_sym->mod_name, sym->module) != 0))))
6593 /* Make sure global procedures don't collide with anything. */
6594 gfc_error ("Binding label '%s' at %L collides with the global "
6595 "entity '%s' at %L", sym->binding_label,
6596 &(sym->declared_at), bind_c_sym->name,
6597 &(bind_c_sym->where));
6600 else if (sym->attr.contained == 0
6601 && (sym->attr.if_source == IFSRC_IFBODY
6602 && sym->attr.flavor == FL_PROCEDURE)
6603 && (bind_c_sym->sym_name != NULL
6604 && strcmp (bind_c_sym->sym_name, sym->name) != 0))
6606 /* Make sure procedures in interface bodies don't collide. */
6607 gfc_error ("Binding label '%s' in interface body at %L collides "
6608 "with the global entity '%s' at %L",
6610 &(sym->declared_at), bind_c_sym->name,
6611 &(bind_c_sym->where));
6614 else if (sym->attr.contained == 0
6615 && (sym->attr.if_source == IFSRC_UNKNOWN))
6616 if ((sym->attr.use_assoc
6617 && (strcmp (bind_c_sym->mod_name, sym->module) != 0))
6618 || sym->attr.use_assoc == 0)
6620 gfc_error ("Binding label '%s' at %L collides with global "
6621 "entity '%s' at %L", sym->binding_label,
6622 &(sym->declared_at), bind_c_sym->name,
6623 &(bind_c_sym->where));
6628 /* Clear the binding label to prevent checking multiple times. */
6629 sym->binding_label[0] = '\0';
6631 else if (bind_c_sym == NULL)
6633 bind_c_sym = gfc_get_gsymbol (sym->binding_label);
6634 bind_c_sym->where = sym->declared_at;
6635 bind_c_sym->sym_name = sym->name;
6637 if (sym->attr.use_assoc == 1)
6638 bind_c_sym->mod_name = sym->module;
6640 if (sym->ns->proc_name != NULL)
6641 bind_c_sym->mod_name = sym->ns->proc_name->name;
6643 if (sym->attr.contained == 0)
6645 if (sym->attr.subroutine)
6646 bind_c_sym->type = GSYM_SUBROUTINE;
6647 else if (sym->attr.function)
6648 bind_c_sym->type = GSYM_FUNCTION;
6656 /* Resolve an index expression. */
6659 resolve_index_expr (gfc_expr *e)
6661 if (gfc_resolve_expr (e) == FAILURE)
6664 if (gfc_simplify_expr (e, 0) == FAILURE)
6667 if (gfc_specification_expr (e) == FAILURE)
6673 /* Resolve a charlen structure. */
6676 resolve_charlen (gfc_charlen *cl)
6685 specification_expr = 1;
6687 if (resolve_index_expr (cl->length) == FAILURE)
6689 specification_expr = 0;
6693 /* "If the character length parameter value evaluates to a negative
6694 value, the length of character entities declared is zero." */
6695 if (cl->length && !gfc_extract_int (cl->length, &i) && i < 0)
6697 gfc_warning_now ("CHARACTER variable has zero length at %L",
6698 &cl->length->where);
6699 gfc_replace_expr (cl->length, gfc_int_expr (0));
6706 /* Test for non-constant shape arrays. */
6709 is_non_constant_shape_array (gfc_symbol *sym)
6715 not_constant = false;
6716 if (sym->as != NULL)
6718 /* Unfortunately, !gfc_is_compile_time_shape hits a legal case that
6719 has not been simplified; parameter array references. Do the
6720 simplification now. */
6721 for (i = 0; i < sym->as->rank; i++)
6723 e = sym->as->lower[i];
6724 if (e && (resolve_index_expr (e) == FAILURE
6725 || !gfc_is_constant_expr (e)))
6726 not_constant = true;
6728 e = sym->as->upper[i];
6729 if (e && (resolve_index_expr (e) == FAILURE
6730 || !gfc_is_constant_expr (e)))
6731 not_constant = true;
6734 return not_constant;
6737 /* Given a symbol and an initialization expression, add code to initialize
6738 the symbol to the function entry. */
6740 build_init_assign (gfc_symbol *sym, gfc_expr *init)
6744 gfc_namespace *ns = sym->ns;
6746 /* Search for the function namespace if this is a contained
6747 function without an explicit result. */
6748 if (sym->attr.function && sym == sym->result
6749 && sym->name != sym->ns->proc_name->name)
6752 for (;ns; ns = ns->sibling)
6753 if (strcmp (ns->proc_name->name, sym->name) == 0)
6759 gfc_free_expr (init);
6763 /* Build an l-value expression for the result. */
6764 lval = gfc_lval_expr_from_sym (sym);
6766 /* Add the code at scope entry. */
6767 init_st = gfc_get_code ();
6768 init_st->next = ns->code;
6771 /* Assign the default initializer to the l-value. */
6772 init_st->loc = sym->declared_at;
6773 init_st->op = EXEC_INIT_ASSIGN;
6774 init_st->expr = lval;
6775 init_st->expr2 = init;
6778 /* Assign the default initializer to a derived type variable or result. */
6781 apply_default_init (gfc_symbol *sym)
6783 gfc_expr *init = NULL;
6785 if (sym->attr.flavor != FL_VARIABLE && !sym->attr.function)
6788 if (sym->ts.type == BT_DERIVED && sym->ts.derived)
6789 init = gfc_default_initializer (&sym->ts);
6794 build_init_assign (sym, init);
6797 /* Build an initializer for a local integer, real, complex, logical, or
6798 character variable, based on the command line flags finit-local-zero,
6799 finit-integer=, finit-real=, finit-logical=, and finit-runtime. Returns
6800 null if the symbol should not have a default initialization. */
6802 build_default_init_expr (gfc_symbol *sym)
6805 gfc_expr *init_expr;
6809 /* These symbols should never have a default initialization. */
6810 if ((sym->attr.dimension && !gfc_is_compile_time_shape (sym->as))
6811 || sym->attr.external
6813 || sym->attr.pointer
6814 || sym->attr.in_equivalence
6815 || sym->attr.in_common
6818 || sym->attr.cray_pointee
6819 || sym->attr.cray_pointer)
6822 /* Now we'll try to build an initializer expression. */
6823 init_expr = gfc_get_expr ();
6824 init_expr->expr_type = EXPR_CONSTANT;
6825 init_expr->ts.type = sym->ts.type;
6826 init_expr->ts.kind = sym->ts.kind;
6827 init_expr->where = sym->declared_at;
6829 /* We will only initialize integers, reals, complex, logicals, and
6830 characters, and only if the corresponding command-line flags
6831 were set. Otherwise, we free init_expr and return null. */
6832 switch (sym->ts.type)
6835 if (gfc_option.flag_init_integer != GFC_INIT_INTEGER_OFF)
6836 mpz_init_set_si (init_expr->value.integer,
6837 gfc_option.flag_init_integer_value);
6840 gfc_free_expr (init_expr);
6846 mpfr_init (init_expr->value.real);
6847 switch (gfc_option.flag_init_real)
6849 case GFC_INIT_REAL_NAN:
6850 mpfr_set_nan (init_expr->value.real);
6853 case GFC_INIT_REAL_INF:
6854 mpfr_set_inf (init_expr->value.real, 1);
6857 case GFC_INIT_REAL_NEG_INF:
6858 mpfr_set_inf (init_expr->value.real, -1);
6861 case GFC_INIT_REAL_ZERO:
6862 mpfr_set_ui (init_expr->value.real, 0.0, GFC_RND_MODE);
6866 gfc_free_expr (init_expr);
6873 mpfr_init (init_expr->value.complex.r);
6874 mpfr_init (init_expr->value.complex.i);
6875 switch (gfc_option.flag_init_real)
6877 case GFC_INIT_REAL_NAN:
6878 mpfr_set_nan (init_expr->value.complex.r);
6879 mpfr_set_nan (init_expr->value.complex.i);
6882 case GFC_INIT_REAL_INF:
6883 mpfr_set_inf (init_expr->value.complex.r, 1);
6884 mpfr_set_inf (init_expr->value.complex.i, 1);
6887 case GFC_INIT_REAL_NEG_INF:
6888 mpfr_set_inf (init_expr->value.complex.r, -1);
6889 mpfr_set_inf (init_expr->value.complex.i, -1);
6892 case GFC_INIT_REAL_ZERO:
6893 mpfr_set_ui (init_expr->value.complex.r, 0.0, GFC_RND_MODE);
6894 mpfr_set_ui (init_expr->value.complex.i, 0.0, GFC_RND_MODE);
6898 gfc_free_expr (init_expr);
6905 if (gfc_option.flag_init_logical == GFC_INIT_LOGICAL_FALSE)
6906 init_expr->value.logical = 0;
6907 else if (gfc_option.flag_init_logical == GFC_INIT_LOGICAL_TRUE)
6908 init_expr->value.logical = 1;
6911 gfc_free_expr (init_expr);
6917 /* For characters, the length must be constant in order to
6918 create a default initializer. */
6919 if (gfc_option.flag_init_character == GFC_INIT_CHARACTER_ON
6920 && sym->ts.cl->length
6921 && sym->ts.cl->length->expr_type == EXPR_CONSTANT)
6923 char_len = mpz_get_si (sym->ts.cl->length->value.integer);
6924 init_expr->value.character.length = char_len;
6925 init_expr->value.character.string = gfc_getmem (char_len+1);
6926 ch = init_expr->value.character.string;
6927 for (i = 0; i < char_len; i++)
6928 *(ch++) = gfc_option.flag_init_character_value;
6932 gfc_free_expr (init_expr);
6938 gfc_free_expr (init_expr);
6944 /* Add an initialization expression to a local variable. */
6946 apply_default_init_local (gfc_symbol *sym)
6948 gfc_expr *init = NULL;
6950 /* The symbol should be a variable or a function return value. */
6951 if ((sym->attr.flavor != FL_VARIABLE && !sym->attr.function)
6952 || (sym->attr.function && sym->result != sym))
6955 /* Try to build the initializer expression. If we can't initialize
6956 this symbol, then init will be NULL. */
6957 init = build_default_init_expr (sym);
6961 /* For saved variables, we don't want to add an initializer at
6962 function entry, so we just add a static initializer. */
6963 if (sym->attr.save || sym->ns->save_all)
6965 /* Don't clobber an existing initializer! */
6966 gcc_assert (sym->value == NULL);
6971 build_init_assign (sym, init);
6974 /* Resolution of common features of flavors variable and procedure. */
6977 resolve_fl_var_and_proc (gfc_symbol *sym, int mp_flag)
6979 /* Constraints on deferred shape variable. */
6980 if (sym->as == NULL || sym->as->type != AS_DEFERRED)
6982 if (sym->attr.allocatable)
6984 if (sym->attr.dimension)
6985 gfc_error ("Allocatable array '%s' at %L must have "
6986 "a deferred shape", sym->name, &sym->declared_at);
6988 gfc_error ("Scalar object '%s' at %L may not be ALLOCATABLE",
6989 sym->name, &sym->declared_at);
6993 if (sym->attr.pointer && sym->attr.dimension)
6995 gfc_error ("Array pointer '%s' at %L must have a deferred shape",
6996 sym->name, &sym->declared_at);
7003 if (!mp_flag && !sym->attr.allocatable
7004 && !sym->attr.pointer && !sym->attr.dummy)
7006 gfc_error ("Array '%s' at %L cannot have a deferred shape",
7007 sym->name, &sym->declared_at);
7015 /* Additional checks for symbols with flavor variable and derived
7016 type. To be called from resolve_fl_variable. */
7019 resolve_fl_variable_derived (gfc_symbol *sym, int no_init_flag)
7021 gcc_assert (sym->ts.type == BT_DERIVED);
7023 /* Check to see if a derived type is blocked from being host
7024 associated by the presence of another class I symbol in the same
7025 namespace. 14.6.1.3 of the standard and the discussion on
7026 comp.lang.fortran. */
7027 if (sym->ns != sym->ts.derived->ns
7028 && sym->ns->proc_name->attr.if_source != IFSRC_IFBODY)
7031 gfc_find_symbol (sym->ts.derived->name, sym->ns, 0, &s);
7032 if (s && (s->attr.flavor != FL_DERIVED
7033 || !gfc_compare_derived_types (s, sym->ts.derived)))
7035 gfc_error ("The type '%s' cannot be host associated at %L "
7036 "because it is blocked by an incompatible object "
7037 "of the same name declared at %L",
7038 sym->ts.derived->name, &sym->declared_at,
7044 /* 4th constraint in section 11.3: "If an object of a type for which
7045 component-initialization is specified (R429) appears in the
7046 specification-part of a module and does not have the ALLOCATABLE
7047 or POINTER attribute, the object shall have the SAVE attribute."
7049 The check for initializers is performed with
7050 has_default_initializer because gfc_default_initializer generates
7051 a hidden default for allocatable components. */
7052 if (!(sym->value || no_init_flag) && sym->ns->proc_name
7053 && sym->ns->proc_name->attr.flavor == FL_MODULE
7054 && !sym->ns->save_all && !sym->attr.save
7055 && !sym->attr.pointer && !sym->attr.allocatable
7056 && has_default_initializer (sym->ts.derived))
7058 gfc_error("Object '%s' at %L must have the SAVE attribute for "
7059 "default initialization of a component",
7060 sym->name, &sym->declared_at);
7064 /* Assign default initializer. */
7065 if (!(sym->value || sym->attr.pointer || sym->attr.allocatable)
7066 && (!no_init_flag || sym->attr.intent == INTENT_OUT))
7068 sym->value = gfc_default_initializer (&sym->ts);
7075 /* Resolve symbols with flavor variable. */
7078 resolve_fl_variable (gfc_symbol *sym, int mp_flag)
7080 int no_init_flag, automatic_flag;
7082 const char *auto_save_msg;
7084 auto_save_msg = "Automatic object '%s' at %L cannot have the "
7087 if (resolve_fl_var_and_proc (sym, mp_flag) == FAILURE)
7090 /* Set this flag to check that variables are parameters of all entries.
7091 This check is effected by the call to gfc_resolve_expr through
7092 is_non_constant_shape_array. */
7093 specification_expr = 1;
7095 if (sym->ns->proc_name
7096 && (sym->ns->proc_name->attr.flavor == FL_MODULE
7097 || sym->ns->proc_name->attr.is_main_program)
7098 && !sym->attr.use_assoc
7099 && !sym->attr.allocatable
7100 && !sym->attr.pointer
7101 && is_non_constant_shape_array (sym))
7103 /* The shape of a main program or module array needs to be
7105 gfc_error ("The module or main program array '%s' at %L must "
7106 "have constant shape", sym->name, &sym->declared_at);
7107 specification_expr = 0;
7111 if (sym->ts.type == BT_CHARACTER)
7113 /* Make sure that character string variables with assumed length are
7115 e = sym->ts.cl->length;
7116 if (e == NULL && !sym->attr.dummy && !sym->attr.result)
7118 gfc_error ("Entity with assumed character length at %L must be a "
7119 "dummy argument or a PARAMETER", &sym->declared_at);
7123 if (e && sym->attr.save && !gfc_is_constant_expr (e))
7125 gfc_error (auto_save_msg, sym->name, &sym->declared_at);
7129 if (!gfc_is_constant_expr (e)
7130 && !(e->expr_type == EXPR_VARIABLE
7131 && e->symtree->n.sym->attr.flavor == FL_PARAMETER)
7132 && sym->ns->proc_name
7133 && (sym->ns->proc_name->attr.flavor == FL_MODULE
7134 || sym->ns->proc_name->attr.is_main_program)
7135 && !sym->attr.use_assoc)
7137 gfc_error ("'%s' at %L must have constant character length "
7138 "in this context", sym->name, &sym->declared_at);
7143 if (sym->value == NULL && sym->attr.referenced)
7144 apply_default_init_local (sym); /* Try to apply a default initialization. */
7146 /* Determine if the symbol may not have an initializer. */
7147 no_init_flag = automatic_flag = 0;
7148 if (sym->attr.allocatable || sym->attr.external || sym->attr.dummy
7149 || sym->attr.intrinsic || sym->attr.result)
7151 else if (sym->attr.dimension && !sym->attr.pointer
7152 && is_non_constant_shape_array (sym))
7154 no_init_flag = automatic_flag = 1;
7156 /* Also, they must not have the SAVE attribute.
7157 SAVE_IMPLICIT is checked below. */
7158 if (sym->attr.save == SAVE_EXPLICIT)
7160 gfc_error (auto_save_msg, sym->name, &sym->declared_at);
7165 /* Reject illegal initializers. */
7166 if (!sym->mark && sym->value)
7168 if (sym->attr.allocatable)
7169 gfc_error ("Allocatable '%s' at %L cannot have an initializer",
7170 sym->name, &sym->declared_at);
7171 else if (sym->attr.external)
7172 gfc_error ("External '%s' at %L cannot have an initializer",
7173 sym->name, &sym->declared_at);
7174 else if (sym->attr.dummy
7175 && !(sym->ts.type == BT_DERIVED && sym->attr.intent == INTENT_OUT))
7176 gfc_error ("Dummy '%s' at %L cannot have an initializer",
7177 sym->name, &sym->declared_at);
7178 else if (sym->attr.intrinsic)
7179 gfc_error ("Intrinsic '%s' at %L cannot have an initializer",
7180 sym->name, &sym->declared_at);
7181 else if (sym->attr.result)
7182 gfc_error ("Function result '%s' at %L cannot have an initializer",
7183 sym->name, &sym->declared_at);
7184 else if (automatic_flag)
7185 gfc_error ("Automatic array '%s' at %L cannot have an initializer",
7186 sym->name, &sym->declared_at);
7193 if (sym->ts.type == BT_DERIVED)
7194 return resolve_fl_variable_derived (sym, no_init_flag);
7200 /* Resolve a procedure. */
7203 resolve_fl_procedure (gfc_symbol *sym, int mp_flag)
7205 gfc_formal_arglist *arg;
7207 if (sym->attr.ambiguous_interfaces && !sym->attr.referenced)
7208 gfc_warning ("Although not referenced, '%s' at %L has ambiguous "
7209 "interfaces", sym->name, &sym->declared_at);
7211 if (sym->attr.function
7212 && resolve_fl_var_and_proc (sym, mp_flag) == FAILURE)
7215 if (sym->ts.type == BT_CHARACTER)
7217 gfc_charlen *cl = sym->ts.cl;
7219 if (cl && cl->length && gfc_is_constant_expr (cl->length)
7220 && resolve_charlen (cl) == FAILURE)
7223 if (!cl || !cl->length || cl->length->expr_type != EXPR_CONSTANT)
7225 if (sym->attr.proc == PROC_ST_FUNCTION)
7227 gfc_error ("Character-valued statement function '%s' at %L must "
7228 "have constant length", sym->name, &sym->declared_at);
7232 if (sym->attr.external && sym->formal == NULL
7233 && cl && cl->length && cl->length->expr_type != EXPR_CONSTANT)
7235 gfc_error ("Automatic character length function '%s' at %L must "
7236 "have an explicit interface", sym->name,
7243 /* Ensure that derived type for are not of a private type. Internal
7244 module procedures are excluded by 2.2.3.3 - ie. they are not
7245 externally accessible and can access all the objects accessible in
7247 if (!(sym->ns->parent
7248 && sym->ns->parent->proc_name->attr.flavor == FL_MODULE)
7249 && gfc_check_access(sym->attr.access, sym->ns->default_access))
7251 gfc_interface *iface;
7253 for (arg = sym->formal; arg; arg = arg->next)
7256 && arg->sym->ts.type == BT_DERIVED
7257 && !arg->sym->ts.derived->attr.use_assoc
7258 && !gfc_check_access (arg->sym->ts.derived->attr.access,
7259 arg->sym->ts.derived->ns->default_access)
7260 && gfc_notify_std (GFC_STD_F2003, "Fortran 2003: '%s' is of a "
7261 "PRIVATE type and cannot be a dummy argument"
7262 " of '%s', which is PUBLIC at %L",
7263 arg->sym->name, sym->name, &sym->declared_at)
7266 /* Stop this message from recurring. */
7267 arg->sym->ts.derived->attr.access = ACCESS_PUBLIC;
7272 /* PUBLIC interfaces may expose PRIVATE procedures that take types
7273 PRIVATE to the containing module. */
7274 for (iface = sym->generic; iface; iface = iface->next)
7276 for (arg = iface->sym->formal; arg; arg = arg->next)
7279 && arg->sym->ts.type == BT_DERIVED
7280 && !arg->sym->ts.derived->attr.use_assoc
7281 && !gfc_check_access (arg->sym->ts.derived->attr.access,
7282 arg->sym->ts.derived->ns->default_access)
7283 && gfc_notify_std (GFC_STD_F2003, "Fortran 2003: Procedure "
7284 "'%s' in PUBLIC interface '%s' at %L "
7285 "takes dummy arguments of '%s' which is "
7286 "PRIVATE", iface->sym->name, sym->name,
7287 &iface->sym->declared_at,
7288 gfc_typename (&arg->sym->ts)) == FAILURE)
7290 /* Stop this message from recurring. */
7291 arg->sym->ts.derived->attr.access = ACCESS_PUBLIC;
7297 /* PUBLIC interfaces may expose PRIVATE procedures that take types
7298 PRIVATE to the containing module. */
7299 for (iface = sym->generic; iface; iface = iface->next)
7301 for (arg = iface->sym->formal; arg; arg = arg->next)
7304 && arg->sym->ts.type == BT_DERIVED
7305 && !arg->sym->ts.derived->attr.use_assoc
7306 && !gfc_check_access (arg->sym->ts.derived->attr.access,
7307 arg->sym->ts.derived->ns->default_access)
7308 && gfc_notify_std (GFC_STD_F2003, "Fortran 2003: Procedure "
7309 "'%s' in PUBLIC interface '%s' at %L "
7310 "takes dummy arguments of '%s' which is "
7311 "PRIVATE", iface->sym->name, sym->name,
7312 &iface->sym->declared_at,
7313 gfc_typename (&arg->sym->ts)) == FAILURE)
7315 /* Stop this message from recurring. */
7316 arg->sym->ts.derived->attr.access = ACCESS_PUBLIC;
7323 if (sym->attr.function && sym->value && sym->attr.proc != PROC_ST_FUNCTION)
7325 gfc_error ("Function '%s' at %L cannot have an initializer",
7326 sym->name, &sym->declared_at);
7330 /* An external symbol may not have an initializer because it is taken to be
7332 if (sym->attr.external && sym->value)
7334 gfc_error ("External object '%s' at %L may not have an initializer",
7335 sym->name, &sym->declared_at);
7339 /* An elemental function is required to return a scalar 12.7.1 */
7340 if (sym->attr.elemental && sym->attr.function && sym->as)
7342 gfc_error ("ELEMENTAL function '%s' at %L must have a scalar "
7343 "result", sym->name, &sym->declared_at);
7344 /* Reset so that the error only occurs once. */
7345 sym->attr.elemental = 0;
7349 /* 5.1.1.5 of the Standard: A function name declared with an asterisk
7350 char-len-param shall not be array-valued, pointer-valued, recursive
7351 or pure. ....snip... A character value of * may only be used in the
7352 following ways: (i) Dummy arg of procedure - dummy associates with
7353 actual length; (ii) To declare a named constant; or (iii) External
7354 function - but length must be declared in calling scoping unit. */
7355 if (sym->attr.function
7356 && sym->ts.type == BT_CHARACTER
7357 && sym->ts.cl && sym->ts.cl->length == NULL)
7359 if ((sym->as && sym->as->rank) || (sym->attr.pointer)
7360 || (sym->attr.recursive) || (sym->attr.pure))
7362 if (sym->as && sym->as->rank)
7363 gfc_error ("CHARACTER(*) function '%s' at %L cannot be "
7364 "array-valued", sym->name, &sym->declared_at);
7366 if (sym->attr.pointer)
7367 gfc_error ("CHARACTER(*) function '%s' at %L cannot be "
7368 "pointer-valued", sym->name, &sym->declared_at);
7371 gfc_error ("CHARACTER(*) function '%s' at %L cannot be "
7372 "pure", sym->name, &sym->declared_at);
7374 if (sym->attr.recursive)
7375 gfc_error ("CHARACTER(*) function '%s' at %L cannot be "
7376 "recursive", sym->name, &sym->declared_at);
7381 /* Appendix B.2 of the standard. Contained functions give an
7382 error anyway. Fixed-form is likely to be F77/legacy. */
7383 if (!sym->attr.contained && gfc_current_form != FORM_FIXED)
7384 gfc_notify_std (GFC_STD_F95_OBS, "CHARACTER(*) function "
7385 "'%s' at %L is obsolescent in fortran 95",
7386 sym->name, &sym->declared_at);
7389 if (sym->attr.is_bind_c && sym->attr.is_c_interop != 1)
7391 gfc_formal_arglist *curr_arg;
7392 int has_non_interop_arg = 0;
7394 if (verify_bind_c_sym (sym, &(sym->ts), sym->attr.in_common,
7395 sym->common_block) == FAILURE)
7397 /* Clear these to prevent looking at them again if there was an
7399 sym->attr.is_bind_c = 0;
7400 sym->attr.is_c_interop = 0;
7401 sym->ts.is_c_interop = 0;
7405 /* So far, no errors have been found. */
7406 sym->attr.is_c_interop = 1;
7407 sym->ts.is_c_interop = 1;
7410 curr_arg = sym->formal;
7411 while (curr_arg != NULL)
7413 /* Skip implicitly typed dummy args here. */
7414 if (curr_arg->sym->attr.implicit_type == 0)
7415 if (verify_c_interop_param (curr_arg->sym) == FAILURE)
7416 /* If something is found to fail, record the fact so we
7417 can mark the symbol for the procedure as not being
7418 BIND(C) to try and prevent multiple errors being
7420 has_non_interop_arg = 1;
7422 curr_arg = curr_arg->next;
7425 /* See if any of the arguments were not interoperable and if so, clear
7426 the procedure symbol to prevent duplicate error messages. */
7427 if (has_non_interop_arg != 0)
7429 sym->attr.is_c_interop = 0;
7430 sym->ts.is_c_interop = 0;
7431 sym->attr.is_bind_c = 0;
7439 /* Resolve the components of a derived type. */
7442 resolve_fl_derived (gfc_symbol *sym)
7445 gfc_dt_list * dt_list;
7448 for (c = sym->components; c != NULL; c = c->next)
7450 if (c->ts.type == BT_CHARACTER)
7452 if (c->ts.cl->length == NULL
7453 || (resolve_charlen (c->ts.cl) == FAILURE)
7454 || !gfc_is_constant_expr (c->ts.cl->length))
7456 gfc_error ("Character length of component '%s' needs to "
7457 "be a constant specification expression at %L",
7459 c->ts.cl->length ? &c->ts.cl->length->where : &c->loc);
7464 if (c->ts.type == BT_DERIVED
7465 && sym->component_access != ACCESS_PRIVATE
7466 && gfc_check_access (sym->attr.access, sym->ns->default_access)
7467 && !c->ts.derived->attr.use_assoc
7468 && !gfc_check_access (c->ts.derived->attr.access,
7469 c->ts.derived->ns->default_access))
7471 gfc_error ("The component '%s' is a PRIVATE type and cannot be "
7472 "a component of '%s', which is PUBLIC at %L",
7473 c->name, sym->name, &sym->declared_at);
7477 if (sym->attr.sequence)
7479 if (c->ts.type == BT_DERIVED && c->ts.derived->attr.sequence == 0)
7481 gfc_error ("Component %s of SEQUENCE type declared at %L does "
7482 "not have the SEQUENCE attribute",
7483 c->ts.derived->name, &sym->declared_at);
7488 if (c->ts.type == BT_DERIVED && c->pointer
7489 && c->ts.derived->components == NULL)
7491 gfc_error ("The pointer component '%s' of '%s' at %L is a type "
7492 "that has not been declared", c->name, sym->name,
7497 if (c->pointer || c->allocatable || c->as == NULL)
7500 for (i = 0; i < c->as->rank; i++)
7502 if (c->as->lower[i] == NULL
7503 || !gfc_is_constant_expr (c->as->lower[i])
7504 || (resolve_index_expr (c->as->lower[i]) == FAILURE)
7505 || c->as->upper[i] == NULL
7506 || (resolve_index_expr (c->as->upper[i]) == FAILURE)
7507 || !gfc_is_constant_expr (c->as->upper[i]))
7509 gfc_error ("Component '%s' of '%s' at %L must have "
7510 "constant array bounds",
7511 c->name, sym->name, &c->loc);
7517 /* Add derived type to the derived type list. */
7518 for (dt_list = gfc_derived_types; dt_list; dt_list = dt_list->next)
7519 if (sym == dt_list->derived)
7522 if (dt_list == NULL)
7524 dt_list = gfc_get_dt_list ();
7525 dt_list->next = gfc_derived_types;
7526 dt_list->derived = sym;
7527 gfc_derived_types = dt_list;
7535 resolve_fl_namelist (gfc_symbol *sym)
7540 /* Reject PRIVATE objects in a PUBLIC namelist. */
7541 if (gfc_check_access(sym->attr.access, sym->ns->default_access))
7543 for (nl = sym->namelist; nl; nl = nl->next)
7545 if (!nl->sym->attr.use_assoc
7546 && !(sym->ns->parent == nl->sym->ns)
7547 && !(sym->ns->parent
7548 && sym->ns->parent->parent == nl->sym->ns)
7549 && !gfc_check_access(nl->sym->attr.access,
7550 nl->sym->ns->default_access))
7552 gfc_error ("NAMELIST object '%s' was declared PRIVATE and "
7553 "cannot be member of PUBLIC namelist '%s' at %L",
7554 nl->sym->name, sym->name, &sym->declared_at);
7558 /* Types with private components that came here by USE-association. */
7559 if (nl->sym->ts.type == BT_DERIVED
7560 && derived_inaccessible (nl->sym->ts.derived))
7562 gfc_error ("NAMELIST object '%s' has use-associated PRIVATE "
7563 "components and cannot be member of namelist '%s' at %L",
7564 nl->sym->name, sym->name, &sym->declared_at);
7568 /* Types with private components that are defined in the same module. */
7569 if (nl->sym->ts.type == BT_DERIVED
7570 && !(sym->ns->parent == nl->sym->ts.derived->ns)
7571 && !gfc_check_access (nl->sym->ts.derived->attr.private_comp
7572 ? ACCESS_PRIVATE : ACCESS_UNKNOWN,
7573 nl->sym->ns->default_access))
7575 gfc_error ("NAMELIST object '%s' has PRIVATE components and "
7576 "cannot be a member of PUBLIC namelist '%s' at %L",
7577 nl->sym->name, sym->name, &sym->declared_at);
7583 for (nl = sym->namelist; nl; nl = nl->next)
7585 /* Reject namelist arrays of assumed shape. */
7586 if (nl->sym->as && nl->sym->as->type == AS_ASSUMED_SHAPE
7587 && gfc_notify_std (GFC_STD_F2003, "NAMELIST array object '%s' "
7588 "must not have assumed shape in namelist "
7589 "'%s' at %L", nl->sym->name, sym->name,
7590 &sym->declared_at) == FAILURE)
7593 /* Reject namelist arrays that are not constant shape. */
7594 if (is_non_constant_shape_array (nl->sym))
7596 gfc_error ("NAMELIST array object '%s' must have constant "
7597 "shape in namelist '%s' at %L", nl->sym->name,
7598 sym->name, &sym->declared_at);
7602 /* Namelist objects cannot have allocatable or pointer components. */
7603 if (nl->sym->ts.type != BT_DERIVED)
7606 if (nl->sym->ts.derived->attr.alloc_comp)
7608 gfc_error ("NAMELIST object '%s' in namelist '%s' at %L cannot "
7609 "have ALLOCATABLE components",
7610 nl->sym->name, sym->name, &sym->declared_at);
7614 if (nl->sym->ts.derived->attr.pointer_comp)
7616 gfc_error ("NAMELIST object '%s' in namelist '%s' at %L cannot "
7617 "have POINTER components",
7618 nl->sym->name, sym->name, &sym->declared_at);
7624 /* 14.1.2 A module or internal procedure represent local entities
7625 of the same type as a namelist member and so are not allowed. */
7626 for (nl = sym->namelist; nl; nl = nl->next)
7628 if (nl->sym->ts.kind != 0 && nl->sym->attr.flavor == FL_VARIABLE)
7631 if (nl->sym->attr.function && nl->sym == nl->sym->result)
7632 if ((nl->sym == sym->ns->proc_name)
7634 (sym->ns->parent && nl->sym == sym->ns->parent->proc_name))
7638 if (nl->sym && nl->sym->name)
7639 gfc_find_symbol (nl->sym->name, sym->ns, 1, &nlsym);
7640 if (nlsym && nlsym->attr.flavor == FL_PROCEDURE)
7642 gfc_error ("PROCEDURE attribute conflicts with NAMELIST "
7643 "attribute in '%s' at %L", nlsym->name,
7654 resolve_fl_parameter (gfc_symbol *sym)
7656 /* A parameter array's shape needs to be constant. */
7658 && (sym->as->type == AS_DEFERRED
7659 || is_non_constant_shape_array (sym)))
7661 gfc_error ("Parameter array '%s' at %L cannot be automatic "
7662 "or of deferred shape", sym->name, &sym->declared_at);
7666 /* Make sure a parameter that has been implicitly typed still
7667 matches the implicit type, since PARAMETER statements can precede
7668 IMPLICIT statements. */
7669 if (sym->attr.implicit_type
7670 && !gfc_compare_types (&sym->ts, gfc_get_default_type (sym, sym->ns)))
7672 gfc_error ("Implicitly typed PARAMETER '%s' at %L doesn't match a "
7673 "later IMPLICIT type", sym->name, &sym->declared_at);
7677 /* Make sure the types of derived parameters are consistent. This
7678 type checking is deferred until resolution because the type may
7679 refer to a derived type from the host. */
7680 if (sym->ts.type == BT_DERIVED
7681 && !gfc_compare_types (&sym->ts, &sym->value->ts))
7683 gfc_error ("Incompatible derived type in PARAMETER at %L",
7684 &sym->value->where);
7691 /* Do anything necessary to resolve a symbol. Right now, we just
7692 assume that an otherwise unknown symbol is a variable. This sort
7693 of thing commonly happens for symbols in module. */
7696 resolve_symbol (gfc_symbol *sym)
7698 int check_constant, mp_flag;
7699 gfc_symtree *symtree;
7700 gfc_symtree *this_symtree;
7704 if (sym->attr.flavor == FL_UNKNOWN)
7707 /* If we find that a flavorless symbol is an interface in one of the
7708 parent namespaces, find its symtree in this namespace, free the
7709 symbol and set the symtree to point to the interface symbol. */
7710 for (ns = gfc_current_ns->parent; ns; ns = ns->parent)
7712 symtree = gfc_find_symtree (ns->sym_root, sym->name);
7713 if (symtree && symtree->n.sym->generic)
7715 this_symtree = gfc_find_symtree (gfc_current_ns->sym_root,
7719 gfc_free_symbol (sym);
7720 symtree->n.sym->refs++;
7721 this_symtree->n.sym = symtree->n.sym;
7726 /* Otherwise give it a flavor according to such attributes as
7728 if (sym->attr.external == 0 && sym->attr.intrinsic == 0)
7729 sym->attr.flavor = FL_VARIABLE;
7732 sym->attr.flavor = FL_PROCEDURE;
7733 if (sym->attr.dimension)
7734 sym->attr.function = 1;
7738 if (sym->attr.procedure && sym->ts.interface
7739 && sym->attr.if_source != IFSRC_DECL)
7741 if (sym->ts.interface->attr.procedure)
7742 gfc_error ("Interface '%s', used by procedure '%s' at %L, is declared "
7743 "in a later PROCEDURE statement", sym->ts.interface->name,
7744 sym->name,&sym->declared_at);
7746 /* Get the attributes from the interface (now resolved). */
7747 if (sym->ts.interface->attr.if_source || sym->ts.interface->attr.intrinsic)
7749 sym->ts.type = sym->ts.interface->ts.type;
7750 sym->ts.kind = sym->ts.interface->ts.kind;
7751 sym->attr.function = sym->ts.interface->attr.function;
7752 sym->attr.subroutine = sym->ts.interface->attr.subroutine;
7753 copy_formal_args (sym, sym->ts.interface);
7755 else if (sym->ts.interface->name[0] != '\0')
7757 gfc_error ("Interface '%s' of procedure '%s' at %L must be explicit",
7758 sym->ts.interface->name, sym->name, &sym->declared_at);
7763 if (sym->attr.flavor == FL_DERIVED && resolve_fl_derived (sym) == FAILURE)
7766 /* Symbols that are module procedures with results (functions) have
7767 the types and array specification copied for type checking in
7768 procedures that call them, as well as for saving to a module
7769 file. These symbols can't stand the scrutiny that their results
7771 mp_flag = (sym->result != NULL && sym->result != sym);
7774 /* Make sure that the intrinsic is consistent with its internal
7775 representation. This needs to be done before assigning a default
7776 type to avoid spurious warnings. */
7777 if (sym->attr.flavor != FL_MODULE && sym->attr.intrinsic)
7779 if (gfc_intrinsic_name (sym->name, 0))
7781 if (sym->ts.type != BT_UNKNOWN && gfc_option.warn_surprising)
7782 gfc_warning ("Type specified for intrinsic function '%s' at %L is ignored",
7783 sym->name, &sym->declared_at);
7785 else if (gfc_intrinsic_name (sym->name, 1))
7787 if (sym->ts.type != BT_UNKNOWN)
7789 gfc_error ("Intrinsic subroutine '%s' at %L shall not have a type specifier",
7790 sym->name, &sym->declared_at);
7796 gfc_error ("Intrinsic '%s' at %L does not exist", sym->name, &sym->declared_at);
7801 /* Assign default type to symbols that need one and don't have one. */
7802 if (sym->ts.type == BT_UNKNOWN)
7804 if (sym->attr.flavor == FL_VARIABLE || sym->attr.flavor == FL_PARAMETER)
7805 gfc_set_default_type (sym, 1, NULL);
7807 if (sym->attr.flavor == FL_PROCEDURE && sym->attr.function)
7809 /* The specific case of an external procedure should emit an error
7810 in the case that there is no implicit type. */
7812 gfc_set_default_type (sym, sym->attr.external, NULL);
7815 /* Result may be in another namespace. */
7816 resolve_symbol (sym->result);
7818 sym->ts = sym->result->ts;
7819 sym->as = gfc_copy_array_spec (sym->result->as);
7820 sym->attr.dimension = sym->result->attr.dimension;
7821 sym->attr.pointer = sym->result->attr.pointer;
7822 sym->attr.allocatable = sym->result->attr.allocatable;
7827 /* Assumed size arrays and assumed shape arrays must be dummy
7831 && (sym->as->type == AS_ASSUMED_SIZE
7832 || sym->as->type == AS_ASSUMED_SHAPE)
7833 && sym->attr.dummy == 0)
7835 if (sym->as->type == AS_ASSUMED_SIZE)
7836 gfc_error ("Assumed size array at %L must be a dummy argument",
7839 gfc_error ("Assumed shape array at %L must be a dummy argument",
7844 /* Make sure symbols with known intent or optional are really dummy
7845 variable. Because of ENTRY statement, this has to be deferred
7846 until resolution time. */
7848 if (!sym->attr.dummy
7849 && (sym->attr.optional || sym->attr.intent != INTENT_UNKNOWN))
7851 gfc_error ("Symbol at %L is not a DUMMY variable", &sym->declared_at);
7855 if (sym->attr.value && !sym->attr.dummy)
7857 gfc_error ("'%s' at %L cannot have the VALUE attribute because "
7858 "it is not a dummy argument", sym->name, &sym->declared_at);
7862 if (sym->attr.value && sym->ts.type == BT_CHARACTER)
7864 gfc_charlen *cl = sym->ts.cl;
7865 if (!cl || !cl->length || cl->length->expr_type != EXPR_CONSTANT)
7867 gfc_error ("Character dummy variable '%s' at %L with VALUE "
7868 "attribute must have constant length",
7869 sym->name, &sym->declared_at);
7873 if (sym->ts.is_c_interop
7874 && mpz_cmp_si (cl->length->value.integer, 1) != 0)
7876 gfc_error ("C interoperable character dummy variable '%s' at %L "
7877 "with VALUE attribute must have length one",
7878 sym->name, &sym->declared_at);
7883 /* If the symbol is marked as bind(c), verify it's type and kind. Do not
7884 do this for something that was implicitly typed because that is handled
7885 in gfc_set_default_type. Handle dummy arguments and procedure
7886 definitions separately. Also, anything that is use associated is not
7887 handled here but instead is handled in the module it is declared in.
7888 Finally, derived type definitions are allowed to be BIND(C) since that
7889 only implies that they're interoperable, and they are checked fully for
7890 interoperability when a variable is declared of that type. */
7891 if (sym->attr.is_bind_c && sym->attr.implicit_type == 0 &&
7892 sym->attr.use_assoc == 0 && sym->attr.dummy == 0 &&
7893 sym->attr.flavor != FL_PROCEDURE && sym->attr.flavor != FL_DERIVED)
7897 /* First, make sure the variable is declared at the
7898 module-level scope (J3/04-007, Section 15.3). */
7899 if (sym->ns->proc_name->attr.flavor != FL_MODULE &&
7900 sym->attr.in_common == 0)
7902 gfc_error ("Variable '%s' at %L cannot be BIND(C) because it "
7903 "is neither a COMMON block nor declared at the "
7904 "module level scope", sym->name, &(sym->declared_at));
7907 else if (sym->common_head != NULL)
7909 t = verify_com_block_vars_c_interop (sym->common_head);
7913 /* If type() declaration, we need to verify that the components
7914 of the given type are all C interoperable, etc. */
7915 if (sym->ts.type == BT_DERIVED &&
7916 sym->ts.derived->attr.is_c_interop != 1)
7918 /* Make sure the user marked the derived type as BIND(C). If
7919 not, call the verify routine. This could print an error
7920 for the derived type more than once if multiple variables
7921 of that type are declared. */
7922 if (sym->ts.derived->attr.is_bind_c != 1)
7923 verify_bind_c_derived_type (sym->ts.derived);
7927 /* Verify the variable itself as C interoperable if it
7928 is BIND(C). It is not possible for this to succeed if
7929 the verify_bind_c_derived_type failed, so don't have to handle
7930 any error returned by verify_bind_c_derived_type. */
7931 t = verify_bind_c_sym (sym, &(sym->ts), sym->attr.in_common,
7937 /* clear the is_bind_c flag to prevent reporting errors more than
7938 once if something failed. */
7939 sym->attr.is_bind_c = 0;
7944 /* If a derived type symbol has reached this point, without its
7945 type being declared, we have an error. Notice that most
7946 conditions that produce undefined derived types have already
7947 been dealt with. However, the likes of:
7948 implicit type(t) (t) ..... call foo (t) will get us here if
7949 the type is not declared in the scope of the implicit
7950 statement. Change the type to BT_UNKNOWN, both because it is so
7951 and to prevent an ICE. */
7952 if (sym->ts.type == BT_DERIVED && sym->ts.derived->components == NULL
7953 && !sym->ts.derived->attr.zero_comp)
7955 gfc_error ("The derived type '%s' at %L is of type '%s', "
7956 "which has not been defined", sym->name,
7957 &sym->declared_at, sym->ts.derived->name);
7958 sym->ts.type = BT_UNKNOWN;
7962 /* Make sure that the derived type has been resolved and that the
7963 derived type is visible in the symbol's namespace, if it is a
7964 module function and is not PRIVATE. */
7965 if (sym->ts.type == BT_DERIVED
7966 && sym->ts.derived->attr.use_assoc
7967 && sym->ns->proc_name->attr.flavor == FL_MODULE)
7971 if (resolve_fl_derived (sym->ts.derived) == FAILURE)
7974 gfc_find_symbol (sym->ts.derived->name, sym->ns, 1, &ds);
7975 if (!ds && sym->attr.function
7976 && gfc_check_access (sym->attr.access, sym->ns->default_access))
7978 symtree = gfc_new_symtree (&sym->ns->sym_root,
7979 sym->ts.derived->name);
7980 symtree->n.sym = sym->ts.derived;
7981 sym->ts.derived->refs++;
7985 /* Unless the derived-type declaration is use associated, Fortran 95
7986 does not allow public entries of private derived types.
7987 See 4.4.1 (F95) and 4.5.1.1 (F2003); and related interpretation
7989 if (sym->ts.type == BT_DERIVED
7990 && sym->ns->proc_name && sym->ns->proc_name->attr.flavor == FL_MODULE
7991 && !sym->ts.derived->attr.use_assoc
7992 && gfc_check_access (sym->attr.access, sym->ns->default_access)
7993 && !gfc_check_access (sym->ts.derived->attr.access,
7994 sym->ts.derived->ns->default_access)
7995 && gfc_notify_std (GFC_STD_F2003, "Fortran 2003: PUBLIC %s '%s' at %L "
7996 "of PRIVATE derived type '%s'",
7997 (sym->attr.flavor == FL_PARAMETER) ? "parameter"
7998 : "variable", sym->name, &sym->declared_at,
7999 sym->ts.derived->name) == FAILURE)
8002 /* An assumed-size array with INTENT(OUT) shall not be of a type for which
8003 default initialization is defined (5.1.2.4.4). */
8004 if (sym->ts.type == BT_DERIVED
8006 && sym->attr.intent == INTENT_OUT
8008 && sym->as->type == AS_ASSUMED_SIZE)
8010 for (c = sym->ts.derived->components; c; c = c->next)
8014 gfc_error ("The INTENT(OUT) dummy argument '%s' at %L is "
8015 "ASSUMED SIZE and so cannot have a default initializer",
8016 sym->name, &sym->declared_at);
8022 switch (sym->attr.flavor)
8025 if (resolve_fl_variable (sym, mp_flag) == FAILURE)
8030 if (resolve_fl_procedure (sym, mp_flag) == FAILURE)
8035 if (resolve_fl_namelist (sym) == FAILURE)
8040 if (resolve_fl_parameter (sym) == FAILURE)
8048 /* Resolve array specifier. Check as well some constraints
8049 on COMMON blocks. */
8051 check_constant = sym->attr.in_common && !sym->attr.pointer;
8053 /* Set the formal_arg_flag so that check_conflict will not throw
8054 an error for host associated variables in the specification
8055 expression for an array_valued function. */
8056 if (sym->attr.function && sym->as)
8057 formal_arg_flag = 1;
8059 gfc_resolve_array_spec (sym->as, check_constant);
8061 formal_arg_flag = 0;
8063 /* Resolve formal namespaces. */
8064 if (sym->formal_ns && sym->formal_ns != gfc_current_ns)
8065 gfc_resolve (sym->formal_ns);
8067 /* Check threadprivate restrictions. */
8068 if (sym->attr.threadprivate && !sym->attr.save && !sym->ns->save_all
8069 && (!sym->attr.in_common
8070 && sym->module == NULL
8071 && (sym->ns->proc_name == NULL
8072 || sym->ns->proc_name->attr.flavor != FL_MODULE)))
8073 gfc_error ("Threadprivate at %L isn't SAVEd", &sym->declared_at);
8075 /* If we have come this far we can apply default-initializers, as
8076 described in 14.7.5, to those variables that have not already
8077 been assigned one. */
8078 if (sym->ts.type == BT_DERIVED
8079 && sym->attr.referenced
8080 && sym->ns == gfc_current_ns
8082 && !sym->attr.allocatable
8083 && !sym->attr.alloc_comp)
8085 symbol_attribute *a = &sym->attr;
8087 if ((!a->save && !a->dummy && !a->pointer
8088 && !a->in_common && !a->use_assoc
8089 && !(a->function && sym != sym->result))
8090 || (a->dummy && a->intent == INTENT_OUT))
8091 apply_default_init (sym);
8096 /************* Resolve DATA statements *************/
8100 gfc_data_value *vnode;
8106 /* Advance the values structure to point to the next value in the data list. */
8109 next_data_value (void)
8112 while (mpz_cmp_ui (values.left, 0) == 0)
8114 if (values.vnode->next == NULL)
8117 values.vnode = values.vnode->next;
8118 mpz_set (values.left, values.vnode->repeat);
8126 check_data_variable (gfc_data_variable *var, locus *where)
8132 ar_type mark = AR_UNKNOWN;
8134 mpz_t section_index[GFC_MAX_DIMENSIONS];
8138 if (gfc_resolve_expr (var->expr) == FAILURE)
8142 mpz_init_set_si (offset, 0);
8145 if (e->expr_type != EXPR_VARIABLE)
8146 gfc_internal_error ("check_data_variable(): Bad expression");
8148 if (e->symtree->n.sym->ns->is_block_data
8149 && !e->symtree->n.sym->attr.in_common)
8151 gfc_error ("BLOCK DATA element '%s' at %L must be in COMMON",
8152 e->symtree->n.sym->name, &e->symtree->n.sym->declared_at);
8155 if (e->ref == NULL && e->symtree->n.sym->as)
8157 gfc_error ("DATA array '%s' at %L must be specified in a previous"
8158 " declaration", e->symtree->n.sym->name, where);
8164 mpz_init_set_ui (size, 1);
8171 /* Find the array section reference. */
8172 for (ref = e->ref; ref; ref = ref->next)
8174 if (ref->type != REF_ARRAY)
8176 if (ref->u.ar.type == AR_ELEMENT)
8182 /* Set marks according to the reference pattern. */
8183 switch (ref->u.ar.type)
8191 /* Get the start position of array section. */
8192 gfc_get_section_index (ar, section_index, &offset);
8200 if (gfc_array_size (e, &size) == FAILURE)
8202 gfc_error ("Nonconstant array section at %L in DATA statement",
8211 while (mpz_cmp_ui (size, 0) > 0)
8213 if (next_data_value () == FAILURE)
8215 gfc_error ("DATA statement at %L has more variables than values",
8221 t = gfc_check_assign (var->expr, values.vnode->expr, 0);
8225 /* If we have more than one element left in the repeat count,
8226 and we have more than one element left in the target variable,
8227 then create a range assignment. */
8228 /* FIXME: Only done for full arrays for now, since array sections
8230 if (mark == AR_FULL && ref && ref->next == NULL
8231 && mpz_cmp_ui (values.left, 1) > 0 && mpz_cmp_ui (size, 1) > 0)
8235 if (mpz_cmp (size, values.left) >= 0)
8237 mpz_init_set (range, values.left);
8238 mpz_sub (size, size, values.left);
8239 mpz_set_ui (values.left, 0);
8243 mpz_init_set (range, size);
8244 mpz_sub (values.left, values.left, size);
8245 mpz_set_ui (size, 0);
8248 gfc_assign_data_value_range (var->expr, values.vnode->expr,
8251 mpz_add (offset, offset, range);
8255 /* Assign initial value to symbol. */
8258 mpz_sub_ui (values.left, values.left, 1);
8259 mpz_sub_ui (size, size, 1);
8261 t = gfc_assign_data_value (var->expr, values.vnode->expr, offset);
8265 if (mark == AR_FULL)
8266 mpz_add_ui (offset, offset, 1);
8268 /* Modify the array section indexes and recalculate the offset
8269 for next element. */
8270 else if (mark == AR_SECTION)
8271 gfc_advance_section (section_index, ar, &offset);
8275 if (mark == AR_SECTION)
8277 for (i = 0; i < ar->dimen; i++)
8278 mpz_clear (section_index[i]);
8288 static try traverse_data_var (gfc_data_variable *, locus *);
8290 /* Iterate over a list of elements in a DATA statement. */
8293 traverse_data_list (gfc_data_variable *var, locus *where)
8296 iterator_stack frame;
8297 gfc_expr *e, *start, *end, *step;
8298 try retval = SUCCESS;
8300 mpz_init (frame.value);
8302 start = gfc_copy_expr (var->iter.start);
8303 end = gfc_copy_expr (var->iter.end);
8304 step = gfc_copy_expr (var->iter.step);
8306 if (gfc_simplify_expr (start, 1) == FAILURE
8307 || start->expr_type != EXPR_CONSTANT)
8309 gfc_error ("iterator start at %L does not simplify", &start->where);
8313 if (gfc_simplify_expr (end, 1) == FAILURE
8314 || end->expr_type != EXPR_CONSTANT)
8316 gfc_error ("iterator end at %L does not simplify", &end->where);
8320 if (gfc_simplify_expr (step, 1) == FAILURE
8321 || step->expr_type != EXPR_CONSTANT)
8323 gfc_error ("iterator step at %L does not simplify", &step->where);
8328 mpz_init_set (trip, end->value.integer);
8329 mpz_sub (trip, trip, start->value.integer);
8330 mpz_add (trip, trip, step->value.integer);
8332 mpz_div (trip, trip, step->value.integer);
8334 mpz_set (frame.value, start->value.integer);
8336 frame.prev = iter_stack;
8337 frame.variable = var->iter.var->symtree;
8338 iter_stack = &frame;
8340 while (mpz_cmp_ui (trip, 0) > 0)
8342 if (traverse_data_var (var->list, where) == FAILURE)
8349 e = gfc_copy_expr (var->expr);
8350 if (gfc_simplify_expr (e, 1) == FAILURE)
8358 mpz_add (frame.value, frame.value, step->value.integer);
8360 mpz_sub_ui (trip, trip, 1);
8365 mpz_clear (frame.value);
8367 gfc_free_expr (start);
8368 gfc_free_expr (end);
8369 gfc_free_expr (step);
8371 iter_stack = frame.prev;
8376 /* Type resolve variables in the variable list of a DATA statement. */
8379 traverse_data_var (gfc_data_variable *var, locus *where)
8383 for (; var; var = var->next)
8385 if (var->expr == NULL)
8386 t = traverse_data_list (var, where);
8388 t = check_data_variable (var, where);
8398 /* Resolve the expressions and iterators associated with a data statement.
8399 This is separate from the assignment checking because data lists should
8400 only be resolved once. */
8403 resolve_data_variables (gfc_data_variable *d)
8405 for (; d; d = d->next)
8407 if (d->list == NULL)
8409 if (gfc_resolve_expr (d->expr) == FAILURE)
8414 if (gfc_resolve_iterator (&d->iter, false) == FAILURE)
8417 if (resolve_data_variables (d->list) == FAILURE)
8426 /* Resolve a single DATA statement. We implement this by storing a pointer to
8427 the value list into static variables, and then recursively traversing the
8428 variables list, expanding iterators and such. */
8431 resolve_data (gfc_data *d)
8434 if (resolve_data_variables (d->var) == FAILURE)
8437 values.vnode = d->value;
8438 if (d->value == NULL)
8439 mpz_set_ui (values.left, 0);
8441 mpz_set (values.left, d->value->repeat);
8443 if (traverse_data_var (d->var, &d->where) == FAILURE)
8446 /* At this point, we better not have any values left. */
8448 if (next_data_value () == SUCCESS)
8449 gfc_error ("DATA statement at %L has more values than variables",
8454 /* 12.6 Constraint: In a pure subprogram any variable which is in common or
8455 accessed by host or use association, is a dummy argument to a pure function,
8456 is a dummy argument with INTENT (IN) to a pure subroutine, or an object that
8457 is storage associated with any such variable, shall not be used in the
8458 following contexts: (clients of this function). */
8460 /* Determines if a variable is not 'pure', ie not assignable within a pure
8461 procedure. Returns zero if assignment is OK, nonzero if there is a
8464 gfc_impure_variable (gfc_symbol *sym)
8468 if (sym->attr.use_assoc || sym->attr.in_common)
8471 if (sym->ns != gfc_current_ns)
8472 return !sym->attr.function;
8474 proc = sym->ns->proc_name;
8475 if (sym->attr.dummy && gfc_pure (proc)
8476 && ((proc->attr.subroutine && sym->attr.intent == INTENT_IN)
8478 proc->attr.function))
8481 /* TODO: Sort out what can be storage associated, if anything, and include
8482 it here. In principle equivalences should be scanned but it does not
8483 seem to be possible to storage associate an impure variable this way. */
8488 /* Test whether a symbol is pure or not. For a NULL pointer, checks the
8489 symbol of the current procedure. */
8492 gfc_pure (gfc_symbol *sym)
8494 symbol_attribute attr;
8497 sym = gfc_current_ns->proc_name;
8503 return attr.flavor == FL_PROCEDURE && (attr.pure || attr.elemental);
8507 /* Test whether the current procedure is elemental or not. */
8510 gfc_elemental (gfc_symbol *sym)
8512 symbol_attribute attr;
8515 sym = gfc_current_ns->proc_name;
8520 return attr.flavor == FL_PROCEDURE && attr.elemental;
8524 /* Warn about unused labels. */
8527 warn_unused_fortran_label (gfc_st_label *label)
8532 warn_unused_fortran_label (label->left);
8534 if (label->defined == ST_LABEL_UNKNOWN)
8537 switch (label->referenced)
8539 case ST_LABEL_UNKNOWN:
8540 gfc_warning ("Label %d at %L defined but not used", label->value,
8544 case ST_LABEL_BAD_TARGET:
8545 gfc_warning ("Label %d at %L defined but cannot be used",
8546 label->value, &label->where);
8553 warn_unused_fortran_label (label->right);
8557 /* Returns the sequence type of a symbol or sequence. */
8560 sequence_type (gfc_typespec ts)
8569 if (ts.derived->components == NULL)
8570 return SEQ_NONDEFAULT;
8572 result = sequence_type (ts.derived->components->ts);
8573 for (c = ts.derived->components->next; c; c = c->next)
8574 if (sequence_type (c->ts) != result)
8580 if (ts.kind != gfc_default_character_kind)
8581 return SEQ_NONDEFAULT;
8583 return SEQ_CHARACTER;
8586 if (ts.kind != gfc_default_integer_kind)
8587 return SEQ_NONDEFAULT;
8592 if (!(ts.kind == gfc_default_real_kind
8593 || ts.kind == gfc_default_double_kind))
8594 return SEQ_NONDEFAULT;
8599 if (ts.kind != gfc_default_complex_kind)
8600 return SEQ_NONDEFAULT;
8605 if (ts.kind != gfc_default_logical_kind)
8606 return SEQ_NONDEFAULT;
8611 return SEQ_NONDEFAULT;
8616 /* Resolve derived type EQUIVALENCE object. */
8619 resolve_equivalence_derived (gfc_symbol *derived, gfc_symbol *sym, gfc_expr *e)
8622 gfc_component *c = derived->components;
8627 /* Shall not be an object of nonsequence derived type. */
8628 if (!derived->attr.sequence)
8630 gfc_error ("Derived type variable '%s' at %L must have SEQUENCE "
8631 "attribute to be an EQUIVALENCE object", sym->name,
8636 /* Shall not have allocatable components. */
8637 if (derived->attr.alloc_comp)
8639 gfc_error ("Derived type variable '%s' at %L cannot have ALLOCATABLE "
8640 "components to be an EQUIVALENCE object",sym->name,
8645 if (sym->attr.in_common && has_default_initializer (sym->ts.derived))
8647 gfc_error ("Derived type variable '%s' at %L with default "
8648 "initialization cannot be in EQUIVALENCE with a variable "
8649 "in COMMON", sym->name, &e->where);
8653 for (; c ; c = c->next)
8657 && (resolve_equivalence_derived (c->ts.derived, sym, e) == FAILURE))
8660 /* Shall not be an object of sequence derived type containing a pointer
8661 in the structure. */
8664 gfc_error ("Derived type variable '%s' at %L with pointer "
8665 "component(s) cannot be an EQUIVALENCE object",
8666 sym->name, &e->where);
8674 /* Resolve equivalence object.
8675 An EQUIVALENCE object shall not be a dummy argument, a pointer, a target,
8676 an allocatable array, an object of nonsequence derived type, an object of
8677 sequence derived type containing a pointer at any level of component
8678 selection, an automatic object, a function name, an entry name, a result
8679 name, a named constant, a structure component, or a subobject of any of
8680 the preceding objects. A substring shall not have length zero. A
8681 derived type shall not have components with default initialization nor
8682 shall two objects of an equivalence group be initialized.
8683 Either all or none of the objects shall have an protected attribute.
8684 The simple constraints are done in symbol.c(check_conflict) and the rest
8685 are implemented here. */
8688 resolve_equivalence (gfc_equiv *eq)
8691 gfc_symbol *derived;
8692 gfc_symbol *first_sym;
8695 locus *last_where = NULL;
8696 seq_type eq_type, last_eq_type;
8697 gfc_typespec *last_ts;
8698 int object, cnt_protected;
8699 const char *value_name;
8703 last_ts = &eq->expr->symtree->n.sym->ts;
8705 first_sym = eq->expr->symtree->n.sym;
8709 for (object = 1; eq; eq = eq->eq, object++)
8713 e->ts = e->symtree->n.sym->ts;
8714 /* match_varspec might not know yet if it is seeing
8715 array reference or substring reference, as it doesn't
8717 if (e->ref && e->ref->type == REF_ARRAY)
8719 gfc_ref *ref = e->ref;
8720 sym = e->symtree->n.sym;
8722 if (sym->attr.dimension)
8724 ref->u.ar.as = sym->as;
8728 /* For substrings, convert REF_ARRAY into REF_SUBSTRING. */
8729 if (e->ts.type == BT_CHARACTER
8731 && ref->type == REF_ARRAY
8732 && ref->u.ar.dimen == 1
8733 && ref->u.ar.dimen_type[0] == DIMEN_RANGE
8734 && ref->u.ar.stride[0] == NULL)
8736 gfc_expr *start = ref->u.ar.start[0];
8737 gfc_expr *end = ref->u.ar.end[0];
8740 /* Optimize away the (:) reference. */
8741 if (start == NULL && end == NULL)
8746 e->ref->next = ref->next;
8751 ref->type = REF_SUBSTRING;
8753 start = gfc_int_expr (1);
8754 ref->u.ss.start = start;
8755 if (end == NULL && e->ts.cl)
8756 end = gfc_copy_expr (e->ts.cl->length);
8757 ref->u.ss.end = end;
8758 ref->u.ss.length = e->ts.cl;
8765 /* Any further ref is an error. */
8768 gcc_assert (ref->type == REF_ARRAY);
8769 gfc_error ("Syntax error in EQUIVALENCE statement at %L",
8775 if (gfc_resolve_expr (e) == FAILURE)
8778 sym = e->symtree->n.sym;
8780 if (sym->attr.protected)
8782 if (cnt_protected > 0 && cnt_protected != object)
8784 gfc_error ("Either all or none of the objects in the "
8785 "EQUIVALENCE set at %L shall have the "
8786 "PROTECTED attribute",
8791 /* Shall not equivalence common block variables in a PURE procedure. */
8792 if (sym->ns->proc_name
8793 && sym->ns->proc_name->attr.pure
8794 && sym->attr.in_common)
8796 gfc_error ("Common block member '%s' at %L cannot be an EQUIVALENCE "
8797 "object in the pure procedure '%s'",
8798 sym->name, &e->where, sym->ns->proc_name->name);
8802 /* Shall not be a named constant. */
8803 if (e->expr_type == EXPR_CONSTANT)
8805 gfc_error ("Named constant '%s' at %L cannot be an EQUIVALENCE "
8806 "object", sym->name, &e->where);
8810 derived = e->ts.derived;
8811 if (derived && resolve_equivalence_derived (derived, sym, e) == FAILURE)
8814 /* Check that the types correspond correctly:
8816 A numeric sequence structure may be equivalenced to another sequence
8817 structure, an object of default integer type, default real type, double
8818 precision real type, default logical type such that components of the
8819 structure ultimately only become associated to objects of the same
8820 kind. A character sequence structure may be equivalenced to an object
8821 of default character kind or another character sequence structure.
8822 Other objects may be equivalenced only to objects of the same type and
8825 /* Identical types are unconditionally OK. */
8826 if (object == 1 || gfc_compare_types (last_ts, &sym->ts))
8827 goto identical_types;
8829 last_eq_type = sequence_type (*last_ts);
8830 eq_type = sequence_type (sym->ts);
8832 /* Since the pair of objects is not of the same type, mixed or
8833 non-default sequences can be rejected. */
8835 msg = "Sequence %s with mixed components in EQUIVALENCE "
8836 "statement at %L with different type objects";
8838 && last_eq_type == SEQ_MIXED
8839 && gfc_notify_std (GFC_STD_GNU, msg, first_sym->name, last_where)
8841 || (eq_type == SEQ_MIXED
8842 && gfc_notify_std (GFC_STD_GNU, msg, sym->name,
8843 &e->where) == FAILURE))
8846 msg = "Non-default type object or sequence %s in EQUIVALENCE "
8847 "statement at %L with objects of different type";
8849 && last_eq_type == SEQ_NONDEFAULT
8850 && gfc_notify_std (GFC_STD_GNU, msg, first_sym->name,
8851 last_where) == FAILURE)
8852 || (eq_type == SEQ_NONDEFAULT
8853 && gfc_notify_std (GFC_STD_GNU, msg, sym->name,
8854 &e->where) == FAILURE))
8857 msg ="Non-CHARACTER object '%s' in default CHARACTER "
8858 "EQUIVALENCE statement at %L";
8859 if (last_eq_type == SEQ_CHARACTER
8860 && eq_type != SEQ_CHARACTER
8861 && gfc_notify_std (GFC_STD_GNU, msg, sym->name,
8862 &e->where) == FAILURE)
8865 msg ="Non-NUMERIC object '%s' in default NUMERIC "
8866 "EQUIVALENCE statement at %L";
8867 if (last_eq_type == SEQ_NUMERIC
8868 && eq_type != SEQ_NUMERIC
8869 && gfc_notify_std (GFC_STD_GNU, msg, sym->name,
8870 &e->where) == FAILURE)
8875 last_where = &e->where;
8880 /* Shall not be an automatic array. */
8881 if (e->ref->type == REF_ARRAY
8882 && gfc_resolve_array_spec (e->ref->u.ar.as, 1) == FAILURE)
8884 gfc_error ("Array '%s' at %L with non-constant bounds cannot be "
8885 "an EQUIVALENCE object", sym->name, &e->where);
8892 /* Shall not be a structure component. */
8893 if (r->type == REF_COMPONENT)
8895 gfc_error ("Structure component '%s' at %L cannot be an "
8896 "EQUIVALENCE object",
8897 r->u.c.component->name, &e->where);
8901 /* A substring shall not have length zero. */
8902 if (r->type == REF_SUBSTRING)
8904 if (compare_bound (r->u.ss.start, r->u.ss.end) == CMP_GT)
8906 gfc_error ("Substring at %L has length zero",
8907 &r->u.ss.start->where);
8917 /* Resolve function and ENTRY types, issue diagnostics if needed. */
8920 resolve_fntype (gfc_namespace *ns)
8925 if (ns->proc_name == NULL || !ns->proc_name->attr.function)
8928 /* If there are any entries, ns->proc_name is the entry master
8929 synthetic symbol and ns->entries->sym actual FUNCTION symbol. */
8931 sym = ns->entries->sym;
8933 sym = ns->proc_name;
8934 if (sym->result == sym
8935 && sym->ts.type == BT_UNKNOWN
8936 && gfc_set_default_type (sym, 0, NULL) == FAILURE
8937 && !sym->attr.untyped)
8939 gfc_error ("Function '%s' at %L has no IMPLICIT type",
8940 sym->name, &sym->declared_at);
8941 sym->attr.untyped = 1;
8944 if (sym->ts.type == BT_DERIVED && !sym->ts.derived->attr.use_assoc
8945 && !gfc_check_access (sym->ts.derived->attr.access,
8946 sym->ts.derived->ns->default_access)
8947 && gfc_check_access (sym->attr.access, sym->ns->default_access))
8949 gfc_error ("PUBLIC function '%s' at %L cannot be of PRIVATE type '%s'",
8950 sym->name, &sym->declared_at, sym->ts.derived->name);
8954 for (el = ns->entries->next; el; el = el->next)
8956 if (el->sym->result == el->sym
8957 && el->sym->ts.type == BT_UNKNOWN
8958 && gfc_set_default_type (el->sym, 0, NULL) == FAILURE
8959 && !el->sym->attr.untyped)
8961 gfc_error ("ENTRY '%s' at %L has no IMPLICIT type",
8962 el->sym->name, &el->sym->declared_at);
8963 el->sym->attr.untyped = 1;
8968 /* 12.3.2.1.1 Defined operators. */
8971 gfc_resolve_uops (gfc_symtree *symtree)
8975 gfc_formal_arglist *formal;
8977 if (symtree == NULL)
8980 gfc_resolve_uops (symtree->left);
8981 gfc_resolve_uops (symtree->right);
8983 for (itr = symtree->n.uop->operator; itr; itr = itr->next)
8986 if (!sym->attr.function)
8987 gfc_error ("User operator procedure '%s' at %L must be a FUNCTION",
8988 sym->name, &sym->declared_at);
8990 if (sym->ts.type == BT_CHARACTER
8991 && !(sym->ts.cl && sym->ts.cl->length)
8992 && !(sym->result && sym->result->ts.cl
8993 && sym->result->ts.cl->length))
8994 gfc_error ("User operator procedure '%s' at %L cannot be assumed "
8995 "character length", sym->name, &sym->declared_at);
8997 formal = sym->formal;
8998 if (!formal || !formal->sym)
9000 gfc_error ("User operator procedure '%s' at %L must have at least "
9001 "one argument", sym->name, &sym->declared_at);
9005 if (formal->sym->attr.intent != INTENT_IN)
9006 gfc_error ("First argument of operator interface at %L must be "
9007 "INTENT(IN)", &sym->declared_at);
9009 if (formal->sym->attr.optional)
9010 gfc_error ("First argument of operator interface at %L cannot be "
9011 "optional", &sym->declared_at);
9013 formal = formal->next;
9014 if (!formal || !formal->sym)
9017 if (formal->sym->attr.intent != INTENT_IN)
9018 gfc_error ("Second argument of operator interface at %L must be "
9019 "INTENT(IN)", &sym->declared_at);
9021 if (formal->sym->attr.optional)
9022 gfc_error ("Second argument of operator interface at %L cannot be "
9023 "optional", &sym->declared_at);
9026 gfc_error ("Operator interface at %L must have, at most, two "
9027 "arguments", &sym->declared_at);
9032 /* Examine all of the expressions associated with a program unit,
9033 assign types to all intermediate expressions, make sure that all
9034 assignments are to compatible types and figure out which names
9035 refer to which functions or subroutines. It doesn't check code
9036 block, which is handled by resolve_code. */
9039 resolve_types (gfc_namespace *ns)
9046 gfc_current_ns = ns;
9048 resolve_entries (ns);
9050 resolve_common_vars (ns->blank_common.head, false);
9051 resolve_common_blocks (ns->common_root);
9053 resolve_contained_functions (ns);
9055 gfc_traverse_ns (ns, resolve_bind_c_derived_types);
9057 for (cl = ns->cl_list; cl; cl = cl->next)
9058 resolve_charlen (cl);
9060 gfc_traverse_ns (ns, resolve_symbol);
9062 resolve_fntype (ns);
9064 for (n = ns->contained; n; n = n->sibling)
9066 if (gfc_pure (ns->proc_name) && !gfc_pure (n->proc_name))
9067 gfc_error ("Contained procedure '%s' at %L of a PURE procedure must "
9068 "also be PURE", n->proc_name->name,
9069 &n->proc_name->declared_at);
9075 gfc_check_interfaces (ns);
9077 gfc_traverse_ns (ns, resolve_values);
9083 for (d = ns->data; d; d = d->next)
9087 gfc_traverse_ns (ns, gfc_formalize_init_value);
9089 gfc_traverse_ns (ns, gfc_verify_binding_labels);
9091 if (ns->common_root != NULL)
9092 gfc_traverse_symtree (ns->common_root, resolve_bind_c_comms);
9094 for (eq = ns->equiv; eq; eq = eq->next)
9095 resolve_equivalence (eq);
9097 /* Warn about unused labels. */
9098 if (warn_unused_label)
9099 warn_unused_fortran_label (ns->st_labels);
9101 gfc_resolve_uops (ns->uop_root);
9105 /* Call resolve_code recursively. */
9108 resolve_codes (gfc_namespace *ns)
9112 for (n = ns->contained; n; n = n->sibling)
9115 gfc_current_ns = ns;
9117 /* Set to an out of range value. */
9118 current_entry_id = -1;
9120 bitmap_obstack_initialize (&labels_obstack);
9121 resolve_code (ns->code, ns);
9122 bitmap_obstack_release (&labels_obstack);
9126 /* This function is called after a complete program unit has been compiled.
9127 Its purpose is to examine all of the expressions associated with a program
9128 unit, assign types to all intermediate expressions, make sure that all
9129 assignments are to compatible types and figure out which names refer to
9130 which functions or subroutines. */
9133 gfc_resolve (gfc_namespace *ns)
9135 gfc_namespace *old_ns;
9137 old_ns = gfc_current_ns;
9142 gfc_current_ns = old_ns;